Zubair Khalid

Virologist/Molecular Biologist | Veterinarian | Bioinformatician

Conventional & Molecular Virology • Vaccine Development • Computational Biology

Dr. Zubair Khalid is a veterinarian and virologist specializing in conventional and molecular virology, vaccine development, and computational biology. Dedicated to advancing animal health through innovative research and multi-omics approaches.

Dr. Zubair Khalid - Veterinarian, Virologist, and Vaccine Development Researcher specializing in Computational Biology, Multi-omics, Animal Health, and Infectious Disease Research

Category: Guides

Geneious Genome Assembly

Geneious genome assembly refers to the process of reconstructing whole or partial genomes from sequencing reads using the Geneious Prime software platform. This guide is for biologists, bioinformaticians, and students who need a practical, source bounded framework to design, execute, and evaluate genome assemblies in Geneious. You will learn core concepts, decision points, a step by step workflow, quality checks, common mistakes, and the limits of interpretation. Throughout, we draw on authoritative resources including the NCBI Bookshelf 1 and the EMBL EBI Training materials 2.

At a Glance

Aspect Key Points
Input data Paired end or long reads in FASTQ, FASTA, or SRA format
Assembly strategies De novo (no reference) or reference guided (mapping then consensus)
Geneious tools Geneious Assembler, De Novo Assembler, Map to Reference
Quality metrics N50, total length, number of contigs, read mapping rate
Common pitfalls Using raw reads without trimming, ignoring GC bias, over trusting short contigs
Interpretation limits Assembly gaps, repeat resolution, and chimeric sequences require careful validation

Core Concepts and Decision Points

Genome assembly in Geneious relies on three fundamental concepts: read pre processing, assembly algorithm selection, and post assembly validation. Before starting, you must decide between de novo and reference guided assembly. De novo assemblies are used when a close reference genome is unavailable or when you aim to discover novel structural variants. Reference guided assemblies are faster and more accurate for organisms with high quality reference genomes 1.

Another critical decision is input read type. Short reads (Illumina) produce high base accuracy but struggle with repeats. Long reads (PacBio, Oxford Nanopore) span repeats but have higher error rates. Geneious supports both, and hybrid approaches often yield the best results. For example, in the complete mitogenome study of Haemaphysalis parva, researchers used a mix of short and long reads to resolve repetitive regions 7.

Read trimming is non negotiable. It removes adapter sequences and low quality bases. Geneious has built in tools, but you can also import pre trimmed reads from external pipelines. The NCBI Sequence Read Archive (SRA) 5 provides raw data for testing and validation.

Practical Workflow for Genome Assembly in Geneious

The following workflow assumes you have sequencing reads and a reference genome (for guided assembly) or no reference (for de novo). We use the Galaxy Training Network 3 as a general framework for reproducible bioinformatics workflows.

Step 1: Data Import and Quality Control

Import your reads into Geneious via File > Import. For FASTQ files, Geneious automatically generates quality statistics. Examine the quality scores, GC content, and read length distribution. Trim low quality regions using Geneious’s “Trim Ends” tool or the BBDuk plugin. Save the trimmed reads as a new list.

Step 2: Choose Assembly Mode

In Geneious Prime, navigate to Tools > Align/Assemble. Select one of:

  • Map to Reference , for reference guided assembly. Choose your reference and select “Generate Consensus Sequence”. The Galaxy Training Network recommends using a high quality reference from NCBI RefSeq 3.
  • De Novo Assemble , for de novo assembly. Choose the correct data type (short reads, long reads, or mixed). Adjust sensitivity settings: higher sensitivity increases memory use but may improve assembly.

Step 3: Run the Assembly

For de novo assembly, set parameters:

  • Minimum overlap identity (default 80%): increase for lower coverage regions.
  • Maximum gap size: adjust based on expected insert sizes.
  • Perform trimming option: enable to remove low quality ends during assembly.

Click OK. The assembly produces contigs and a consensus sequence. Save the result as an assembled genome object.

Step 4: Post Assembly Cleaning

Remove contigs shorter than 500 bp unless you have a specific reason to keep them. Use the “Filter Contigs” option. You may also extract circular sequences (e.g., plasmids, mitogenomes) by using the “Find Repeats” tool. The Galaxy Training Network includes a module on handling bacterial genomes that is directly applicable 3.

Step 5: Annotation (Optional)

Geneious includes plugins for Prokka and other annotation tools. For automated annotation, you can export your assembly to FASTA and use external services like RAST. If annotating within Geneious, verify each feature against databases using BLAST.

Quality Checks and Validation

After assembly, you must evaluate its quality. The Bioconductor project 4 offers R packages (e.g., QuasR, rpt) that can compute metrics from Geneious output BAM files.

Key Metrics

  • N50: the length of the shortest contig at which 50% of the total assembly length is covered. Higher N50 indicates a more contiguous assembly.
  • Total assembly length: should approximate the expected genome size.
  • Number of contigs: fewer is better, but many short contigs may indicate repeats or chimeras.
  • Read mapping rate: map your trimmed reads back to the assembly using Geneious. A rate above 90% suggests high completeness. Lower rates may indicate missing regions or contamination.

Validation Against Known Sequences

If a close reference exists, align your assembly to it using Mauve or the Geneious aligner. Check for large scale rearrangements or missing genes. You can also verify conserved single copy genes. The EMBL EBI Training materials recommend using BUSCO outside of Geneious for universal ortholog completeness 2.

Experimental Confirmation

For critical regions (e.g., antibiotic resistance genes, novel genes), consider PCR or Sanger sequencing. In the study of Pseudomonas aeruginosa ST2407, genomic analysis confirmed the chromosomal location of blaNDM 1 through additional long read sequencing 11.

Common Mistakes and How to Avoid Them

  1. Skipping read trimming: Low quality bases cause mismatches and gaps. Always trim and review QC reports.
  2. Using default parameters without justification: Adjust settings based on your data type. Long reads require different overlap settings than short reads. For example, the duck reovirus study in Cairina moschata used high sensitivity for long read assembly 8.
  3. Ignoring GC bias: Extreme GC content can cause under representation. Check the GC plot. If skewed, consider a correction method or use an assembler that handles GC bias natively.
  4. Over merger of contigs: Geneious can join contigs based on mate pair information but may create chimeras. Inspect joins manually, especially in repetitive regions.
  5. Not documenting assembly parameters: Reproducibility requires you to record the exact settings. Use the Geneious history or save a report.

Limits of Interpretation

A genome assembly is a computational reconstruction, not a perfect representation. The following limits are inherent:

  • Repeat resolution: Tandem and interspersed repeats collapse into single consensus copies, losing true copy number. Long reads help but do not eliminate the problem. The telomeric motif study in Nineta flava highlighted how non canonical repeats can be missed by standard assemblers 9.
  • Strain heterogeneity: In mixed samples or polyploid genomes, the assembly may produce a consensus that does not represent any single molecule. The Staphylococcus pseudintermedius study showed that co infection with multiple strains could produce ambiguous contigs 6.
  • Misassemblies: Chimeric contigs formed by joining unrelated regions can lead to false gene fusions or inversions. Validation with long reads or optical maps is essential.
  • Missing sequences: Small plasmids, very high GC regions, or extremely low coverage areas may be absent. The influenza study in Sri Lanka reported that intra host diversity was partially lost due to low coverage 10.
  • Annotation dependence: Downstream interpretation relies on gene prediction tools, which can be inaccurate. Always manually inspect genes of interest.

Frequently Asked Questions

Q1: Can I use Geneious for assembly of viral genomes?
Yes, Geneious works well for small genomes like viruses. For RNA viruses, use the “Map to Reference” mode with a primer scheme. Always include replicate coverage checks.

Q2: What is the fastest assembly method in Geneious for bacteria?
The Geneious De Novo Assembler with “High Sensitivity/Medium Speed” preset is a good starting point. For large bacterial genomes, you may want to use the SPAdes plugin (third party) and import the result.

Q3: How do I handle mate pair or long jump libraries?
Import them as paired reads and change the insert size in the “Set Paired Read Distances” tool. This information helps resolve repeats. The Galaxy Training Network has a tutorial on mate pair assembly 3.

Q4: My assembly has thousands of small contigs. What should I do?
First, filter contigs shorter than 200 bp. Then check for contamination using tools like Kraken or BLAST against NCBI nt. If contamination is ruled out, try a different assembler or increase sequencing depth.

References and Further Reading

Related Articles