TITLE: Phylogenetic Analysis of Complete Bovine Coronavirus Genome Sequences
AUTHORS: Aspen M Workman, Tara   G. McDaneld, Gregory P P Harhay, Subha Das, John   Dustin Loy, Benjamin   M. Hause

[DESCRIPTION]
Bovine coronavirus (BCoV) has spilled over to many species, including humans, where the host range variant coronavirus OC43 is endemic. Balance of the opposing activities of the surface spike (S) and hemagglutinin esterase (HE) glycoproteins control virion avidity which is critical for interspecies transmission and host adaptation.  Here, 78 genomes were sequenced directly from clinical samples collected between 2013 and 2022 from cattle in 12 states, primarily in the Midwestern U.S.  Relatively little genetic diversity was observed, with genomes having >98% nucleotide identity. Eleven isolates collected between 2020 and 2022 from four states (Nebraska, Colorado, California, and Wisconsin) contained a 12-nucleotide insertion in the receptor-binding domain (RBD) of the HE gene identical to one recently reported in China, and a single genome from Nebraska collected in 2020 contained a novel 12-nucleotide deletion in the HE gene RBD. Isogenic HE proteins containing either the insertion or deletion in the HE RBD maintained esterase activity and the ability to bind bovine submaxillary mucin, a substrate enriched in the receptor 9-O-acetylated-sialic acid, despite modeling that predicted structural changes in the HE R3 loop critical for receptor binding.  The emergence of BCoV with structural variants in the RBD raises the possibility of further interspecies transmission.

[STEPS]
SECTION: MAFFT Alignment
3. Align  using MAFFT v7.450

The sequences were aligned using the MAFFT accuracy  methods   ---globalpair as defined in the mafft manual page  below 

 


Keep in mind ...

The commands below generate files that contain both STDOUT (mafft output to the "terminal") as well as the multifasta alignment file (.afa). The STDOUT is found at the beginning of the file and the alignments follow the STDOUT
 All .afa_out files were split in a text editor into two separate .out and .afa files for downstream analysis.  
The input FASTA file headers contain description information; this header information should be stripped out to facilitate readable, compact leaf names/identifiers. Use the following command to strip out the description text.  
 thethe


SECTION: Tree Building
4.  RaXML Tee Building


Notes  below are cut from the RaXML manual to put the command line parameters used in context
 


raxmlHPC reports 192_Genomes_Final_6_14_22.jtt_100.globalpair.SEQID_Only.afa has  3438 DNA alignment patterns 

 

Given the directions above from the RaXML manual, use 8 threads for raxmlHPC tree building

 end

Output :

    
  
      

Use RAxML_bipartitionsBranchLabels.192_Genomes_Final_6_14_22.jtt_100.globalpair.SEQID_Only.T8.N2000 file to create phylogenetic tree with bootstrap support values provided as branch labels.


Visualize this tree of 192 Bovine coronavirus genomes, MAFFT Global Pair Alignment, PAM 100 in the Interactive Tree of Life 


    
The pdf suitable for download is    



SECTION: MAFFT Alignment
3.2. Align 192 genomes  using the default PAM (JTT) 200 substitution matrix 
 
Outputs :    


This  was split into    &  

 
Use sed command to strip out the description from the header in FASTA


 Create SEQID_Only file   




SECTION: MAFFT Alignment
3.4. Align SPIKE using the default PAM (JTT) 200 substitution matrix 
 
Outputs :    


This  was split into   &   


Use sed command to strip out the description from the header in FASTA  

 Create SEQID_Only file   

SECTION: MAFFT Alignment
3.6. Align HE using the default PAM (JTT) 200 substitution matrix 
 
Outputs :    

This  was split into    &   


Use sed command to strip out the description from the header in FASTA

Create SEQID_Only file   

SECTION: Software
1. 
 

SECTION: MAFFT Alignment
3.3. Align 192 BCV genomes  using the  PAM (JTT) 100 substitution matrix 
 
Outputs :    


This  was split into    &  

 
Use sed command to strip out the description from the header in FASTA


 Create SEQID_Only file   

SECTION: MAFFT Alignment
3.5. Align SPIKE  using the  PAM (JTT) 100 substitution matrix 
 
Outputs :    


This  was split into  &   


Use sed command to strip out the description from the header in FASTA

 Create SEQID_Only file   

SECTION: MAFFT Alignment
3.7. Align HE using the default PAM (JTT) 100 substitution matrix 
 
Outputs :    


This  was split into:    &   

Use sed command to strip out the description from the header in FASTA


 Create SEQID_Only file  

SECTION: MAFFT Alignment
3.1. The input FASTA files are:

       

SECTION: Software
2.