Comprehensive Workflow for the Genome-wide Identification and Expression Meta-analysis of the ATL E3 Ubiquitin Ligase Gene Family in Grapevine

Podsumowanie

This article describes the procedure for the identification and characterization of a gene family in grapevine applied to the family of Arabidopsis Tóxicos in Levadura (ATL) E3 ubiquitin ligases.

Streszczenie

Classification and nomenclature of genes in a family can significantly contribute to the description of the diversity of encoded proteins and to the prediction of family functions based on several features, such as the presence of sequence motifs or of particular sites for post-translational modification and the expression profile of family members in different conditions. This work describes a detailed protocol for gene family characterization. Here, the procedure is applied to the characterization of the Arabidopsis Tóxicos in Levadura (ATL) E3 ubiquitin ligase family in grapevine. The methods include the genome-wide identification of family members, the characterization of gene localization, structure, and duplication, the analysis of conserved protein motifs, the prediction of protein localization and phosphorylation sites as well as gene expression profiling across the family in different datasets. Such procedure, which could be extended to further analyses depending on experimental purposes, could be applied to any gene family in any plant species for which genomic data are available, and it provides valuable information to identify interesting candidates for functional studies, giving insights into the molecular mechanisms of plant adaptation to their environment.

Wprowadzenie

During the last decade, much research has been carried out in grapevine genomics. Grapevine is a recognized economically relevant crop, which has become a model for research on fruit development and on the responses of woody plants to biotic and abiotic stresses. In this context, the release of the Vitis vinifera cv. PN40024 genome in 2007¹ and its updated version in 2011² led to a rapid accumulation of "Omics"-scale data and to a burst of high-throughput studies. Based on the published sequence data, the comprehensive analysis of a given gene family (generally composed of proteins sharing conserved motifs, structural and/or functional similarities and evolutionary relationships), can now be performed to uncover its molecular functions, evolution, and gene expression profiles. These analyses can contribute to understanding how gene families control physiological processes at a genome-wide level.

Many aspects of the plant life cycle are regulated by ubiquitin-mediated degradation of key proteins, which require a fine-tuned turnover to ensure regular cellular processes. Important components of the ubiquitin-mediated degradation process are the E3 ubiquitin ligases, which are responsible for system flexibility, thanks to the recruitment of specific targets³. Accordingly, these enzymes represent a huge gene family, with around 1,400 E3 ligase-encoding genes predicted in Arabidopsis thaliana genome⁴, each E3 ubiquitin ligase acting for the ubiquitination of specific target proteins. Despite the importance of substrate-specific ubiquitination in cellular regulation in plants, little is known about how the ubiquitination pathway is regulated and target proteins have been identified only in a few cases. The deciphering of such specificity and regulation mechanisms relies first on the identification and characterization of the different components of the system, in particular the E3 ligases. Among ubiquitin ligases, the ATL subfamily is characterized by 91 members identified in A. thaliana displaying a RING-H2 finger domain^5,6, some of them playing a role in defense and hormone responses⁷.

The first crucial step to define the members of a new gene family is the precise definition of the family features, such as consensus motifs, key domains, and protein sequence characteristics. Indeed, the reliable retrieval of all gene family members based on BLAST analysis requires some mandatory sequence characteristics, in particular protein domains responsible for protein function/activity, serving as protein signature. This can be facilitated by previous characterization of the same gene family in other plant species or achieved by analyzing different genes putatively belonging to the same family in different plant species, to isolate common sequences. The family members can then be individually named following common rules settled by international consortia for a given plant species. In grapevine, for instance, such procedure is subjected to the recommendations of the Super-Nomenclature Committee for Grape Gene Annotation (sNCGGa), establishing the construction of a phylogenetic tree including V. vinifera and A. thaliana gene family members to allow gene annotation based on nucleotide sequences⁸.

Chromosome localization of family members and gene duplication survey allow highlighting the presence of whole-genome or tandem duplicated genes. Such information appears useful to unravel putative gene functions, since it might show functional redundancy or reveal different situations, i.e., non-functionalization, neo-functionalization, or sub-functionalization⁹. Both neo- and sub-functionalization are important events that create genetic novelty, providing new cellular components for plant adaptation to changing environments¹⁰. In particular, duplications of ancestral genes and production of new genes were very frequent during the evolution of the grapevine genome and newly formed genes originating from proximal and tandem duplications in grapevine were more likely to produce new functions¹¹.

Another key factor in deciphering gene family function is the transcriptomic profile. The availability of public databases giving access to a huge amount of transcriptomic data can be thus exploited to assign putative functions to gene family members using large-scale in silico expression analyses. Indeed, the peculiar expression of some genes in specific plant organs or in response to certain stresses can give some hints regarding the putative roles of the corresponding proteins in defined conditions, and give support to hypotheses about possible sub-functionalization of duplicated genes to respond to different challenges. For that purpose, it is important to consider several datasets: these can be already available gene expression matrixes, such as the genome-wide transcriptomic atlas of grapevine organs and developmental stages¹², or can be built ad hoc by retrieving transcriptomic datasets for the particular plant species subjected to defined stresses. Moreover, a simple approach using two matrices, one with pairwise similarity data and the other one with pairwise co-expression coefficients can be applied to evaluate the relationships between sequence similarity and expression patterns within a gene family.

The aim of this work is to provide a global approach, defining gene structure, conserved protein motifs, chromosomal location, gene duplications, and expression patterns, as well the prediction of protein localization and phosphorylation sites, to attain an exhaustive characterization of a gene family in plants. Such a comprehensive approach is applied here to the characterization of the ATL E3 ubiquitin ligase family in grapevine. According to the emerging role of ATL subfamily members in regulating key cellular processes⁷, this work can well assist the identification of strong candidates for functional studies, and eventually unravel the molecular mechanisms governing the adaptation of this important crop to its environment.

Protokół

1. Identification of Putative ATL Gene Family Member(s)

1. PSI-BLAST web version
   1. Open the BLAST web page¹³ and click on the protein BLAST section.
   2. In the "Enter Query sequence" field, enter the amino acid sequence of the protein (here VIT_05s0077g01970) that will be used as the probe to identify the other family members.
      NOTE: A good representative protein should be used (a protein displaying all the important features that characterize the family).
   3. In the field "Choose search set", select the "Reference protein" database (refseq_protein) and the organism of interest (V. vinifera - taxid:29760).
   4. In the field "Program selection", select PSI-BLAST algorithm and click the BLAST button to run the analysis.
      NOTE: By clicking on the "Algorithm parameters" it is possible to adjust some advanced parameters (Max target sequences, Scoring matrix, PSI-BLAST threshold, etc.).
   5. The first BLAST round retrieves all the sequences displaying relevant matches with the query (e-value above the selected threshold - by default 0.005; 0.001 in this experiment). Unselect all the entries, which clearly do not belong to the family under examination by clicking on the tick in the "select for PSI-BLAST" column and run the second PSI-BLAST iteration by clicking the BLAST button as in step 1.1.4.
   6. Newly identified sequences are highlighted in yellow. Unselect the clearly wrong retrieved hits and uncover further iterations as described in step 1.1.5.
   7. Continue with iterations until the algorithm does not find any relevant entry or it reaches convergence (no new entries are found). Download the list of putative gene family members for further analyses. Visually inspect the retrieved hits in each iteration to avoid the presence of false positives.
2. PSI-BLAST standalone version
   1. Download the standalone version of BLAST by clicking the "download BLAST" button on the BLAST home page¹³.
      NOTE: The standalone BLAST software is a command line version of the web interface described before. It enables executing the PSI-BLAST search against a custom local or remote database. Moreover, it allows searching with a pre-defined Position Specific Score Matrix (PSSM).

2. Manual Inspection of the PSI-BLAST-identified Family Members

1. Multiple alignment
   1. Collect the amino acidic sequences previously identified in a FASTA-formatted file and upload it into the MEGA software¹⁴ to proceed with the multiple alignment.
   2. Open the MEGA software, click the "Align" button, click "Edit/Build Alignment", click "Create a new alignment", click "Protein".
   3. Click "Edit" from the alignment menu and "Insert Sequence from File". Browse for the FASTA file created before and confirm the upload of all the surveyed sequences.
   4. Click "Alignment" from the alignment menu and "Align by MUSCLE". Use default parameters, click "Compute" button, and wait for the completion of the multiple alignment.
   5. Visually inspect the multiple alignment to exclude incorrectly predicted family members. The canonical CxxC(13x)PxCxHxxHxxCxxxW(7x)CxxCW motif, (in particular the presence of the proline residue before the third cysteine), is the key feature required to define the ATL family members.
2. Analysis of specific LOGO
   1. Submit the definitive list of family members (96 grapevine sequences fulfill the requirements to be considered ATL) to the Multiple Em for Motif Elicitation (MEME)¹⁵ to define conserved motifs across the family.
   2. From the MEME home page, click the "MEME" button, and complete the "Data Submission Form" with particular information regarding the family of interest.
   3. Use MEME analysis to confirm the presence of the two expected motifs within the grapevine ATL family members, i.e., the RING-H2 and the GLD motifs.
3. Alternatively, perform steps 2.1 and 2.2 simultaneously using the bioinformatics software suite (see Table of Materials).
   1. Upload FASTA file (see step 2.1.1) into the suite. Select "File" from the menu, then "Import" and click "From file". Browse the FASTA file and click "Open".
   2. Select all the imported sequences in the list and click on "Align/Assemble" button in the toolbar, then click "Pairwise Multiple Alignment". Select "Muscle alignment" and click "OK" to launch the alignment using default parameters.
   3. To visualize the LOGO of the alignment, click on "Graphs" → "options" and select "Sequence Logo".

3. Analysis of Protein Physical Parameters and Domains

1. As the definition of the different physical parameters of the surveyed family members is important to have a comprehensive description of the family, submit the list of family members to specific web tools.
   1. For isoelectric point (pI) and molecular weight (kDa), use the ProtParam tool¹⁶ on the Expasy website with default parameters.
   2. For protein subcellular localization, use different tools to obtain a more reliable prediction such as ngLOC v1.0¹⁷ with default settings, targetP v1.1¹⁸ with default settings, and protein prowler subcellular localization v1.2¹⁹ with a cut-off of probability of 0.5. For phosphorylation sites, use the MUsite v1.0 web tool²⁰ with default parameters.
2. Investigate additional protein domains in family members.
   1. Open the Pfam database webpage²¹, select "Sequence search" tool, submit protein sequences in the query box, and click "Go" to run the analysis.
      NOTE: Each protein sequence is analyzed individually. An e-value of 1.0 in the default setting allows discriminating between significant and non-significant hits.
   2. Open the TMHMM Server²² from the Center for Biological Sequence Analysis to investigate the presence of putative transmembrane regions. Paste all protein sequences simultaneously in the query box (or alternatively upload a text file including all protein sequences in FASTA format) and click "Submit" to run the analysis.
   3. Analyze proteins lacking predicted transmembrane domains, according to TMHMM (step 3.2.2), with ProtScale tool to identify putative hydrophobic regions. Open ProtScale webpage²³. Paste each protein sequence in the query box and select "Hphob. / Kyte & Doolittle" as amino acid scale. Click "Submit" to run the analysis.

4. Chromosomal Distribution, Duplications, and Exon-intron Organization

1. Map the ATL family members on the chromosomes based on the information retrieved from the Grapevine Genome CRIBI Biotech Center website²⁴.
   1. Browse the PhenoGram website homepage²⁵. Write the "Input File" as a tab-delimited text file with the specific features of the genes to be mapped on the chromosomes, according to the exhaustive guidelines and examples regarding the compilation of the provided file following the path "Phenogram" → "Documentation" → "Options" → "Input file".
   2. Write the "Title" of the work. Select the genome to be drawn. For genomes not implemented in the software, such as the grapevine genome, select "other" in the drop-down menu. Write the genome file according to the guidelines and examples provided, following the path "Phenogram" → "Documentation" → "Options" → "Genome", and upload it.
   3. Use default parameters of "Phenotype spacing", "Phenotype color", "Image format", or select alternatives in the respective menus, and click "Plot" to obtain the visualization of the genes on the chromosomes.
2. Define the duplication state of the family members using the MCScanX software²⁶.
   1. Download and unzip a copy of MCscanX on a local machine running command lines 1 (Supplementary File 1). Enter the MCscanX folder and create the required executables running command lines 2 (Supplementary File 1).
      NOTE: Installation of MCscanX is known to fail on some Linux 64 bit machines due to an issue regarding the function chdir. If an error message is returned related to this function upon the make command execution, the command lines 3 (Supplementary File 1) should be run and the command "make" should be attempted afterwards.
   2. Download the V. vinifera proteins and the annotation file running command lines 4 (Supplementary File 1).
      NOTE: The grapevine annotation file needs to be unzipped and the single chromosomes information cat in a unique file by running command lines 5 (Supplementary File 1).
   3. Run an "all versus all" blastp search using the V. vinifera protein file as both the query and the subject.
   4. Create a searchable blast database using the V. vinifera protein file running command lines 6 (Supplementary File 1). Perform the blastp search by using the V. vinifera proteins file as a query against the database created previously by running command lines 7 (Supplementary File 1).
   5. Convert the annotation file in a suitable format for MCScanX. Run command lines 8 (Supplementary File 1) to download the custom perl script parseMSCanXgff.pl. Perform the analysis running command lines 9 (Supplementary File 1).
      NOTE: A file vitis.gff is generated that holds gene coordinates in the following format:
      sp# gene starting position ending position
      where "sp" is a two-letter code for the species (Vv for grapevine) whereas "#" is the name of the scaffold. Note that the provided custom perl script is suitable for most conversion, although some code modification may be required in some specific cases due to the diversity of the information provided in the available annotation file.
   6. Launch MCScanX running command lines 10 (Supplementary File 1).
      NOTE: The "vitis" is the prefix of both the annotation and the blast output file. This represents a compulsory requirement for the software to run.
   7. Analyze MCScanX results. MCScanX produces one text file "vitis.collinearity", which contains collinear blocks. Such a file can be inspected by any text editor (see example output 1 Supplementary File 1).
      NOTE: A "mcscaxOutput.html" directory is generated that contains html files featuring multiple alignments of collinear blocks against each reference chromosome. These files can be inspected through a web browser.
   8. Classify paralogous genes based on their relative positions in chromosomes running command lines 11 (Supplementary File 1).
      NOTE: Paralogous gene classification is described in Supplementary Table II. The generated output file "vitis.gene_type" contains all origin information with a simple tab delimited format.
   9. Perform enrichment analysis to evaluate whether the gene family has prevalently originated by a specific mechanism running command lines 12 (Supplementary File 1).
      NOTE: File "vitis.gene_type" is generated at step 4.2.8, whereas file "gene_family_file" represents a one line text file in which the name of the family (e.g., ATL_genes) is followed by the locus names for the all the genes belonging to the family separated by a tab. The applied statistical test for enrichment is a Fisher exact test and the p-values of different origins are stored in the file "outputFile.txt".
3. Visualize the exon-intron organization of the genes using Interactive Tree Of Life (iTOL)²⁷, an on-line tool for the display, annotation, and management of phylogenetic trees.
   1. Upload a phylogenetic tree in the "Upload" section of the iTOL website. The tree is built according to Section 5 below. For each family member gene, retrieve gene structure prediction from the V1 annotation of the grapevine genome (CRIBI website cited above). Calculate the length (in bp) of putative exons, introns, and untranslated regions (UTRs).
   2. Use the "Protein domains" dataset for graphical visualization of the exon-intron pattern. Write a plain text file including calculated lengths according to the specifications provided following the path "Help" → "Help pages" → "Dataset types" → "Protein domains" in the iTOL website²⁷. Using "Protein domains" dataset, the "rectangle (RE)" and the "rectangle gap (GP)" shapes represent the exon and the UTRs, respectively.

5. Phylogenetic Analysis and Nomenclature

1. Analyze the relationships among ATL family members through the construction of a high quality phylogenetic tree and the definition of a family nomenclature.
   1. For a grapevine gene family, follow the rules established by the Grapevine Super Nomenclature Committee⁸.
   2. Retrieve A. thaliana ATL sequences, required as reference for grapevine gene nomenclature⁸, from the UniProt database²⁸ .
   3. Write a FASTA file including all nucleotide sequences of grapevine and A. thaliana gene family members to be included in the phylogenetic analysis. The nucleotide sequences allow the maximum of variability among family members (compared to protein sequences).
2. Phylogenetic tree
   NOTE: The use of the Phylogeny.fr ²⁹ pipeline is recommended to get a high quality phylogenetic tree, but not mandatory.
   1. Browse the Phylogeny.fr homepage²⁹, and select the "Phylogeny analysis" pipeline.
      NOTE: "One Click" is suitable in most of the cases, but if needed it is possible to select specific advanced settings ("Advanced") or even a fully customized analysis ("A la Carte"; see step 5.2.5).
   2. Write the "Name of the analysis", upload the FASTA file created previously (step 5.2.1, and click "Submit" to run the analysis.
   3. Alternatively, if the procedure described above (steps 5.2.1, 5.2.2) results in an error message, complete each step of the Phylogeny suite pipeline individually, as follows.
      1. From the MUSCLE software homepage³⁰, upload the FASTA file in "STEP 1", select "Pearson/FASTA" as "Output format" in "STEP 2", and click "Submit" in "STEP 3" to align query sequences.
      2. Click "Download alignment file" and save as FASTA file for further steps.
      3. Process the alignment FASTA file to eliminate poorly aligned positions using Gblocks Server tool³¹. Upload the alignment FASTA file, select "DNA" as "Type of sequence" and chose the option(s) of stringency that best fits with the analysis (e.g., for grapevine ATL gene family select all the three options proposed for "less stringent selection" because of high sequence divergence). Click "Get blocks" to run the analysis.
      4. Click "Resulting alignment" at the bottom of the output page and save the results as a new FASTA file.
      5. From the Phylogeny.fr homepage²⁹, select "A la Carte" as "Phylogeny analysis" pipeline. Then, deselect "Multiple alignment" and "Alignment curation". Click "Create workflow", upload the Gblocks-curated FASTA file (step 5.2.5.4), select "Bootstrapping procedure" with default parameters in "Settings", and click "Submit" to run the analysis.
   4. Collapse poorly supported branches (i.e., bootstrap values < 70%) by clicking "Collapse branches" in the "Select and action" section and download the final results in the Newick format to further analyses.
3. Assign a gene name based on the phylogeny.
   1. Review the phylogenetic tree to evaluate the reliability of the tree structure by uploading it into the iTOL suite cited above (section 4.3).
   2. Assign manually a gene name to each family member. In the case of one-to-one orthologues, assign the Arabidopsis-like name (e.g., AtATL3 → VviATL3). Differentiate grapevine genes (two or more) deriving from a single Arabidopsis homolog with the same phylogenetic distance using numbers, or letters if the Arabidopsis gene ends with a number (e.g., AtATL23 → VviATL23a, VviATL23b).
   3. In the case of one-to-many or many-to-many orthologues, assign a new gene name composed of the Arabidopsis-like name (here, "ATL") paired with a number higher than the highest number already used for both V. vinifera and Arabidopsis (e.g., VviATL83).
   4. Complete the nomenclature of the newly defined family descending from the top to the bottom of the phylogenetic tree.

6. Grapevine Organ and Stage Expression Profiling

1. Generate the working data matrix containing expression data for the family members.
   1. Download the V. vinifera cv. Corvina gene expression Atlas datamatrix from the link distributed on the ResearchGate platform³². This file contains the RMA normalized expression values to be used in following steps.
   2. Extract the expression values for each family gene from the Atlas datamatrix and write a "working datamatrix" containing the same header row as the Atlas datamatrix. Save the "working datamatrix" as a tab-delimited text file.
2. Perform the hierarchical bi-clustered analysis using Multi Experiment Viewer (MeV) software.
   1. Download and install MeV software³³.
   2. Upload the "working datamatrix" (step 6.1.2) following the path "File" → "Load Data" → "Browse" and select the text file. Select "Single-color Array" and remove the tick from "Load Annotation" when an automatic annotation is not provided. Select the upper-leftmost expression value of the expression table preview and click the "Load" button.
   3. Adjust the data applying Log2 transformation ("Adjust Data" → "Log Transformations" → "Log2 Transform") and Gene/Row normalization ("Adjust Data" → "Gene/Row Adjustments" → "Median Center Gene/Row"). Set the proper scale limit ("Display" → "Set Color Scale Limits").
   4. Calculate the Hierarchical Clustering following the path "Analysis" → "Clustering" → "HCL". Select "Optimize Gene Leaf Order" and "Optimize Sample Leaf Order" in "Ordering Optimization field", "Pearson Correlation" in the "Distance Matrix Selection" field, and "Average linkage clustering" in the "Linkage Method Selection" field. Then, click "OK" to run the analysis.
   5. View the results in the "Analysis Results" → "HCL" menu on the left panel of the window. Export the heat map by clicking "Save Image" in the "File" menu.

7. Expression Profiling in Response to Biotic and Abiotic Stresses

1. Repeat step 6.1 with the GSE accession ID obtained from respective publications and studies investigating biotic and abiotic stress on grapevine. For example, experiments providing the transcriptome profile of grapevine berries infected with the fungal pathogen Botrytis cinerea using the NimbleGen Grape Whole-genome microarray can be browsed with GSE ID of GSE52586. Repeat steps 6.1.1 and 6.1.2.
2. Search the NCBI Sequence Reads Archive³⁴ with the SRA/BioProject ID (e.g., SRP055458 or PRJNA275778 for "grapevine flower shading" experiments) and download all associated raw sequence reads. RNA-seq datasets from many different studies are processed using a single pipeline for consistency.
   1. Briefly, trim raw sequence FASTQ reads (single- and pair-end) and filter quality with Trimmomatic³⁵. Use an AVGQUAL and MINLEN filter of 20 and 40, respectively and all parameters default.
   2. Index the 12X grapevine reference genome¹ using Bowtie2³⁶. Download the 12X grapevine reference genome (e.g., bowtie2-build) before running bowtie2 command.
   3. Obtain count matrix tables with htseq-count³⁷ using the grapevine V1 gene model annotation (GFF/GTF) file.
3. Perform differential gene expression (re-)analysis in R³⁸ with limma³⁹ libraries for RMA-normalized matrices and DESeq2⁴⁰ libraries for count matrix tables obtained from steps 7.1.1 and 7.2.1, respectively.
   1. Perform a standard "two-group" comparison (i.e., "treatment"/"control"). Ensure that the design matrix/groupings of "controls" and "treatment" conditions are properly specified.
      NOTE: A typical design for microarray differential expression analysis (GSE52586) to compare EL-33 berries infected with Botrytis cinerea against control (healthy) berries at the same development stage with limma running command lines 13 is shown in Supplementary File 1. A typical design for RNA-seq differential expression analysis (SRP055458 or PRJNA275778) to compare flower (at 7 days after cap-fall) under shade treatment against the control with DESeq2 running command lines 14 is shown in Supplementary File 1.
   2. Obtain the lists of differentially expressed genes (DEG) in each contrast, for limma, use the functions lmFit(), followed by eBayes(), and then by topTable()functions, while for DESeq2, use the DESeqDataSetFromMatrix(), DESeq(), and results() functions. Below, a typical workflow to be followed.
      1. For microarray differential expression analysis, see command lines 15 (Supplementary File 1). For RNA-seq differential expression analysis see command lines 16 (Supplementary File 1). Repeat the above steps for all other contrasts with different appropriate design schema (See examples in step 7.3.1)
4. From the lists of DEGs generated, extract all rows that do not correspond to ATL V1 accession, retain columns containing the log2 Fold Change (Treatment/Control) > |0.5| and adjusted p-values (FDR) < 0.05, and merge them accordingly into a matrix table, whether a study falls into "abiotic" or "biotic/pathogen interaction" compendia.
5. Construct the hierarchical clustered heatmaps (abiotic and biotic compendia) in R using the libraries gplots.
   NOTE: Calling the heatmap.2 function constructs the heatmap along with row dendrograms from the respective matrix tables. Additional arguments using cellnote function helps to distinguish differentially expressed (log2FC > 0.5, FDR < 0.05) ATL genes in each comparison across a large range of experimental conditions by a * symbol. Apply the typical workflow in R running command lines 17 (Supplementary File 1) or alternatively, repeat steps 6.2.2 to 6.2.5 to construct the heatmaps using MeV software.

8. Analysis of the Relationships Between Paralogous Sequence Divergence and Gene Co-expression

1. Construct the matrix containing pairwise similarity. The elements of the similarity matrix are the values of sequence similarity calculated from the pairwise protein alignments.
   1. Use the EMBOSS needle web server⁴¹ with default settings to make pairwise sequence alignments and save as text file. Open the output text file and remove all comment lines, together with column and row names to generate a file called "similarityTable.txt".
      NOTE: Such a table features a line for each ATL gene reporting the similarity values calculated in each of the pairwise alignment. The order of the loci in rows and columns is the same so that a symmetric matrix is generated with respect of the diagonal values.
2. Construct the matrix with co-expression data by calculating the Pearson correlation coefficient. The following procedure requires R and the perl module PDL.
   1. Download the expression values for the 96 ATL genes running command lines 18 (Supplementary File 1) within a terminal. Perform a co-expression analysis by using a custom perl script that can be downloaded by running command lines 19 (Supplementary File 1). Such script will calculate the Pearson correlation coefficient between pairs of ATL loci as previously reported.
   2. Launch the script running command lines 20 (Supplementary File 1) and follow the output instructions. The script will produce an output file (namely "coexpressionTable.txt") containing a co-expression matrix featuring the same locus names order of matrix obtained in step 8.1 (this ordering is essential to run the Mantel test, see below).
3. Perform a Mantel test between the data matrices obtained in steps 8.1 and 8.2. After entering the R environment (run command "R" from within a terminal), load the ade4 library using the following command: library(ade4)
   1. Run the Mantel test by loading the two data matrices and performing the statistics running command lines 21 (Supplementary File 1), with "nrep" representing the number of permutations. The test consists of calculating the correlation between the elements of these matrices, permuting the matrices and then calculating the same test statistic again.
      NOTE: All the obtained values of the statistic test are used to build a reference distribution of the statistic test, which will be used to calculate a p-value to test for significance. The number of permutations defines the precision with which the p-value can be obtained.

Wyniki

The VIT_05s0077g01970 gene, identified as the most similar to A. thaliana ATL2 (At3g16720) through a BLASTp search, was used as probe to survey the ATL family members in the grapevine genome (V. vinifera cv Pinot Noir PN40024). The PSI-BLAST analysis converged after a few cycles revealing a list of putative genes belonging to the grapevine ATL gene family (Figure 1A). The presence of the canonical RING-H2 domain for each candidate was evalua...

Dyskusje

In the genomic era, many gene families have been deeply characterized in several plant species. This information is preliminary to functional studies and provide a frame to investigate further the role of different members in a family. In this context, there is also a need for a nomenclature system allowing to uniquely identify each member in a family, avoiding the redundancy and confusions that may arise when names are assigned independently to different genes by different research groups.

Af...

Materiały

Name	Company	Catalog Number	Comments
Personal computer
Basic Local Alignment Search Tool (BLAST)			https://blast.ncbi.nlm.nih.gov/Blast.cgi
Molecular Evolutionary Genetics Analysis (MEGA)			http://www.megasoftware.net/
Motif-based sequence analysis tools (MEME)			http://meme-suite.org/
Geneious	Biomatters Limited		http://www.geneious.com/
ProtParam Tool			http://web.expasy.org/protparam/
ngLOC			http://genome.unmc.edu/ngLOC/index.html
TargetP v1.1 Server			http://www.cbs.dtu.dk/services/TargetP/
Protein Prowler			http://bioinf.scmb.uq.edu.au:8080/pprowler_webapp_1-2/
MUsite			http://musite.sourceforge.net/
Pfam			http://pfam.xfam.org/
TMHMM Server v. 2.0			http://www.cbs.dtu.dk/services/TMHMM/
ProtScale			http://web.expasy.org/protscale/
Grape Genome Database (CRIBI)			http://genomes.cribi.unipd.it/grape/
PhenoGram			http://visualization.ritchielab.psu.edu/phenograms/plot
MCScanX			http://chibba.pgml.uga.edu/mcscan2/
Interactive Tree Of Life (iTOL)			http://itol.embl.de/
UniProt			http://www.uniprot.org/
Phylogeny.fr			http://www.phylogeny.fr/index.cgi
MUSCLE			http://www.ebi.ac.uk/Tools/msa/muscle/
Gblocks Server			http://molevol.cmima.csic.es/castresana/Gblocks_server.html
Vitis vinifera cv. Corvina gene expression Atlas datamatrix			https://www.researchgate.net/publication/273383414_54sample_ datamatrix_geneIDs_Fasoli2012
Multi Experiment Viewer (MeV)			http://mev.tm4.org/#/welcome
Sequence Read Archive (SRA)			https://www.ncbi.nlm.nih.gov/sra
R			https://www.r-project.org/
EMBOSS Needle (EMBL-EBI)			http://www.ebi.ac.uk/Tools/psa/emboss_needle/