Bdelloid rotifers are small animals well known for their ability to survive in harsh conditions. They can endure extreme dehydration by entering a state of dormancy and then revive when water becomes available again. Scientists have always been interested in the degree to which these animals resist radiation, since radiation can severely damage DNA by causing both of its strands to break. Interestingly, Adineta vaga's ability to survive radiation may be connected to the way in which it tolerates drying out, as dehydration can also damage DNA. The rotifer probably evolved strong DNA protection and repair mechanisms primarily as a response to its severe environment. Its resistance to radiation is merely a byproduct of adapting to the cycles of drying and rehydration.

In order to find out more about this ability, the researchers subjected A. vaga to proton radiation and then examined what effect it had on its chromosomes. The radiation caused serious damage, resulting in hundreds of breaks in the DNA. For the vast majority of organisms, such extensive chromosome damage would be disastrous; it's similar to going through an instruction manual in a shredder and ending up with a large number of separate fragments. The scientists believed that the rotifers would repair their chromosomes before producing offspring. However, what they discovered was unexpected: some of the broken pieces of chromosome did survive and were passed on to the next generations, where the repair process continued.

The discovery took on even greater interest when the scientists examined the offspring of the irradiated rotifers. Rather than observing the same DNA damage in all the generations, they found different sized sections missing and evidence that some of the genetic regions which had been lost were being slowly restored. The researchers proposed that the rotifers use a repair mechanism known as Break-Induced Homologous Extension Repair, or BIHER. In simple terms, a broken piece of DNA can locate a matching sequence on the other chromosome and use it as a guide to rebuild itself. If the damage is beyond what can be repaired in one generation, the repair process can continue in the next.

What makes the rotifer's method of repairing DNA so unusual is that, normally, we consider DNA repair to take place immediately after damage has occurred, with the cell rapidly repairing the damaged DNA. In A. vaga, however, the repair process can be a lot longer. Certain pieces of chromosome survive the initial damage, are inherited, and then serve as the basis for further repairs in later generations. It's as if one generation begins mending a broken book, the next one continues the work, and another eventually completes the job until all the missing parts have been restored.

One reason why this might work is that bdelloid rotifers have holocentric chromosomes. Rather than having a single point, as most chromosomes do, which controls movement during cell division, the ability to move is distributed along the entire length of their chromosomes. As a result, fragments of chromosomes are sometimes able to retain the parts that are necessary for them to move and divide, and this may account for how broken pieces of chromosome can last long enough to be passed on and then repaired.

This finding could have implications that go beyond those concerning rotifers. Scientists are interested in understanding how genomes cope with severe chromosome damage, for example in the case of chromothripsis, where chromosomes break and are then reassembled in new configurations. Investigating A. vaga could help us to learn how cells manage large genetic issues. Moreover, its high resistance to radiation might also prove useful in research into extreme environments, such as the radiation hazards encountered on long space missions. However, these are ideas for future investigation; at present the rotifer's DNA repair mechanism cannot be used to protect people from radiation.

What is perhaps the most interesting aspect of this discovery is that the rotifer's remarkable ability could have developed for a very simple reason. Since an animal living in moss constantly undergoes periods of drying out and then getting wet again, the same mechanisms that protect its DNA when it is dehydrated might also enable it to survive radiation which would damage most other animals. Therefore, the rotifer probably did not evolve solely in order to resist radiation; instead, it became extremely good at coping with the difficult conditions it encounters, and radiation resistance was just a byproduct.

The case of Adineta vaga demonstrates that survival does not always mean avoiding damage; on some occasions it involves being able to recover from damage which appears to be impossible to repair. This small animal can survive having its chromosomes broken into many pieces, retain some of those pieces, pass them on to its offspring, and continue to repair them over successive generations. Although we generally regard DNA as something fragile that must always remain intact, the rotifer indicates an alternative possibility: perhaps a genome does not have to be perfect in order to survive, since sometimes the pieces last long enough for life to reassemble them, bit by bit, generation after generation.

Ream more at https://www.science.org/doi/10.1126/sciadv.aee0891

Followings are the quick way to check if a value exist in an array.

do_something if 'flour' ~~ @ingredients   # ~~ operand.   BEWARE: it is broken.

do_something if grep {$_ eq 'flour'} @ingredients # grep (slower than 'any')

do_something if any {$_ eq 'flour'} @ingredients # List::MoreUtils / Util::Any

do_something if any(@ingredients) eq 'flour'   # use syntax 'junction';

do_something if @ingredients->contains('flour')   # added with autobox

What is Genome Assembly?

Genome assembly involves piecing together short DNA sequence reads generated by sequencing platforms (e.g., Illumina, PacBio, Oxford Nanopore) into longer, contiguous sequences called contigs. This can be performed as:

• De Novo Assembly: Without a reference genome.
• Reference-Guided Assembly: Using a reference genome to guide the assembly process.

Step 1: Preparing Your Data

Before starting the assembly, ensure that your raw sequencing data is high quality.

1. Input Data

   • Short Reads: Illumina sequencing generates short, accurate reads ideal for scaffolding.
   • Long Reads: PacBio and Nanopore sequencing provide long reads for resolving repetitive regions.

2. Quality Control (QC)
   Use tools like FastQC or MultiQC to assess the quality of your reads:

   fastqc reads.fastq multiqc .

   Look for issues like low-quality bases, adapter contamination, or overrepresented sequences.

3. Read Trimming and Filtering
   Trim low-quality bases and adapters using Trimmomatic or Cutadapt:

   trimmomatic PE reads_R1.fastq reads_R2.fastq trimmed_R1.fastq trimmed_R2.fastq \ ILLUMINACLIP:adapters.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:20 MINLEN:36

Step 2: Choosing an Assembly Strategy

Select an assembly strategy based on your data type:

• Short-Read Assemblers:

  • SPAdes: Popular for microbial genomes.
  • Velvet: Fast for smaller genomes.

• Long-Read Assemblers:

  • Canu: Ideal for long-read datasets.
  • Flye: Versatile for small and large genomes.

• Hybrid Assemblers:

  • MaSuRCA: Combines short and long reads.
  • Unicycler: Optimized for bacterial genomes.

Step 3: Running the Assembly

3.1. SPAdes (Short-Read Assembly)

SPAdes is an excellent choice for small genomes, such as bacteria.

The output includes assembled contigs (contigs.fasta) and scaffolds (scaffolds.fasta).

3.2. Canu (Long-Read Assembly)

Canu is designed for high-error long reads from PacBio or Nanopore.

The output will be in canu_output/genome.contigs.fasta.

3.3. Hybrid Assembly with Unicycler

Unicycler combines short and long reads for improved assemblies.

Step 4: Assessing Assembly Quality

After assembly, evaluate its quality using the following tools:

1. QUAST
   QUAST generates assembly statistics, such as N50, genome size, and GC content:

   quast contigs.fasta -o quast_output

2. BUSCO
   BUSCO checks genome completeness by identifying conserved genes:

   busco -i contigs.fasta -o busco_output -l fungi_odb10 -m genome

3. Assembly Graph Visualization
   Visualize assembly graphs with Bandage:

   Bandage load assembly_graph.gfa

Step 5: Post-Assembly Steps

1. Polishing
   Improve assembly accuracy using tools like Pilon (for short reads) or Racon (for long reads).

   racon long_reads.fasta mapped_reads.sam contigs.fasta > polished_contigs.fasta

2. Scaffolding
   Link contigs into scaffolds using tools like SSPACE or Opera-LG if required.

3. Annotation
   Annotate the assembled genome using Prokka for prokaryotes or Maker for eukaryotes.

   prokka --outdir annotation_output --prefix genome contigs.fasta

Step 6: Sharing and Archiving

1. Submit to Public Repositories
   Share your assembly in databases like NCBI GenBank, ENA, or DDBJ.

2. Metadata Preparation
   Include detailed metadata for your submission, such as organism name, sequencing platform, and coverage.

Best Practices

• Always perform quality checks at each stage to ensure data integrity.
• Use multiple tools to cross-validate results when working with complex genomes.
• Document parameters and software versions for reproducibility.

Conclusion

Genome assembly is a powerful process that transforms raw sequencing data into a coherent representation of an organism’s genome. By following this step-by-step guide, you can successfully assemble genomes and uncover valuable biological insights. Whether you’re assembling a microbial genome or tackling the complexities of a eukaryotic genome, these tools and strategies will set you on the path to success.

• nodes, which are labeled by sequences and ids
• edges, which connect two nodes via either of their respective ends
• paths, describe genomes, sequence alignments, and annotations (such as gene models and transcripts) as walks through nodes connected by edges

Address of the bookmark: https://github.com/vgteam/vg

http://www.i-programmer.info/news/181-algorithms/4537-a-new-dna-sequence-search-compressive-genomics.html

http://en.wikipedia.org/wiki/Compression_of_Genomic_Re-Sequencing_Data

http://www.nature.com/nbt/journal/v30/n7/full/nbt.2241.html

http://bioinformatics.oxfordjournals.org/content/29/13/i283.full

http://groups.csail.mit.edu/cb/cast/

Address of the bookmark: https://urgi.versailles.inra.fr/Tools/REPET/TEannot-tuto

http://katjulab.com/research.html

Address of the bookmark: https://gvolante.riken.jp/analysis.html

GenBank release 257.0 (8/15/2023) is now available on the NCBI FTP site. This release has 25.10 trillion bases and 3.69 billion records.

The current release has:

• 246,119,175 traditional records containing 2,112,058,517,945 base pairs of sequence data
• 2,631,493,489 WGS records containing 22,294,446,104,543 base pairs of sequence data
• 686,271,945 bulk-oriented TSA records containing 646,176,166,908 base pairs of sequence data
• 124,421,006 bulk-oriented TLS records containing 48,289,699,026 base pairs of sequence data