OrthoMCL (http://orthomcl.org/orthomcl/) was used to construct orthologous gene families between P. somniferum (CHM and HN1) and 12 other plant species. MUSCLE44 was utilized to construct multiple sequence alignments of 35 single-copy orthologs among 14 species. RAxML software45 (version 7.2.3) was carried out to construct the maximum likelihood tree with the PROTGAMMAAUTO model by using the sequence alignments with A. trichopoda as an outgroup (Supplementary Fig. S9). The MCMCTree program of PAML (http://abacus.gene.ucl.ac.uk/software/paml.html) was applied to estimate divergence time using CDS alignments transformed from protein alignments. Five calibration values were chosen from the Time Tree website (http://www.timetree.org). Expansions and contractions of orthologous gene families were determined using CAFÉ 2.2 (Computational Analysis of gene Family Evolution)46 (Supplementary Fig. S11). For whole-genome duplication (WGD) analysis, the syntenic region between and within P. somniferum (Pso), M. cordata (Mco), and A. coerulea (Aco) was determined by MCscanX47 based on the all-to-all BLASTP results. The protein sequences of homologous gene pairs in the syntenic region were extracted and aligned using the MUSCLE program44. Subsequently, the protein sequence alignments were converted into CDS files, and 4DTv values were calculated based on the CDS alignments, accompanying the correction of the HKY model.

Note: The content above has been extracted from a research article, so it may not display correctly.

Please log in to submit your questions online.
Your question will be posted on the Bio-101 website. We will send your questions to the authors of this protocol and Bio-protocol community members who are experienced with this method. you will be informed using the email address associated with your Bio-protocol account.

We use cookies on this site to enhance your user experience. By using our website, you are agreeing to allow the storage of cookies on your computer.