Cutoffs and k-mers: implications from a transcriptome study in allopolyploid plants.

Area of Science:

• Plant genomics
• Bioinformatics
• Evolutionary biology

Background:

• Transcriptome analysis is vital for studying non-model plants.
• Whole genome duplication (WGD) and allopolyploidy in plants create challenges for gene assembly.
• Divergent homeologous genes complicate de novo and reference-based transcriptome assembly.

Purpose of the Study:

• To develop a successful strategy for assembling transcriptomes from allopolyploid plant species.
• To address challenges posed by homeologous gene divergence in transcriptome assembly.
• To improve gene assembly for evolutionary and ecological studies in non-model plants.

Main Methods:

• Systematic analysis of transcriptome assembly parameters.
• Testing 19 different coverage cutoffs and 20 different k-mer sizes.
• Utilizing 75 bp Illumina sequencing reads for assembly.

Main Results:

• No single parameter set worked for all genes; optimal values varied.
• Assembly success depended on specific k-mer size and coverage combinations.
• Gene expression levels and similarity influenced optimal assembly parameters.

Conclusions:

• Simultaneous consideration of k-mer size and coverage is crucial for allopolyploid transcriptomes.
• This approach maximizes the assembly of full-length expressed sequence tags (ESTs).
• It effectively avoids chimeric assemblies of homeologous and paralogous gene copies.