Plant Reactome Knowledgebase: empowering plant pathway exploration and OMICS data analysis

Plant Reactome (https://plantreactome.gramene.org) is a freely accessible, comprehensive plant pathway knowledgebase. It provides curated reference pathways from rice (Oryza sativa) and gene-orthology-based pathway projections to 129 additional species, spanning single-cell photoautotrophs, non-vascular plants, and higher plants, thus encompassing a wide-ranging taxonomic diversity. Currently, Plant Reactome houses a collection of 339 reference pathways, covering metabolic and transport pathways, hormone signaling, genetic regulations of developmental processes, and intricate transcriptional networks that orchestrate a plant's response to abiotic and biotic stimuli. Beyond being a mere repository, Plant Reactome serves as a dynamic data discovery platform. Users can analyze and visualize omics data, such as gene expression, gene-gene interaction, proteome, and metabolome data, all within the rich context of plant pathways. Plant Reactome is dedicated to fostering data interoperability, upholding global data standards, and embracing the tenets of the Findable, Accessible, Interoperable and Re-usable (FAIR) data policy.

Reference genome of the nutrition-rich orphan crop chia (Salvia hispanica) and its implications for future breeding

Chia (Salvia hispanica L.) is one of the most popular nutrition-rich foods and pseudocereal crops of the family Lamiaceae. Chia seeds are a rich source of proteins, polyunsaturated fatty acids (PUFAs), dietary fibers, and antioxidants. In this study, we present the assembly of the chia reference genome, which spans 303.6 Mb and encodes 48,090 annotated protein-coding genes. Our analysis revealed that ~42% of the chia genome harbors repetitive content, and identified ~3 million single nucleotide polymorphisms (SNPs) and 15,380 simple sequence repeat (SSR) marker sites. By investigating the chia transcriptome, we discovered that ~44% of the genes undergo alternative splicing with a higher frequency of intron retention events. Additionally, we identified chia genes associated with important nutrient content and quality traits, such as the biosynthesis of PUFAs and seed mucilage fiber (dietary fiber) polysaccharides. Notably, this is the first report of in-silico annotation of a plant genome for protein-derived small bioactive peptides (biopeptides) associated with improving human health. To facilitate further research and translational applications of this valuable orphan crop, we have developed the Salvia genomics database (SalviaGDB), accessible at https://salviagdb.org.

Genotype and phenotype data standardization, utilization and integration in the big data era for agricultural sciences

Large-scale genotype and phenotype data have been increasingly generated to identify genetic markers, understand gene function and evolution and facilitate genomic selection. These datasets hold immense value for both current and future studies, as they are vital for crop breeding, yield improvement and overall agricultural sustainability. However, integrating these datasets from heterogeneous sources presents significant challenges and hinders their effective utilization. We established the Genotype-Phenotype Working Group in November 2021 as a part of the AgBioData Consortium (https://www.agbiodata.org) to review current data types and resources that support archiving, analysis and visualization of genotype and phenotype data to understand the needs and challenges of the plant genomic research community. For 2021-22, we identified different types of datasets and examined metadata annotations related to experimental design/methods/sample collection, etc. Furthermore, we thoroughly reviewed publicly funded repositories for raw and processed data as well as secondary databases and knowledgebases that enable the integration of heterogeneous data in the context of the genome browser, pathway networks and tissue-specific gene expression. Based on our survey, we recommend a need for (i) additional infrastructural support for archiving many new data types, (ii) development of community standards for data annotation and formatting, (iii) resources for biocuration and (iv) analysis and visualization tools to connect genotype data with phenotype data to enhance knowledge synthesis and to foster translational research. Although this paper only covers the data and resources relevant to the plant research community, we expect that similar issues and needs are shared by researchers working on animals. Database URL: https://www.agbiodata.org.

Exploring Pan-Genomes: An Overview of Resources and Tools for Unraveling Structure, Function, and Evolution of Crop Genes and Genomes

The availability of multiple sequenced genomes from a single species made it possible to explore intra- and inter-specific genomic comparisons at higher resolution and build clade-specific pan-genomes of several crops. The pan-genomes of crops constructed from various cultivars, accessions, landraces, and wild ancestral species represent a compendium of genes and structural variations and allow researchers to search for the novel genes and alleles that were inadvertently lost in domesticated crops during the historical process of crop domestication or in the process of extensive plant breeding. Fortunately, many valuable genes and alleles associated with desirable traits like disease resistance, abiotic stress tolerance, plant architecture, and nutrition qualities exist in landraces, ancestral species, and crop wild relatives. The novel genes from the wild ancestors and landraces can be introduced back to high-yielding varieties of modern crops by implementing classical plant breeding, genomic selection, and transgenic/gene editing approaches. Thus, pan-genomic represents a great leap in plant research and offers new avenues for targeted breeding to mitigate the impact of global climate change. Here, we summarize the tools used for pan-genome assembly and annotations, web-portals hosting plant pan-genomes, etc. Furthermore, we highlight a few discoveries made in crops using the pan-genomic approach and future potential of this emerging field of study.

Biocuration of a Transcription Factors Network Involved in Submergence Tolerance during Seed Germination and Coleoptile Elongation in Rice (Oryza sativa)

Modeling biological processes and genetic-regulatory networks using in silico approaches provides a valuable framework for understanding how genes and associated allelic and genotypic differences result in specific traits. Submergence tolerance is a significant agronomic trait in rice; however, the gene-gene interactions linked with this polygenic trait remain largely unknown. In this study, we constructed a network of 57 transcription factors involved in seed germination and coleoptile elongation under submergence. The gene-gene interactions were based on the co-expression profiles of genes and the presence of transcription factor binding sites in the promoter region of target genes. We also incorporated published experimental evidence, wherever available, to support gene-gene, gene-protein, and protein-protein interactions. The co-expression data were obtained by re-analyzing publicly available transcriptome data from rice. Notably, this network includes OSH1, OSH15, OSH71, Sub1B, ERFs, WRKYs, NACs, ZFP36, TCPs, etc., which play key regulatory roles in seed germination, coleoptile elongation and submergence response, and mediate gravitropic signaling by regulating OsLAZY1 and/or IL2. The network of transcription factors was manually biocurated and submitted to the Plant Reactome Knowledgebase to make it publicly accessible. We expect this work will facilitate the re-analysis/re-use of OMICs data and aid genomics research to accelerate crop improvement.

The Gene Ontology knowledgebase in 2023

The Gene Ontology (GO) knowledgebase (http://geneontology.org) is a comprehensive resource concerning the functions of genes and gene products (proteins and noncoding RNAs). GO annotations cover genes from organisms across the tree of life as well as viruses, though most gene function knowledge currently derives from experiments carried out in a relatively small number of model organisms. Here, we provide an updated overview of the GO knowledgebase, as well as the efforts of the broad, international consortium of scientists that develops, maintains, and updates the GO knowledgebase. The GO knowledgebase consists of three components: (1) the GO-a computational knowledge structure describing the functional characteristics of genes; (2) GO annotations-evidence-supported statements asserting that a specific gene product has a particular functional characteristic; and (3) GO Causal Activity Models (GO-CAMs)-mechanistic models of molecular "pathways" (GO biological processes) created by linking multiple GO annotations using defined relations. Each of these components is continually expanded, revised, and updated in response to newly published discoveries and receives extensive QA checks, reviews, and user feedback. For each of these components, we provide a description of the current contents, recent developments to keep the knowledgebase up to date with new discoveries, and guidance on how users can best make use of the data that we provide. We conclude with future directions for the project.

Plant Reactome and PubChem: The Plant Pathway and (Bio)Chemical Entity Knowledgebases

Plant Reactome (https://plantreactome.gramene.org) and PubChem ( https://pubchem.ncbi.nlm.nih.gov ) are two reference data portals and resources for curated plant pathways, small molecules, metabolites, gene products, and macromolecular interactions. Plant Reactome knowledgebase, a conceptual plant pathway network, is built by biocuration and integrating (bio)chemical entities, gene products, and macromolecular interactions. It provides manually curated pathways for the reference species Oryza sativa (rice) and gene orthology-based projections that extend pathway knowledge to 106 plant species. Currently, it hosts 320 reference pathways for plant metabolism, hormone signaling, transport, genetic regulation, plant organ development and differentiation, and biotic and abiotic stress responses. In addition to the pathway browsing and search functions, the Plant Reactome provides the analysis tools for pathway comparison between reference and projected species, pathway enrichment in gene expression data, and overlay of gene-gene interaction data on pathways. PubChem, a popular reference database of (bio)chemical entities, provides information on small molecules and other types of chemical entities, such as siRNAs, miRNAs, lipids, carbohydrates, and chemically modified nucleotides. The data in PubChem is collected from hundreds of data sources, including Plant Reactome. This chapter provides a brief overview of the Plant Reactome and the PubChem knowledgebases, their association to other public resources providing accessory information, and how users can readily access the contents.

Ensembl Genomes 2022: an expanding genome resource for non-vertebrates

Ensembl Genomes (https://www.ensemblgenomes.org) provides access to non-vertebrate genomes and analysis complementing vertebrate resources developed by the Ensembl project (https://www.ensembl.org). The two resources collectively present genome annotation through a consistent set of interfaces spanning the tree of life presenting genome sequence, annotation, variation, transcriptomic data and comparative analysis. Here, we present our largest increase in plant, metazoan and fungal genomes since the project's inception creating one of the world's most comprehensive genomic resources and describe our efforts to reduce genome redundancy in our Bacteria portal. We detail our new efforts in gene annotation, our emerging support for pangenome analysis, our efforts to accelerate data dissemination through the Ensembl Rapid Release resource and our new AlphaFold visualization. Finally, we present details of our future plans including updates on our integration with Ensembl, and how we plan to improve our support for the microbial research community. Software and data are made available without restriction via our website, online tools platform and programmatic interfaces (available under an Apache 2.0 license). Data updates are synchronised with Ensembl's release cycle.

Plant lectins and their many roles: Carbohydrate-binding and beyond

Lectins are ubiquitous proteins that reversibly bind to specific carbohydrates and, thus, serve as readers of the sugar code. In photosynthetic organisms, lectin family proteins play important roles in capturing and releasing photosynthates via an endogenous lectin cycle. Often, lectin proteins consist of one or more lectin domains in combination with other types of domains. This structural diversity of lectins is the basis for their current classification, which is consistent with their diverse functions in cell signaling associated with growth and development, as well as in the plant's response to biotic, symbiotic, and abiotic stimuli. Furthermore, the lectin family shows evolutionary expansion that has distinct clade-specific signatures. Although the function(s) of many plant lectin family genes are unknown, studies in the model plant Arabidopsis thaliana have provided insights into their diverse roles. Here, we have used a biocuration approach rooted in the critical review of scientific literature and information available in the public genomic databases to summarize the expression, localization, and known functions of lectins in Arabidopsis. A better understanding of the structure and function of lectins is expected to aid in improving agricultural productivity through the manipulation of candidate genes for breeding climate-resilient crops, or by regulating metabolic pathways by applications of plant growth regulators.

Chia (Salvia hispanica) Gene Expression Atlas Elucidates Dynamic Spatio-Temporal Changes Associated With Plant Growth and Development

Chia (Salvia hispanica L.), now a popular superfood and a pseudocereal, is one of the richest sources of dietary nutrients such as protein, fiber, and polyunsaturated fatty acids (PUFAs). At present, the genomic and genetic information available in the public domain for this crop are scanty, which hinders an understanding of its growth and development and genetic improvement. We report an RNA-sequencing (RNA-Seq)-based comprehensive transcriptome atlas of Chia sampled from 13 tissue types covering vegetative and reproductive growth stages. We used ~355 million high-quality reads of total ~394 million raw reads from transcriptome sequencing to generate de novo reference transcriptome assembly and the tissue-specific transcript assemblies. After the quality assessment of the merged assemblies and implementing redundancy reduction methods, 82,663 reference transcripts were identified. About 65,587 of 82,663 transcripts were translated into 99,307 peptides, and we were successful in assigning InterPro annotations to 45,209 peptides and gene ontology (GO) terms to 32,638 peptides. The assembled transcriptome is estimated to have the complete sequence information for ~86% of the genes found in the Chia genome. Furthermore, the analysis of 53,200 differentially expressed transcripts (DETs) revealed their distinct expression patterns in Chia's vegetative and reproductive tissues; tissue-specific networks and developmental stage-specific networks of transcription factors (TFs); and the regulation of the expression of enzyme-coding genes associated with important metabolic pathways. In addition, we identified 2,411 simple sequence repeats (SSRs) as potential genetic markers from the transcripts. Overall, this study provides a comprehensive transcriptome atlas, and SSRs, contributing to building essential genomic resources to support basic research, genome annotation, functional genomics, and molecular breeding of Chia.

Beyond gene ontology (GO): using biocuration approach to improve the gene nomenclature and functional annotation of rice S-domain kinase subfamily

The S-domain subfamily of receptor-like kinases (SDRLKs) in plants is poorly characterized. Most members of this subfamily are currently assigned gene function based on the S-locus Receptor Kinase from Brassica that acts as the female determinant of self-incompatibility (SI). However, Brassica like SI mechanisms does not exist in most plants. Thus, automated Gene Ontology (GO) pipelines are not sufficient for functional annotation of SDRLK subfamily members and lead to erroneous association with the GO biological process of SI. Here, we show that manual bio-curation can help to correct and improve the gene annotations and association with relevant biological processes. Using publicly available genomic and transcriptome datasets, we conducted a detailed analysis of the expansion of the rice (Oryza sativa) SDRLK subfamily, the structure of individual genes and proteins, and their expression.The 144-member SDRLK family in rice consists of 82 receptor-like kinases (RLKs) (67 full-length, 15 truncated),12 receptor-like proteins, 14 SD kinases, 26 kinase-like and 10 GnK2 domain-containing kinases and RLKs. Except for nine genes, all other SDRLK family members are transcribed in rice, but they vary in their tissue-specific and stress-response expression profiles. Furthermore, 98 genes show differential expression under biotic stress and 98 genes show differential expression under abiotic stress conditions, but share 81 genes in common.Our analysis led to the identification of candidate genes likely to play important roles in plant development, pathogen resistance, and abiotic stress tolerance. We propose a nomenclature for 144 SDRLK gene family members based on gene/protein conserved structural features, gene expression profiles, and literature review. Our biocuration approach, rooted in the principles of findability, accessibility, interoperability and reusability, sets forth an example of how manual annotation of large-gene families can fill in the knowledge gap that exists due to the implementation of automated GO projections, thereby helping to improve the quality and contents of public databases.

Gramene 2021: harnessing the power of comparative genomics and pathways for plant research

Gramene (http://www.gramene.org), a knowledgebase founded on comparative functional analyses of genomic and pathway data for model plants and major crops, supports agricultural researchers worldwide. The resource is committed to open access and reproducible science based on the FAIR data principles. Since the last NAR update, we made nine releases; doubled the genome portal's content; expanded curated genes, pathways and expression sets; and implemented the Domain Informational Vocabulary Extraction (DIVE) algorithm for extracting gene function information from publications. The current release, #63 (October 2020), hosts 93 reference genomes-over 3.9 million genes in 122 947 families with orthologous and paralogous classifications. Plant Reactome portrays pathway networks using a combination of manual biocuration in rice (320 reference pathways) and orthology-based projections to 106 species. The Reactome platform facilitates comparison between reference and projected pathways, gene expression analyses and overlays of gene-gene interactions. Gramene integrates ontology-based protein structure-function annotation; information on genetic, epigenetic, expression, and phenotypic diversity; and gene functional annotations extracted from plant-focused journals using DIVE. We train plant researchers in biocuration of genes and pathways; host curated maize gene structures as tracks in the maize genome browser; and integrate curated rice genes and pathways in the Plant Reactome.

Gramene 2021: harnessing the power of comparative genomics and pathways for plant research

Gramene (http://www.gramene.org), a knowledgebase founded on comparative functional analyses of genomic and pathway data for model plants and major crops, supports agricultural researchers worldwide. The resource is committed to open access and reproducible science based on the FAIR data principles. Since the last NAR update, we made nine releases; doubled the genome portal's content; expanded curated genes, pathways and expression sets; and implemented the Domain Informational Vocabulary Extraction (DIVE) algorithm for extracting gene function information from publications. The current release, #63 (October 2020), hosts 93 reference genomes-over 3.9 million genes in 122 947 families with orthologous and paralogous classifications. Plant Reactome portrays pathway networks using a combination of manual biocuration in rice (320 reference pathways) and orthology-based projections to 106 species. The Reactome platform facilitates comparison between reference and projected pathways, gene expression analyses and overlays of gene-gene interactions. Gramene integrates ontology-based protein structure-function annotation; information on genetic, epigenetic, expression, and phenotypic diversity; and gene functional annotations extracted from plant-focused journals using DIVE. We train plant researchers in biocuration of genes and pathways; host curated maize gene structures as tracks in the maize genome browser; and integrate curated rice genes and pathways in the Plant Reactome.

Ensembl Genomes 2020-enabling non-vertebrate genomic research

Ensembl Genomes (http://www.ensemblgenomes.org) is an integrating resource for genome-scale data from non-vertebrate species, complementing the resources for vertebrate genomics developed in the context of the Ensembl project (http://www.ensembl.org). Together, the two resources provide a consistent set of interfaces to genomic data across the tree of life, including reference genome sequence, gene models, transcriptional data, genetic variation and comparative analysis. Data may be accessed via our website, online tools platform and programmatic interfaces, with updates made four times per year (in synchrony with Ensembl). Here, we provide an overview of Ensembl Genomes, with a focus on recent developments. These include the continued growth, more robust and reproducible sets of orthologues and paralogues, and enriched views of gene expression and gene function in plants. Finally, we report on our continued deeper integration with the Ensembl project, which forms a key part of our future strategy for dealing with the increasing quantity of available genome-scale data across the tree of life.

Plant Reactome: a knowledgebase and resource for comparative pathway analysis

Plant Reactome (https://plantreactome.gramene.org) is an open-source, comparative plant pathway knowledgebase of the Gramene project. It uses Oryza sativa (rice) as a reference species for manual curation of pathways and extends pathway knowledge to another 82 plant species via gene-orthology projection using the Reactome data model and framework. It currently hosts 298 reference pathways, including metabolic and transport pathways, transcriptional networks, hormone signaling pathways, and plant developmental processes. In addition to browsing plant pathways, users can upload and analyze their omics data, such as the gene-expression data, and overlay curated or experimental gene-gene interaction data to extend pathway knowledge. The curation team actively engages researchers and students on gene and pathway curation by offering workshops and online tutorials. The Plant Reactome supports, implements and collaborates with the wider community to make data and tools related to genes, genomes, and pathways Findable, Accessible, Interoperable and Re-usable (FAIR).

Gramene 2018: unifying comparative genomics and pathway resources for plant research

Gramene (http://www.gramene.org) is a knowledgebase for comparative functional analysis in major crops and model plant species. The current release, #54, includes over 1.7 million genes from 44 reference genomes, most of which were organized into 62,367 gene families through orthologous and paralogous gene classification, whole-genome alignments, and synteny. Additional gene annotations include ontology-based protein structure and function; genetic, epigenetic, and phenotypic diversity; and pathway associations. Gramene's Plant Reactome provides a knowledgebase of cellular-level plant pathway networks. Specifically, it uses curated rice reference pathways to derive pathway projections for an additional 66 species based on gene orthology, and facilitates display of gene expression, gene-gene interactions, and user-defined omics data in the context of these pathways. As a community portal, Gramene integrates best-of-class software and infrastructure components including the Ensembl genome browser, Reactome pathway browser, and Expression Atlas widgets, and undergoes periodic data and software upgrades. Via powerful, intuitive search interfaces, users can easily query across various portals and interactively analyze search results by clicking on diverse features such as genomic context, highly augmented gene trees, gene expression anatomograms, associated pathways, and external informatics resources. All data in Gramene are accessible through both visual and programmatic interfaces.

Involving community in genes and pathway curation

Biocuration plays a crucial role in building databases and complex systems-level platforms required for processing, annotating and analyzing ‘Big Data’ in biology. However, biocuration efforts cannot keep pace with a dramatic increase in the production of omics data; this presents one of the bottlenecks in genomics. In two pathway curation jamborees, Plant Reactome curators tested strategies for introducing researchers to pathway curation tools, harnessing biologists’ expertise in curating plant pathways and developing a network of community biocurators. We summarize the strategy, workflow and outcomes of these exercises, and discuss the role of community biocuration in advancing databases and genomic resources.

AgBioData consortium recommendations for sustainable genomics and genetics databases for agriculture

The future of agricultural research depends on data. The sheer volume of agricultural biological data being produced today makes excellent data management essential. Governmental agencies, publishers and science funders require data management plans for publicly funded research. Furthermore, the value of data increases exponentially when they are properly stored, described, integrated and shared, so that they can be easily utilized in future analyses. AgBioData (https://www.agbiodata.org) is a consortium of people working at agricultural biological databases, data archives and knowledgbases who strive to identify common issues in database development, curation and management, with the goal of creating database products that are more Findable, Accessible, Interoperable and Reusable. We strive to promote authentic, detailed, accurate and explicit communication between all parties involved in scientific data. As a step toward this goal, we present the current state of biocuration, ontologies, metadata and persistence, database platforms, programmatic (machine) access to data, communication and sustainability with regard to data curation. Each section describes challenges and opportunities for these topics, along with recommendations and best practices.

Variant Effect Prediction Analysis Using Resources Available at Gramene Database

The goal of Gramene database ( www.gramene.org ) is to empower the plant research community in conducting comparative genomics studies across model plants and crops by employing a phylogenetic framework and orthology-based projections. Gramene database (release #49) provides resources for comparative plant genomics including well-annotated plant genomes (39 complete reference genomes and six partial genomes), genetic or structural variation data for 14 plant species, pathways for 58 plant species, and gene expression data for 14 species including Arabidopsis, rice, maize, soybean, wheat, etc. (fetched from EBI-EMBL Gene Expression Atlas database). Gramene also facilitates visualization and analysis of user-defined data in the context of species-specific Genome Browsers or pathways. This chapter describes basic navigation for Gramene users and illustrates how they can use the genome section to analyze the gene expression and nucleotide variation data generated in their labs. This includes (1) upload and display of genomic data onto a Genome Browser track, (2) analysis of variation data using online Variant Effect Predictor (VEP) tool for smaller data sets, and (3) the use of the stand-alone Perl scripts and command line protocols for variant effect prediction on larger data sets.

Towards an open grapevine information system

Viticulture, like other fields of agriculture, is currently facing important challenges that will be addressed only through sustained, dedicated and coordinated research. Although the methods used in biology have evolved tremendously in recent years and now involve the routine production of large data sets of varied nature, in many domains of study, including grapevine research, there is a need to improve the findability, accessibility, interoperability and reusability (FAIR-ness) of these data. Considering the heterogeneous nature of the data produced, the transnational nature of the scientific community and the experience gained elsewhere, we have formed an open working group, in the framework of the International Grapevine Genome Program (www.vitaceae.org), to construct a coordinated federation of information systems holding grapevine data distributed around the world, providing an integrated set of interfaces supporting advanced data modeling, rich semantic integration and the next generation of data mining tools. To achieve this goal, it will be critical to develop, implement and adopt appropriate standards for data annotation and formatting. The development of this system, the GrapeIS, linking genotypes to phenotypes, and scientific research to agronomical and oeneological data, should provide new insights into grape biology, and allow the development of new varieties to meet the challenges of biotic and abiotic stress, environmental change, and consumer demand.

Gramene Database: Navigating Plant Comparative Genomics Resources

Gramene (http://www.gramene.org) is an online, open source, curated resource for plant comparative genomics and pathway analysis designed to support researchers working in plant genomics, breeding, evolutionary biology, system biology, and metabolic engineering. It exploits phylogenetic relationships to enrich the annotation of genomic data and provides tools to perform powerful comparative analyses across a wide spectrum of plant species. It consists of an integrated portal for querying, visualizing and analyzing data for 44 plant reference genomes, genetic variation data sets for 12 species, expression data for 16 species, curated rice pathways and orthology-based pathway projections for 66 plant species including various crops. Here we briefly describe the functions and uses of the Gramene database.

Plant Reactome: a resource for plant pathways and comparative analysis

Plant Reactome (http://plantreactome.gramene.org/) is a free, open-source, curated plant pathway database portal, provided as part of the Gramene project. The database provides intuitive bioinformatics tools for the visualization, analysis and interpretation of pathway knowledge to support genome annotation, genome analysis, modeling, systems biology, basic research and education. Plant Reactome employs the structural framework of a plant cell to show metabolic, transport, genetic, developmental and signaling pathways. We manually curate molecular details of pathways in these domains for reference species Oryza sativa (rice) supported by published literature and annotation of well-characterized genes. Two hundred twenty-two rice pathways, 1025 reactions associated with 1173 proteins, 907 small molecules and 256 literature references have been curated to date. These reference annotations were used to project pathways for 62 model, crop and evolutionarily significant plant species based on gene homology. Database users can search and browse various components of the database, visualize curated baseline expression of pathway-associated genes provided by the Expression Atlas and upload and analyze their Omics datasets. The database also offers data access via Application Programming Interfaces (APIs) and in various standardized pathway formats, such as SBML and BioPAX.

FragariaCyc: A Metabolic Pathway Database for Woodland Strawberry Fragaria vesca

FragariaCyc is a strawberry-specific cellular metabolic network based on the annotated genome sequence of Fragaria vesca L. ssp. vesca, accession Hawaii 4. It was built on the Pathway-Tools platform using MetaCyc as the reference. The experimental evidences from published literature were used for supporting/editing existing entities and for the addition of new pathways, enzymes, reactions, compounds, and small molecules in the database. To date, FragariaCyc comprises 66 super-pathways, 488 unique pathways, 2348 metabolic reactions, 3507 enzymes, and 2134 compounds. In addition to searching and browsing FragariaCyc, researchers can compare pathways across various plant metabolic networks and analyze their data using Omics Viewer tool. We view FragariaCyc as a resource for the community of researchers working with strawberry and related fruit crops. It can help understanding the regulation of overall metabolism of strawberry plant during development and in response to diseases and abiotic stresses. FragariaCyc is available online at http://pathways.cgrb.oregonstate.edu.

Gramene 2016: comparative plant genomics and pathway resources

Gramene (http://www.gramene.org) is an online resource for comparative functional genomics in crops and model plant species. Its two main frameworks are genomes (collaboration with Ensembl Plants) and pathways (The Plant Reactome and archival BioCyc databases). Since our last NAR update, the database website adopted a new Drupal management platform. The genomes section features 39 fully assembled reference genomes that are integrated using ontology-based annotation and comparative analyses, and accessed through both visual and programmatic interfaces. Additional community data, such as genetic variation, expression and methylation, are also mapped for a subset of genomes. The Plant Reactome pathway portal (http://plantreactome.gramene.org) provides a reference resource for analyzing plant metabolic and regulatory pathways. In addition to ∼ 200 curated rice reference pathways, the portal hosts gene homology-based pathway projections for 33 plant species. Both the genome and pathway browsers interface with the EMBL-EBI's Expression Atlas to enable the projection of baseline and differential expression data from curated expression studies in plants. Gramene's archive website (http://archive.gramene.org) continues to provide previously reported resources on comparative maps, markers and QTL. To further aid our users, we have also introduced a live monthly educational webinar series and a Gramene YouTube channel carrying video tutorials.

The floral transcriptome of Eucalyptus grandis

As a step toward functional annotation of genes required for floral initiation and development within the Eucalyptus genome, we used short read sequencing to analyze transcriptomes of floral buds from early and late developmental stages, and compared these with transcriptomes of diverse vegetative tissues, including leaves, roots, and stems. A subset of 4807 genes (13% of protein-coding genes) were differentially expressed between floral buds of either stage and vegetative tissues. A similar proportion of genes were differentially expressed among all tissues. A total of 479 genes were differentially expressed between early and late stages of floral development. Gene function enrichment identified 158 gene ontology classes that were overrepresented in floral tissues, including 'pollen development' and 'aromatic compound biosynthetic process'. At least 40 floral-dominant genes lacked functional annotations and thus may be novel floral transcripts. We analyzed several genes and gene families in depth, including 49 putative biomarkers of floral development, the MADS-box transcription factors, 'S-domain'-receptor-like kinases, and selected gene family members with phosphatidylethanolamine-binding protein domains. Expanded MADS-box gene subfamilies in Eucalyptus grandis included SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1), SEPALLATA (SEP) and SHORT VEGETATIVE PHASE (SVP) Arabidopsis thaliana homologs. These data provide a rich resource for functional and evolutionary analysis of genes controlling eucalypt floral development, and new tools for breeding and biotechnology.

De novo transcriptome assembly and analyses of gene expression during photomorphogenesis in diploid wheat Triticum monococcum

BACKGROUND

Triticum monococcum (2n) is a close ancestor of T. urartu, the A-genome progenitor of cultivated hexaploid wheat, and is therefore a useful model for the study of components regulating photomorphogenesis in diploid wheat. In order to develop genetic and genomic resources for such a study, we constructed genome-wide transcriptomes of two Triticum monococcum subspecies, the wild winter wheat T. monococcum ssp. aegilopoides (accession G3116) and the domesticated spring wheat T. monococcum ssp. monococcum (accession DV92) by generating de novo assemblies of RNA-Seq data derived from both etiolated and green seedlings.

PRINCIPAL FINDINGS

The de novo transcriptome assemblies of DV92 and G3116 represent 120,911 and 117,969 transcripts, respectively. We successfully mapped ∼90% of these transcripts from each accession to barley and ∼95% of the transcripts to T. urartu genomes. However, only ∼77% transcripts mapped to the annotated barley genes and ∼85% transcripts mapped to the annotated T. urartu genes. Differential gene expression analyses revealed 22% more light up-regulated and 35% more light down-regulated transcripts in the G3116 transcriptome compared to DV92. The DV92 and G3116 mRNA sequence reads aligned against the reference barley genome led to the identification of ∼500,000 single nucleotide polymorphism (SNP) and ∼22,000 simple sequence repeat (SSR) sites.

CONCLUSIONS

De novo transcriptome assemblies of two accessions of the diploid wheat T. monococcum provide new empirical transcriptome references for improving Triticeae genome annotations, and insights into transcriptional programming during photomorphogenesis. The SNP and SSR sites identified in our analysis provide additional resources for the development of molecular markers.

VitisCyc: a metabolic pathway knowledgebase for grapevine (Vitis vinifera)

We have developed VitisCyc, a grapevine-specific metabolic pathway database that allows researchers to (i) search and browse the database for its various components such as metabolic pathways, reactions, compounds, genes and proteins, (ii) compare grapevine metabolic networks with other publicly available plant metabolic networks, and (iii) upload, visualize and analyze high-throughput data such as transcriptomes, proteomes, metabolomes etc. using OMICs-Viewer tool. VitisCyc is based on the genome sequence of the nearly homozygous genotype PN40024 of Vitis vinifera "Pinot Noir" cultivar with 12X v1 annotations and was built on BioCyc platform using Pathway Tools software and MetaCyc reference database. Furthermore, VitisCyc was enriched for plant-specific pathways and grape-specific metabolites, reactions and pathways. Currently VitisCyc harbors 68 super pathways, 362 biosynthesis pathways, 118 catabolic pathways, 5 detoxification pathways, 36 energy related pathways and 6 transport pathways, 10,908 enzymes, 2912 enzymatic reactions, 31 transport reactions and 2024 compounds. VitisCyc, as a community resource, can aid in the discovery of candidate genes and pathways that are regulated during plant growth and development, and in response to biotic and abiotic stress signals generated from a plant's immediate environment. VitisCyc version 3.18 is available online at http://pathways.cgrb.oregonstate.edu.

Gramene 2013: comparative plant genomics resources

Gramene (http://www.gramene.org) is a curated online resource for comparative functional genomics in crops and model plant species, currently hosting 27 fully and 10 partially sequenced reference genomes in its build number 38. Its strength derives from the application of a phylogenetic framework for genome comparison and the use of ontologies to integrate structural and functional annotation data. Whole-genome alignments complemented by phylogenetic gene family trees help infer syntenic and orthologous relationships. Genetic variation data, sequences and genome mappings available for 10 species, including Arabidopsis, rice and maize, help infer putative variant effects on genes and transcripts. The pathways section also hosts 10 species-specific metabolic pathways databases developed in-house or by our collaborators using Pathway Tools software, which facilitates searches for pathway, reaction and metabolite annotations, and allows analyses of user-defined expression datasets. Recently, we released a Plant Reactome portal featuring 133 curated rice pathways. This portal will be expanded for Arabidopsis, maize and other plant species. We continue to provide genetic and QTL maps and marker datasets developed by crop researchers. The project provides a unique community platform to support scientific research in plant genomics including studies in evolution, genetics, plant breeding, molecular biology, biochemistry and systems biology.

Socio-economic and ethnic inequalities in diabetes retinal screening

OBJECTIVE

We aimed to quantify socio-economic and ethnic inequalities in diabetes retinal screening.

METHODS

Data were analysed for the retinal screening programme for three South London boroughs for the 18-month period to February 2009. Sight-threatening diabetic retinopathy (STDR) was defined as the occurrence of diabetic maculopathy, severe non-proliferative or proliferative diabetic retinopathy. Odds ratios were adjusted for sex, age group, duration and type of diabetes, self-reported ethnicity and deprivation quintile by participant postal code.

RESULTS

There were 76 351 records obtained but, after excluding duplicate and ineligible records, data were analysed for 59 495 records from 31 484 subjects. There were 7026 (22%) subjects called for appointments who were not screened in the period, with 24 458 (78%) having one or more screening episodes. Non-attendance for screening was highest in young adults aged 18-34 years (32%) and in those aged 85 years or greater (28%). In the most deprived quintile, non-attendance was 23% compared with 21% in the least deprived quintile [odds ratio (OR) 1.37, 95% confidence interval (CI) 1.16-1.61, P < 0.001]. There were 2819 (11.5%) participants with STDR, including 10.8% in the least deprived quintile and 12.2% in the most deprived quintile (OR 1.10, 95% CI 0.95-1.16, P = 0.196). Compared with white Europeans (9.4%), STDR was higher in Africans (15.2%) and African Caribbeans (14.7%), resulting from a higher frequency of diabetic maculopathy.

CONCLUSION

Socio-economic inequality in diabetes retinal screening may be smaller than reported in earlier studies. This study suggested an increased frequency of diabetic maculopathy among participants of African origins.

Structural modules for receptor dimerization in the S-locus receptor kinase extracellular domain

The highly polymorphic S-locus receptor kinase (SRK) is the stigma determinant of specificity in the self-incompatibility response of the Brassicaceae. SRK spans the plasma membrane of stigma epidermal cells, and it is activated in an allele-specific manner on binding of its extracellular region (eSRK) to its cognate pollen coat-localized S-locus cysteine-rich (SCR) ligand. SRK, like several other receptor kinases, forms dimers in the absence of ligand. To identify domains in SRK that mediate ligand-independent dimerization, we assayed eSRK for self-interaction in yeast. We show that SRK dimerization is mediated by two regions in eSRK, primarily by a C-terminal region inferred by homology modeling/fold recognition techniques to assume a PAN_APPLE-like structure, and secondarily by a region containing a signature sequence of the S-domain gene family, which might assume an EGF-like structure. We also show that eSRK exhibits a marked preference for homodimerization over heterodimerization with other eSRK variants and that this preference is mediated by a small, highly variable region within the PAN_APPLE domain. Thus, the extensive polymorphism exhibited by the eSRK not only determines differential affinity toward the SCR ligand, as has been assumed thus far, but also underlies a previously unrecognized allelic specificity in SRK dimerization. We propose that preference for SRK homodimerization explains the codominance exhibited by a majority of SRKs in the typically heterozygous stigmas of self-incompatible plants, whereas an increased propensity for heterodimerization combined with reduced affinity of heterodimers for cognate SCRs might underlie the dominant-recessive or mutual weakening relationships exhibited by some SRK allelic pairs.

The menD and menE homologs code for 2-succinyl-6-hydroxyl-2,4-cyclohexadiene-1-carboxylate synthase and O-succinylbenzoic acid-CoA synthase in the phylloquinone biosynthetic pathway of Synechocystis sp. PCC 6803

The genome of the cyanobacterium Synechocystis sp. PCC 6803 contains genes identified as menD and menE, homologs of Escherichia coli genes that code for 2-succinyl-6-hydroxyl-2,4-cyclohexadiene-1-carboxylate (SHCHC) synthase and O-succinylbenzoic acid-CoA ligase in the menaquinone biosynthetic pathway. In cyanobacteria, the product of this pathway is 2-methyl-3-phytyl-1,4-naphthoquinone (phylloquinone), a molecule used exclusively as an electron transfer cofactor in Photosystem (PS) I. The menD(-) and menE(-) strains were generated, and both were found to lack phylloquinone. Hence, no alternative pathways exist in cyanobacteria to produce O-succinylbenzoyl-CoA. Q-band EPR studies of photoaccumulated quinone anion radical and optical kinetic studies of the P700(+) [F(A)/F(B)](-) backreaction indicate that in the mutant strains, plastoquinone-9 functions as the electron transfer cofactor in the A(1) site of PS I. At a light intensity of 40 microE m(-2) s(-1), the menD(-) and menE(-) mutant strains grew photoautotrophically and photoheterotrophically, but with doubling times slower than the wild type. Both of which are sensitive to high light intensities. Low-temperature fluorescence studies show that in the menD(-) and menE(-) mutants, the ratio of PS I to PS II is reduced relative to the wild type. Whole-chain electron transfer rates in the menD(-) and menE(-) mutant cells are correspondingly higher on a chlorophyll basis. The slower growth rate and high-light sensitivity of the menD(-) and menE(-) mutants are therefore attributed to a lower content of PS I per cell.

Interactions among COX1, COX2, and COX3 mRNA-specific translational activator proteins on the inner surface of the mitochondrial inner membrane of Saccharomyces cerevisiae

The core of the cytochrome c oxidase complex is composed of its three largest subunits, Cox1p, Cox2p, and Cox3p, which are encoded in mitochondrial DNA of Saccharomyces cerevisiae and inserted into the inner membrane from the inside. Mitochondrial translation of the COX1, COX2, and COX3 mRNAs is activated mRNA specifically by the nuclearly coded proteins Pet309p, Pet111p, and the concerted action of Pet54p, Pet122p, and Pet494p, respectively. Because the translational activators recognize sites in the 5'-untranslated leaders of these mRNAs and because untranslated mRNA sequences contain information for targeting their protein products, the activators are likely to play a role in localizing translation. Herein, we report physical associations among the mRNA-specific translational activator proteins, located on the matrix side of the inner membrane. These interactions, detected by coimmune precipitation and by two-hybrid experiments, suggest that the translational activator proteins could be organized on the surface of the inner membrane such that synthesis of Cox1p, Cox2p, and Cox3p would be colocalized in a way that facilitates assembly of the core of the cytochrome c oxidase complex. In addition, we found interactions between Nam1p/Mtf2p and the translational activators, suggesting an organized delivery of mitochondrial mRNAs to the translation system.

Diurnal regulation of plastid genes in Populus deltoides

Light regulates leaf and chloroplast development, together with overall chloroplast gene expression at various levels. Plants respond to diurnal and seasonal changes in light by changing expression of photosynthesis genes and metabolism. In Populus deltoides, a deciduous tree species, leaf development begins in the month of March and leaf maturation is attained by summer, which is subsequently followed by autumnal senescence and fall. In the present study, diurnal changes in the steady state transcript levels of plastid genes were examined in the fully developed leaves during summer season. Our results show that steady state level of the psaA/B, psbA, psbEFLJ and petA transcripts showed differential accumulation during diurnal cycle in summer. However, there was no significant change in the pigment composition during the day/night cycle. Our studies suggest that the diurnal regulation of steady state mRNA accumulation may play a crucial role during daily adjustments in plants life with rapidly changing light irradiance and temperature.

Targeted inactivation of the psaK1, psaK2 and psaM genes encoding subunits of Photosystem I in the cyanobacterium Synechocystis sp. PCC 6803

PsaK and PsaM are small, integral membrane proteins, which are associated with the Photosystem I complexes of cyanobacteria. The complete genome sequence of Synechocystis sp. PCC 6803 has revealed the presence of two unlinked psaK genes: psaK1 (ssr0390) and psaK2 (sll0629). To investigate structural and functional roles of the PsaK1, PsaK2 and PsaM polypeptides in Synechocystis sp. PCC 6803, we generated targeted mutants that lack the functional psaK1, psaK2 or psaM genes. Inactivation of psaK1, psaK2 or psaM did not affect photoautotrophic growth, photosynthetic activity and accumulation of other subunits of the Photosystem I complex. The psaK1 (-), psaK2 (-) and psaK1 (-) psaK2 (-) mutants showed normal levels of Photosystem I trimers, whereas the lack of PsaM resulted in a 75% reduction in the recovery of trimers compared to the wild type. A 6.2 kDa polypeptide was observed in the Photosystem I preparations from the wild type, but not from the psaK2 (-) strain, suggesting the presence of PsaK2 in the Photosystem I complexes. Using reverse-transcription and polymerase chain reaction, we confirmed the expression of the psaK2 gene in Synechocystis sp. PCC 6803. To conclude, both psaK1 and psaK2 are expressed in Synechocystis sp. PCC 6803 and the absence of both proteins results in only a small reduction in Photosystem I electron transport. The PsaM subunit is required for the formation of stable Photosystem I trimers.

Subgroup-A Rous sarcoma virus-induced growth stimulation of chick embryos infected via the chorioallantoic membrane

Chicks that hatch from eggs containing group specific antigen (gs antigen) of lymphoid leukosis virus (LLV) subgroups, grow poorly. In our laboratory for more precise identification of LLV-of subgroup A (LLV-A) resistant and susceptible genotypes by progeny testing, the chorioallantoic membrane (CAM) assay in complemented by liver tumour (LT) assay, wherein Rous sarcoma virus (RSV) of subgroup A (homologous to LLV-A) was used. The present study was conducted in a light breed (White Leghorn) and also in a heavy breed (Rhode Island Red) to ascertain the effect of infection on embryonic growth by RSV subgroup A. Mean relative body weight (rbw) of infected LT negative chicks of either breed exceeded the control highly significantly (P < 0.01) by 2%. However, neither the dose of virus inoculated per embryo, nor egg size influenced the relative body weight of day old chicks (P > 0.05). No difference in relative body weight of LT positive and control chicks was observed.

The genetic control of susceptibility to subgroup D RNA tumour virus infection in commercial breeds of chickens

The genetic control of resistance to RSV (RAV-50) in chick embryos inoculated via the chorioallantoic membrane route was studied in three breeds of chickens. The relative susceptibility to subgroup D infection was 0.01, 0.02 and 0.003 in broilers, White Leghorns and Bantam x Australorp crosses, respectively. The distribution pattern of putative dr gene frequencies were 0.80, 0.68 and 0.90, respectively, in the same three breeds. Liver tumours resulting from the virus infection were observed in a very few birds about 28 days after hatching.

Characterization of the orf31-petG gene cluster from the plastid genome of Populus deltoides

The orf31-petG gene cluster is located approximately 1.2 kb away from the psbEFLJ operon in the chloroplast genome of Populus deltoides. The orf31 (ycf7) encodes an unidentified polypeptide while the petG gene encodes subunit V of an important component, cytochrome b6/f complex, involved in photosynthetic electron transport. We have determined the nucleotide sequence of the orf31-petG gene cluster from the plastid genome of a tree, Populus deltoides. Our sequence analysis suggests that these genes possess high homology with the published sequences of these genes from other plants. Northern analysis suggests development dependent transcription of the orf31-petG cluster in leaves.

Expression of tumour phenotypes, hatchability and liver tumour mortality in two routes of inoculation in fowl

A total of 350 and 200 eleven-day-olc embryos (pooled breeds) of twelve hatch replicates were inoculated with pseudotype of Bryan high titre, RSV(RAV-49) of subgroup C viz CAM (chorioallantoic membrane) and YS (yolk sac) route, respectively. An increase in hatchability (about 16%) and decrease in the incidence of CAM(+) [71%] and LT(+) [47%] phenotypes was noticed when inoculation was done via YS route as compared to the inoculation via CAM routes. A delay in LT(+) mortality was also recorded in YS route of infection. Chi-square analysis within a route basis indicated highly significant contingency (P < 0.01) in association of CAM infection phenotypes and LT incidence phenotypes for CAM route of infection in contrast to the YS route of infection.

Genetic control of cellular infection by subgroups A and C RNA tumour viruses in guinea fowl

An investigation was carried out in guinea fowl to determine their susceptibility to infection by Rous sarcoma viruses of subgroups A and C. A standard dose of each subgroup virus was inoculated into 14-day-old embryos via the chorioallantoic membrane (CAM). On the 10th day after inoculation, 50% of the embryonic chorioallantoic membranes were harvested to assess their infection status (CAM(+) or (-)), while the rest were allowed to hatch. The hatchabilities of the embryos inoculated with subgroups A and C were about 50% and 57%, respectively. The relative sensitivities of guinea fowl to infection by viruses of subgroups A and C were observed to be 0.220 and 0.003, respectively, as compared to chickens (1.00). Mortality due to subgroup A virus-induced liver tumours (LT) was 54% and four phenotypic subclasses, namely CAM(+) LT(+), CAM(+) LT(-), CAM(-) LT(+) and CAM(-) LT(-), were observed in guinea fowl as in chickens. However, a higher incidence (31%) of conversely associated phenotypes, i.e. CAM(+) LT(-) and CAM(-) LT(+), were observed in guinea fowl. Mortality caused by subgroup A virus-induced liver tumours was first observed in inoculated guinea fowl keets during the 3rd week after hatching, and 93% of the mortality occurred within 6 weeks. The peak mortality occurred in the 4th week after hatching. The target organs for transformation were considered to be the liver and spleen because of the equal incidence of tumours in these organs. Males and females were equally likely to die from liver tumours.(ABSTRACT TRUNCATED AT 250 WORDS)

Genetic susceptibility of indigenous chicks to subgroup A Rous sarcoma virus inoculated via the chorioallantoic membrane

An investigation was made using chicks of two Indian indigenous breeds of fowl, Kadaknath and Aseel, to ascertain genetic resistance to infection by Rous sarcoma virus of subgroup A. A standard inoculation dose of 0.2 ml virus containing 1000 pock forming units ml-1 was injected via the chorioallantoic membrane (CAM) into the 11-day-old embryos that were subsequently hatched. The sensitivity of the two indigenous breeds was compared with the highly susceptible exotic White Leghorn (WL) strain maintained in the laboratory. The Kadaknath breed was about three-fold and Assel, about six-fold less sensitive than the WL strain, indicating superiority of the indigenous breeds over the exotic breed of fowl. Most of the CAM-susceptible chicks died of liver tumour (LT) and most of the CAM-resistant chicks survived. However, conversely associated tumour phenotype subclass chicks, i.e. CAM-susceptible LT-negative chicks that survived and CAM-resistant LT-positive chicks that died, occurred consistently in the three breeds of fowl. Nevertheless, the overall survival potential of Kadaknath chicks measured up to 8 weeks post-hatching was greater than that of Aseel chicks. Neither transformation of embryonic tissue prior to hatching nor the visceral metastasis including liver conformed with the degree of CAM-infection as measured by number of pocks on CAMs.

Genetics of post-hatching survival potential of Australorp chicks infected as embryos by subgroup A Rous sarcoma virus: further support to 4-allele genetic model

Embryos (II day-old) of Australorp breed were inoculated via chorioallantoic membrane (CAM) with subgroup A Rous sarcoma virus, and hatched subsequently. The post-hatch survival period in chicks was recorded upto the last chick that died by virus-induced liver tumour, which had a range from 3 to 50 days with an average of 13 +/- 8.7 days. The survival potential of progency tested Australorp parents selected on the basis of negative CAM-infection and those selected on uninoculated embryos, differed significantly (P less than 0.01) while maintaining an inverse relationship between liver tumour mortality and degrees of infection of CAMs. The homozygous susceptibles lacking either ar1 or ar2 or both alleles of the tva (tumour virus a) locus died within 7 days of post-hatching, supporting thereby 4-allele genetic model of tva locus recently proposed for the control of LT- and CAM-infection phenotypes.

Induced liver tumour deaths by subgroup A Rous sarcoma virus in chicks inoculated via chorioallantoic membrane, a genetic marker

A study was made using two strains of light breed (White Leghorn strains, A and B) and four heavy beeds (Rhode Island Red, New Hampshire, Australorp, Columbian) to evaluate the breed difference in survival potential of chicks that were infected as 11-day-old embryos via chorioallantoic membranes (CAMs) with a subgroup A Rous sarcoma virus. Of the 1185 chicks hatched over multiple hatch-replicates, 845 chicks died rapidly of a fibrosarcomatous liver tumour (LT) with a peak mortality about 74% attained by the second week, post-hatch, in the heavy breeds and more than 90% by the second week in the light breed. The breeds did not differ in induced LT mortality when the chicks hatched from eggs that had at least 25 pock counts on CAMs, apparently genetically susceptible, i.e. 25 biologically active virus particles were enough to induce an unpreventable fatal LT. However, low pock-count on CAMs did not act as a pointer for predicting genetic resistance to infection because about 23% of chicks developed from eggs that had no pocks on CAMs, apparently genetically resistant, also died of LT, requiring further studies.

Cellular immunity in chicks expressing fibrosarcomatous liver tumor following embryo infection with subgroup A Rous sarcoma virus

Bryan standard strain of Rous sarcoma virus (BS-RSV) of subgroup A was inoculated into heavy and light breeds of chicken embryos via chorioallantoic membrane (CAM) to ascertain cell-mediated immune response, as measured by a leukocyte migration inhibition (LMI) test. Chicks hatched from eggs with pock-positive CAMs were more likely to develop liver tumors than those hatched from eggs with pock-negative CAMs. Chicks that developed tumors usually had a positive cell-mediated immune response, and those that were negative for liver tumor were negative, based on the LMI test.