GAL-Cluster Conservation

Overview

The GAL1, GAL7, and GAL10 genes encode enzymes in the Leloir pathway for galactose metabolism. This project systematically located these genes across yeast species, compared their protein-level conservation, and investigated how the surrounding chromosomal architecture has evolved: revealing a contrast between gene-level conservation and large-scale genomic rearrangement.

Genomic Locations

Gene coordinates were identified using NCBI Genome Annotation pages and confirmed with assembly-specific filtering.

All three genes co-localize to a single chromosome in both species, consistent with their role as a co-regulated cluster.

Protein Sequence Similarity (C. albicans vs. S. cerevisiae)

Pairwise alignments were performed with the UniProt Align Tool (Clustal Omega) and verified with EMBOSS Needle (BLOSUM62).

GAL7 was selected for a broader cross-Candida BLAST analysis. A UniProt BLAST search against Candida species returned 27 results with 13 highly significant hits (E-value = 0.0). Multiple-sequence alignment of these hits showed protein identity ranging from ~79% to 100%, indicating strong conservation of GAL7 function across the genus.

PIPELINE

    #!/bin/bash
    module load minimap2

    minimap2 -ax splice Unknown-Yeast-Genome.fasta GAL_mRNA.fasta > GAL_alignment.sam

    #!/bin/bash
    module load samtools 

    samtools view -O BAM -o gal_alignments.bam GAL_alignments.sam
    samtools sort -O BAM -o gal_alignments.sort.bam gal_alignments.bam
    samtools index -b gal_alignments.sort.bam
    samtools view gal_alignments.sort.bam | cut -f1,2,3,4

    #!/bin/bash
    module load samtools
    module load anaconda3/2023.09
    source activate mummer4

    samtools faidx cAlb_genome.fna "NC_032089.1" > CAlb_Chr1.fasta
    samtools faidx Unknown-Yeast-Genome.fasta "scaffold_1" > UnknownYeast_scaffold1.fasta
    nucmer -t 2 --mum -p chr1_vs_scaffold1 CAlb_Chr1.fasta UnknownYeast_scaffold1.fasta
    mummerplot -t png --filter -R CAlb_Chr1.fasta -Q UnknownYeast_scaffold1.fasta \
        -p chr1_vs_scaffold1_plot chr1_vs_scaffold1.delta

    #!/bin/bash
    module load anaconda3/2023.09
    source activate mummer4

    nucmer -t 4 --mum -p mum_alignment cAlb_genome.fna Unknown-Yeast-Genome.fasta
    mummerplot -t png --filter -Q Unknown-Yeast-Genome.fasta -R cAlb_genome.fna \
        -p genome_plot mum_alignment.delta

    #!/usr/bin/env python3
    import csv

    with open("out.coords") as infile, open("circa_links.csv", "w", newline="") as outfile:
        
        # skip 4 header lines
        for _ in range(4):
            infile.readline()
        
        writer = csv.writer(outfile)
        writer.writerow(["chromosome1", "start1", "end1", "chromosome2", "start2", "end2"])
        
        for line in infile:
            fields = line.strip().split("\t")
            if len(fields) < 9:
                continue
            s1, e1, s2, e2 = fields[0], fields[1], fields[2], fields[3]
            chr1, chr2 = fields[-2], fields[-1]
            writer.writerow([chr1, s1, e1, chr2, s2, e2])

    #!/bin/bash

    # Combine both into one file
    awk '{print $1","$2}' cAlb_genome.fna.fai > chromosome_lengths.csv
    awk '{print $1","$2}' Unknown-Yeast-Genome.fasta.fai >> chromosome_lengths.csv

    # Add header
    sed -i '1s/^/chromosome,length\n/' chromosome_lengths.csv

Key Findings:

Mapping to the unknown Candida draft genome confirmed that the GAL cluster order (GAL1 → GAL10 → GAL7) is preserved on scaffold_1, mirroring the arrangement in C. albicans.

The dot-plot revealed a striking contrast: while gene order within the GAL cluster is conserved, the cluster sits within a region showing a large chromosomal inversion relative to C. albicans. At the whole-genome level, scaffolds 2 and 4 show complete inversions relative to C. albicans chromosomes 2 and 8, and multiple scaffolds align to chromosome 1: evidence of translocations or chromosomal fragmentation.

GAL7’s protein sequence is highly conserved (79–100% identity across Candida), while the surrounding genomic architecture has diverged substantially. This contrast highlights how functional constraint at the sequence level can persist even as chromosomal structure undergoes significant reorganization.

Skills & Tools:

• NCBI Genome Annotation, UniProt BLAST & Align, EMBOSS Needle
• Minimap2 (splice-aware alignment), SAMtools, MUMmer4 / nucmer, mummerplot
• Circa 2.0 (Circos synteny plots)
• Python scripting (coordinate parsing, CSV formatting)
• HPC / SLURM job scheduling (Centaurus cluster)
• BAM/SAM processing, dot-plot interpretation

References:

1. Altschul et al. (1990). Basic local alignment search tool. J. Mol. Biol., 215(3), 403–410.
2. Krzywinski et al. (2009). Circos: An information aesthetic for comparative genomics. Genome Research, 19(9), 1639–1645.
3. Kurtz et al. (2004). Versatile and open software for comparing large genomes. Genome Biology, 5(2), R12.
4. Li, H. (2018). Minimap2: Pairwise alignment for nucleotide sequences. Bioinformatics, 34(18), 3094–3100.
5. Li et al. (2009). The Sequence Alignment/Map format and SAMtools. Bioinformatics, 25(16), 2078–2079.