Ultra-fast genome comparison for large-scale genomic experiments

Chromosome-scale pearl millet genomes reveal CLAMT1b as key determinant of strigolactone pattern and Striga susceptibility

The yield of pearl millet, a resilient cereal crop crucial for African food security, is severely impacted by the root parasitic weed Striga hermonthica, which requires host-released hormones, called strigolactones (SLs), for seed germination. Herein, we identify four SLs present in the Striga-susceptible line SOSAT-C88-P10 (P10) but absent in the resistant 29Aw (Aw). We generate chromosome-scale genome assemblies, including four gapless chromosomes for each line. The Striga-resistant Aw lacks a 0.7 Mb genome segment containing two putative CARLACTONOIC ACID METHYLTRANSFERASE1 (CLAMT1) genes, which may contribute to SL biosynthesis. Functional assays show that P10CLAMT1b produces the SL-biosynthesis intermediate methyl carlactonoate (MeCLA) and that MeCLA is the precursor of P10-specific SLs. Screening a diverse pearl millet panel confirms the pivotal role of the CLAMT1 section for SL diversity and Striga susceptibility. Our results reveal a reason for Striga susceptibility in pearl millet and pave the way for generating resistant lines through marker-assisted breeding or direct genetic modification.

Avirulent UG10 Entamoeba histolytica mutant derived from HM-1:IMSS strain shows limited genome variability and aberrant 5-methyl cytosine genomic distribution

Entamoeba histolytica, an intestinal parasite of global significance, poses substantial health risks with its associated high morbidity and mortality rates. Despite the current repertoire of molecular tools for the study of gene function in, the regulatory mechanisms governing its pathogenicity remain largely unexplored. This knowledge gap underscores the need to elucidate key genetic determinants orchestrating cellular functions critical to its virulence. Previously, our group generated an avirulent strain, termed UG10, with the same genetic background as the HM1:IMSS strain. UG10 strain, despite showing normal expression levels of well-known virulence factors, was unable to perform in-vitro and in-vivo activities related to amoebic virulence. In this study, we aimed to uncover the genome-wide modifications that rendered the avirulent phenotype of the UG10 strain through whole-genome sequencing. As a complementary approach, we conducted Methylated DNA Immunoprecipitation coupled with sequencing (MeDIP-seq) analysis on both the highly virulent HM1:IMSS strain and the low-virulence UG10 strain to uncover the genome-wide methylation profile. These dual methodologies revealed two aspects of the UG10 avirulent strain. One is the random integration of fragments from the ribosomal gene cluster and tRNA genes, ranging from 120 to 400 bp; and secondly, a clear, enriched methylation profile in the coding and non-coding strand relative to the start codon sequence in genes encoding small GTPases, which is associated with the previously described avirulent phenotype. This study provides the foundation to explore other genetic and epigenetic regulatory circuitries in E. histolytica and novel targets to understand the pathogenic mechanism of this parasite.

A comprehensive analysis integrating phenotypic assessment uncovering thornless cultivar lineages in Aralia elata

Aralia elata is an Araliaceae woody plant species found in Northeastern Asia. To understand how genetic pools are distributed for A.elata clones, we were to analyze the population structure of A.elata cultivars and identify how these are correlated with thorn-related phenotype which determines the utility of A.elata. We found that the de novo assembled genome of 'Yeongchun' shared major genomic compartments with the public A.elata genome assembled from the wild-type from China. To identify the population structure of the 32 Korean and Japanese cultivars, we identified 44 SSR markers and revealed three main sub-clusters using ΔK analysis with one isolated cultivar. Machine-learning based clustering with thorn-related phenotype correlated moderately with population structure based on SSR analysis suggested multi-layered genetic regulation of thorn-related phenotypes. Thus, we revealed genetic lineage of A.elata and uncovered isolated cultivar which can provide new genetic material for further breeding.

Chromosome-scale assembly of the wild wheat relative Aegilops umbellulata

Wild wheat relatives have been explored in plant breeding to increase the genetic diversity of bread wheat, one of the most important food crops. Aegilops umbellulata is a diploid U genome-containing grass species that serves as a genetic reservoir for wheat improvement. In this study, we report the construction of a chromosome-scale reference assembly of Ae. umbellulata accession TA1851 based on corrected PacBio HiFi reads and chromosome conformation capture. The total assembly size was 4.25 Gb with a contig N50 of 17.7 Mb. In total, 36,268 gene models were predicted. We benchmarked the performance of hifiasm and LJA, two of the most widely used assemblers using standard and corrected HiFi reads, revealing a positive effect of corrected input reads. Comparative genome analysis confirmed substantial chromosome rearrangements in Ae. umbellulata compared to bread wheat. In summary, the Ae. umbellulata assembly provides a resource for comparative genomics in Triticeae and for the discovery of agriculturally important genes.

Insights from a chum salmon (Oncorhynchus keta) genome assembly regarding whole-genome duplication and nucleotide variation influencing gene function

Chum salmon are ecologically important to Pacific Ocean ecosystems and commercially important to fisheries. To improve the genetic resources available for this species, we sequenced and assembled the genome of a male chum salmon using Oxford Nanopore read technology and the Flye genome assembly software (contig N50: ∼2 Mbp, complete BUSCOs: ∼98.1%). We also resequenced the genomes of 59 chum salmon from hatchery sources to better characterize the genome assembly and the diversity of nucleotide variants impacting phenotype variation. With genomic sequences from a doubled haploid individual, we were able to identify regions of the genome assembly that have been collapsed due to high sequence similarity between homeologous (duplicated) chromosomes. The homeologous chromosomes are relics of an ancient salmonid-specific genome duplication. These regions were enriched with genes whose functions are related to the immune system and responses to toxins. From analyzing nucleotide variant annotations of the resequenced genomes, we were also able to identify genes that have increased levels of variants thought to moderately impact gene function. Genes related to the immune system and the detection of chemical stimuli (olfaction) had increased levels of these variants based on a gene ontology enrichment analysis. The tandem organization of many of the enriched genes raises the question of why they have this organization.

Einkorn genomics sheds light on history of the oldest domesticated wheat

Einkorn (Triticum monococcum) was the first domesticated wheat species, and was central to the birth of agriculture and the Neolithic Revolution in the Fertile Crescent around 10,000 years ago1,2. Here we generate and analyse 5.2-Gb genome assemblies for wild and domesticated einkorn, including completely assembled centromeres. Einkorn centromeres are highly dynamic, showing evidence of ancient and recent centromere shifts caused by structural rearrangements. Whole-genome sequencing analysis of a diversity panel uncovered the population structure and evolutionary history of einkorn, revealing complex patterns of hybridizations and introgressions after the dispersal of domesticated einkorn from the Fertile Crescent. We also show that around 1% of the modern bread wheat (Triticum aestivum) A subgenome originates from einkorn. These resources and findings highlight the history of einkorn evolution and provide a basis to accelerate the genomics-assisted improvement of einkorn and bread wheat.

Multi-Omics Pipeline and Omics-Integration Approach to Decipher Plant's Abiotic Stress Tolerance Responses

The present day's ongoing global warming and climate change adversely affect plants through imposing environmental (abiotic) stresses and disease pressure. The major abiotic factors such as drought, heat, cold, salinity, etc., hamper a plant's innate growth and development, resulting in reduced yield and quality, with the possibility of undesired traits. In the 21st century, the advent of high-throughput sequencing tools, state-of-the-art biotechnological techniques and bioinformatic analyzing pipelines led to the easy characterization of plant traits for abiotic stress response and tolerance mechanisms by applying the 'omics' toolbox. Panomics pipeline including genomics, transcriptomics, proteomics, metabolomics, epigenomics, proteogenomics, interactomics, ionomics, phenomics, etc., have become very handy nowadays. This is important to produce climate-smart future crops with a proper understanding of the molecular mechanisms of abiotic stress responses by the plant's genes, transcripts, proteins, epigenome, cellular metabolic circuits and resultant phenotype. Instead of mono-omics, two or more (hence 'multi-omics') integrated-omics approaches can decipher the plant's abiotic stress tolerance response very well. Multi-omics-characterized plants can be used as potent genetic resources to incorporate into the future breeding program. For the practical utility of crop improvement, multi-omics approaches for particular abiotic stress tolerance can be combined with genome-assisted breeding (GAB) by being pyramided with improved crop yield, food quality and associated agronomic traits and can open a new era of omics-assisted breeding. Thus, multi-omics pipelines together are able to decipher molecular processes, biomarkers, targets for genetic engineering, regulatory networks and precision agriculture solutions for a crop's variable abiotic stress tolerance to ensure food security under changing environmental circumstances.

Population genomics reveals mechanisms and dynamics of de novo expressed open reading frame emergence in Drosophila melanogaster

Novel genes are essential for evolutionary innovations and differ substantially even between closely related species. Recently, multiple studies across many taxa showed that some novel genes arise de novo, that is, from previously noncoding DNA. To characterize the underlying mutations that allowed de novo gene emergence and their order of occurrence, homologous regions must be detected within noncoding sequences in closely related sister genomes. So far, most studies do not detect noncoding homologs of de novo genes because of incomplete assemblies and annotations, and long evolutionary distances separating genomes. Here, we overcome these issues by searching for de novo expressed open reading frames (neORFs), the not-yet fixed precursors of de novo genes that emerged within a single species. We sequenced and assembled genomes with long-read technology and the corresponding transcriptomes from inbred lines of Drosophila melanogaster, derived from seven geographically diverse populations. We found line-specific neORFs in abundance but few neORFs shared by lines, suggesting a rapid turnover. Gain and loss of transcription is more frequent than the creation of ORFs, for example, by forming new start and stop codons. Consequently, the gain of ORFs becomes rate limiting and is frequently the initial step in neORFs emergence. Furthermore, transposable elements (TEs) are major drivers for intragenomic duplications of neORFs, yet TE insertions are less important for the emergence of neORFs. However, highly mutable genomic regions around TEs provide new features that enable gene birth. In conclusion, neORFs have a high birth-death rate, are rapidly purged, but surviving neORFs spread neutrally through populations and within genomes.

Genome-Wide Pathway Exploration of the Epidermidibacterium keratini EPI-7T

Functional cosmetics industries using skin microbiome screening and beneficial materials isolated from key microorganisms are receiving increasing attention. Since Epidermidibacterium keratini EPI-7T was first discovered in human skin, previous studies have confirmed that it can produce a new pyrimidine compound, 1,1'-biuracil, having anti-aging effects on human skin. Therefore, we conducted genomic analyses to judge the use value of E. keratini EPI-7T and provide up-to-date information. Whole-genome sequencing analysis of E. keratini EPI-7T was performed to generate new complete genome and annotation information. E. keratini EPI-7T genome was subjected to comparative genomic analysis with a group of closely-related strains and skin flora strains through bioinformatic analysis. Furthermore, based on annotation information, we explored metabolic pathways for valuable substances that can be used in functional cosmetics. In this study, the whole-genome sequencing (WGS) and annotation results of E. keratini EPI-7T were improved, and through comparative analysis, it was confirmed that the E. keratini EPI-7T has more metabolite-related genes than comparison strains. In addition, we annotated the vital genes for biosynthesis of 20 amino acids, orotic acid, riboflavin (B2) and chorismate. In particular, we were able to prospect that orotic acid could accumulate inside E. keratini EPI-7T under uracil-enriched conditions. Therefore, through a genomics approach, this study aims to provide genetic information for the hidden potential of E. keratini EPI-7T and the strain development and biotechnology utilization to be conducted in further studies.

A chromosome-scale genome sequence of sudangrass (Sorghum sudanense) highlights the genome evolution and regulation of dhurrin biosynthesis

Sudangrass is more similar to US commercial sorghums than to cultivated sorghums from Africa sequence-wise and contain significantly lower dhurrin than sorghums. CYP79A1 is linked to dhurrin content in sorghum. Sudangrass [Sorghum sudanense (Piper) Stapf] is a hybrid between grain sorghum and its wild relative S. bicolor ssp. verticilliflorum and is grown as a forage crop due to its high biomass production and low dhurrin content compared to sorghum. In this study, we sequenced the sudangrass genome and showed that the assembled genome was 715.95 Mb with 35,243 protein-coding genes. Phylogenetic analysis with whole genome proteomes demonstrated that the sudangrass genome was more similar to US commercial sorghums than to its wild relatives and cultivated sorghums from Africa. We confirmed that at seedling stage, sudangrass accessions contained significantly lower dhurrin as measured by hydrocyanic acid potential (HCN-p) than cultivated sorghum accessions. Genome-wide association study identified a QTL most tightly associated with HCN-p and the linked SNPs were located in the 3' UTR of Sobic.001G012300 which encodes CYP79A1, the enzyme that catalyzes the first step of dhurrin biosynthesis. As in other grasses such as maize and rice, we also found that copia/gypsy long terminal repeat (LTR) retrotransposons were more abundant in cultivated than in wild sorghums, implying that crop domestication in the grasses was accompanied by increased copia/gypsy LTR retrotransposon insertions in the genomes.

Robertsonian fusion and centromere repositioning contributed to the formation of satellite-free centromeres during the evolution of zebras

Centromeres are epigenetically specified by the histone H3 variant CENP-A and typically associated to highly repetitive satellite DNA. We previously discovered natural satellite-free neocentromeres in Equus caballus and E. asinus. Here, through ChIP-seq with an anti-CENP-A antibody, we found an extraordinarily high number of centromeres lacking satellite DNA in the zebras E. burchelli (15 of 22) and E. grevyi (13 of 23), demonstrating that the absence of satellite DNA at the majority of centromeres is compatible with genome stability and species survival and challenging the role of satellite DNA in centromere function. Nine satellite-free centromeres are shared between the two species in agreement with their recent separation. We assembled all centromeric regions and improved the reference genome of E. burchelli. Sequence analysis of the CENP-A binding domains revealed that they are LINE-1 and AT-rich with four of them showing DNA amplification. In the two zebras, satellite-free centromeres emerged from centromere repositioning or following Robertsonian fusion. In five chromosomes, the centromeric function arose near the fusion points, which are located within regions marked by traces of ancestral pericentromeric sequences. Therefore, besides centromere repositioning, Robertsonian fusions are an important source of satellite-free centromeres during evolution. Finally, in one case, a satellite-free centromere was seeded on an inversion breakpoint. At eleven chromosomes, whose primary constrictions seemed to be associated to satellite repeats by cytogenetic analysis, satellite-free neocentromeres were instead located near the ancestral inactivated satellite-based centromeres, therefore, the centromeric function has shifted away from a satellite repeat containing locus to a satellite-free new position.

Chromosome-level genome assembly of Plazaster borealis sheds light on the morphogenesis of multiarmed starfish and its regenerative capacity

Background

Plazaster borealis has a unique morphology, displaying multiple arms with a clear distinction between disk and arms, rather than displaying pentaradial symmetry, a remarkable characteristic of echinoderms. Herein we report the first chromosome-level reference genome of P. borealis and an essential tool to further investigate the basis of the divergent morphology.

Findings

In total, 57.76 Gb of a long read and 70.83 Gb of short-read data were generated to assemble a de novo 561-Mb reference genome of P. borealis, and Hi-C sequencing data (57.47 Gb) were used for scaffolding into 22 chromosomal scaffolds comprising 92.38% of the genome. The genome completeness estimated by BUSCO was 98.0% using the metazoan set, indicating a high-quality assembly. Through the comparative genome analysis, we identified evolutionary accelerated genes known to be involved in morphogenesis and regeneration, suggesting their potential role in shaping body pattern and capacity of regeneration.

Conclusion

This first chromosome-level genome assembly of P. borealis provides fundamental insights into echinoderm biology, as well as the genomic mechanism underlying its unique morphology and regeneration.

Chromosome-Scale, Haplotype-Resolved Genome Assembly of Non-Sex-Reversal Females of Swamp Eel Using High-Fidelity Long Reads and Hi-C Data

The Asian swamp eel (Monopterus albus) is an excellent model species for studying sex change and chromosome evolution. M. albus is also widely reared in East Asia and South-East Asia because of its great nutritional value. The low fecundity of this species (about 300 eggs per fish) greatly hinders fries production and breeding programs. Interestingly, about 3–5% of the eels could remain as females for 3 years and lay more than 3,000 eggs per fish, which are referred to as non-sex-reversal (NSR) females. Here, we presented a new chromosome-level genome assembly of such NSR females using Illumina, HiFi, and Hi-C sequencing technologies. The new assembly (Mal.V2_NSR) is 838.39 Mb in length, and the N50 of the contigs is 49.8 Mb. Compared with the previous assembly obtained using the continuous long-read sequencing technology (Mal.V1_CLR), we found a remarkable increase of continuity in the new assembly Mal.V2_NSR with a 20-times longer contig N50. Chromosomes 2 and 12 were assembled into a single contig, respectively. Meanwhile, two highly contiguous haplotype assemblies were also obtained, with contig N50 being 14.54 and 12.13 Mb, respectively. BUSCO and Merqury analyses indicate completeness and high accuracy of these three assemblies. A comparative genomic analysis revealed substantial structural variations (SVs) between Mal.V2_NSR and Mal.V1_CLR and two phased haplotype assemblies, as well as whole chromosome fusion events when compared with the zig-zag eel. Additionally, our newly obtained assembly provides a genomic view of sex-related genes and a complete landscape of the MHC genes. Therefore, these high-quality genome assemblies would provide great help for future breeding works of the swamp eel, and it is a valuable new reference for genetic and genomic studies of this species.

Evolutionary diversification of the canonical Wnt signaling effector TCF/LEF in chordates

Wnt signaling is essential during animal development and regeneration, but also plays an important role in diseases such as cancer and diabetes. The canonical Wnt signaling pathway is one of the most conserved signaling cascades in the animal kingdom, with the T‐cell factor/lymphoid enhancer factor (TCF/LEF) proteins being the major mediators of Wnt/β‐catenin‐regulated gene expression. In comparison with invertebrates, vertebrates possess a high diversity of TCF/LEF family genes, implicating this as a possible key change to Wnt signaling at the evolutionary origin of vertebrates. However, the precise nature of this diversification is only poorly understood. The aim of this study is to clarify orthology, paralogy, and isoform relationships within the TCF/LEF gene family within chordates via in silico comparative study of TCF/LEF gene structure, molecular phylogeny, and gene synteny. Our results support the notion that the four TCF/LEF paralog subfamilies in jawed vertebrates (gnathostomes) evolved via the two rounds of whole‐genome duplications that occurred during early vertebrate evolution. Importantly, gene structure comparisons and synteny analysis of jawless vertebrate (cyclostome) TCFs suggest that a TCF7L2‐like form of gene structure is a close proxy for the ancestral vertebrate structure. In conclusion, we propose a detailed evolutionary path based on a new pre‐whole‐genome duplication vertebrate TCF gene model. This ancestor gene model highlights the chordate and vertebrate innovations of TCF/LEF gene structure, providing the foundation for understanding the role of Wnt/β‐catenin signaling in vertebrate evolution.

Chromosome-level genome assembly of the shuttles hoppfish, Periophthalmus modestus

Background

The shuttles hoppfish (mudskipper), Periophthalmus modestus, is one of the mudskippers, which are the largest group of amphibious teleost fishes, which are uniquely adapted to live on mudflats. Because mudskippers can survive on land for extended periods by breathing through their skin and through the lining of the mouth and throat, they were evaluated as a model for the evolutionary sea-land transition of Devonian protoamphibians, ancestors of all present tetrapods.

Results

A total of 39.6, 80.2, 52.9, and 33.3 Gb of Illumina, Pacific Biosciences, 10X linked, and Hi-C data, respectively, was assembled into 1,419 scaffolds with an N50 length of 33 Mb and BUSCO score of 96.6%. The assembly covered 117% of the estimated genome size (729 Mb) and included 23 pseudo-chromosomes anchored by a Hi-C contact map, which corresponded to the top 23 longest scaffolds above 20 Mb and close to the estimated one. Of the genome, 43.8% were various repetitive elements such as DNAs, tandem repeats, long interspersed nuclear elements, and simple repeats. Ab initio and homology-based gene prediction identified 30,505 genes, of which 94% had homology to the 14 Actinopterygii transcriptomes and 89% and 85% to Pfam familes and InterPro domains, respectively. Comparative genomics with 15 Actinopterygii species identified 59,448 gene families of which 12% were only in P. modestus.

Conclusions

We present the high quality of the first genome assembly and gene annotation of the shuttles hoppfish. It will provide a valuable resource for further studies on sea-land transition, bimodal respiration, nitrogen excretion, osmoregulation, thermoregulation, vision, and mechanoreception.

Role of TlyA in the Biology of Uncultivable Mycobacteria

TlyA proteins are related to distinct functions in a diverse spectrum of bacterial pathogens including mycobacterial spp. There are several annotated proteins function as hemolysin or pore forming molecules that play an important role in the virulence of pathogenic organisms. Many studies reported the dual activity of mycobacterial TlyA as 'hemolysin' and 'S-adenosylmethionine dependent rRNA methylase'. To act as a hemolysin, a sequence must have a signal sequence and transmembrane segment which helps the protein to enter the extracellular environment. Interestingly, the mycobacterial tlyA has neither a traditional signal sequences of general/sec/tat pathways nor any transmembrane segments are present. Still it can reach the extracellular milieu with the help of non-classical signal mechanisms. Also, retention of tlyA in cultivable mycobacterial pathogens (such as Mycobacterium tuberculosis and M. marinum) as well as uncultivated mycobacterial pathogens despite their extreme reductive evolution (such as M. leprae, M. lepromatosis and M. uberis) suggests its crucial role in evolutionary biology of pathogenic mycobacteria. Numerous virulence factors have been characterised from the uncultivable mycobacteria but the information of TlyA protein is still limited in terms of molecular and structural characterisation. The genomic insights offered by comparative analysis of TlyA sequences and its conserved domains reveal its pore forming activity which further confirms its role as a virulence protein, particularly in uncultivable mycobacteria. Therefore, this review presents a comparative analysis of mycobacterial TlyA family by sequence homology and alignment to improve our understanding of this unconventional hemolysin and RNA methyltransferase TlyA of uncultivable mycobacteria.

The pink salmon genome: Uncovering the genomic consequences of a two-year life cycle

Pink salmon (Oncorhynchus gorbuscha) adults are the smallest of the five Pacific salmon native to the western Pacific Ocean. Pink salmon are also the most abundant of these species and account for a large proportion of the commercial value of the salmon fishery worldwide. A two-year life history of pink salmon generates temporally isolated populations that spawn either in even-years or odd-years. To uncover the influence of this genetic isolation, reference genome assemblies were generated for each year-class and whole genome re-sequencing data was collected from salmon of both year-classes. The salmon were sampled from six Canadian rivers and one Japanese river. At multiple centromeres we identified peaks of Fst between year-classes that were millions of base-pairs long. The largest Fst peak was also associated with a million base-pair chromosomal polymorphism found in the odd-year genome near a centromere. These Fst peaks may be the result of a centromere drive or a combination of reduced recombination and genetic drift, and they could influence speciation. Other regions of the genome influenced by odd-year and even-year temporal isolation and tentatively under selection were mostly associated with genes related to immune function, organ development/maintenance, and behaviour.

Telomere-to-telomere genome assembly of asparaginase-producing Trichoderma simmonsii

BACKGROUND

Trichoderma is a genus of fungi in the family Hypocreaceae and includes species known to produce enzymes with commercial use. They are largely found in soil and terrestrial plants. Recently, Trichoderma simmonsii isolated from decaying bark and decorticated wood was newly identified in the Harzianum clade of Trichoderma. Due to a wide range of applications in agriculture and other industries, genomes of at least 12 Trichoderma spp. have been studied. Moreover, antifungal and enzymatic activities have been extensively characterized in Trichoderma spp. However, the genomic information and bioactivities of T. simmonsii from a particular marine-derived isolate remain largely unknown. While we screened for asparaginase-producing fungi, we observed that T. simmonsii GH-Sj1 strain isolated from edible kelp produced asparaginase. In this study, we report a draft genome of T. simmonsii GH-Sj1 using Illumina and Oxford Nanopore technologies. Furthermore, to facilitate biotechnological applications of this species, RNA-sequencing was performed to elucidate the transcriptional profile of T. simmonsii GH-Sj1 in response to asparaginase-rich conditions.

RESULTS

We generated ~ 14 Gb of sequencing data assembled in a ~ 40 Mb genome. The T. simmonsii GH-Sj1 genome consisted of seven telomere-to-telomere scaffolds with no sequencing gaps, where the N50 length was 6.4 Mb. The total number of protein-coding genes was 13,120, constituting ~ 99% of the genome. The genome harbored 176 tRNAs, which encode a full set of 20 amino acids. In addition, it had an rRNA repeat region consisting of seven repeats of the 18S-ITS1-5.8S-ITS2-26S cluster. The T. simmonsii genome also harbored 7 putative asparaginase-encoding genes with potential medical applications. Using RNA-sequencing analysis, we found that 3 genes among the 7 putative genes were significantly upregulated under asparaginase-rich conditions.

CONCLUSIONS

The genome and transcriptome of T. simmonsii GH-Sj1 established in the current work represent valuable resources for future comparative studies on fungal genomes and asparaginase production.

Pyramiding dsRNAs increases phytonematode tolerance in cotton plants

Host-derived suppression of nematode essential genes decreases reproduction of Meloidogyne incognita in cotton. Root-knot nematodes (RKN) represent one of the most damaging plant-parasitic nematode genera worldwide. RNAi-mediated suppression of essential nematode genes provides a novel biotechnological strategy for the development of sustainable pest-control methods. Here, we used a Host Induced Gene Silencing (HIGS) approach by stacking dsRNA sequences into a T-DNA construct to target three essential RKN genes: cysteine protease (Mi-cpl), isocitrate lyase (Mi-icl), and splicing factor (Mi-sf), called dsMinc1, driven by the pUceS8.3 constitutive soybean promoter. Transgenic dsMinc1-T₄ plants infected with Meloidogyne incognita showed a significant reduction in gall formation (57-64%) and egg masses production (58-67%), as well as in the estimated reproduction factor (60-78%), compared with the susceptible non-transgenic cultivar. Galls of the RNAi lines are smaller than the wild-type (WT) plants, whose root systems exhibited multiple well-developed root swellings. Transcript levels of the three RKN-targeted genes decreased 13- to 40-fold in nematodes from transgenic cotton galls, compared with those from control WT galls. Finally, the development of non-feeding males in transgenic plants was 2-6 times higher than in WT plants, indicating a stressful environment for nematode development after RKN gene silencing. Data strongly support that HIGS of essential RKN genes is an effective strategy to improve cotton plant tolerance. This study presents the first application of dsRNA sequences to target multiple genes to promote M. incognita tolerance in cotton without phenotypic penalty in transgenic plants.

Unraveling Genome Evolution Throughout Visual Analysis: The XCout Portal

Due to major breakthroughs in sequencing technologies throughout the last decades, the time and cost per sequencing experiment have reduced drastically, overcoming the data generation barrier during the early genomic era. Such a shift has encouraged the scientific community to develop new computational methods that are able to compare large genomic sequences, thus enabling large-scale studies of genome evolution. The field of comparative genomics has proven itself invaluable for studying the evolutionary mechanisms and the forces driving genome evolution. In this line, a full genome comparison study between 2 species requires a quadratic number of comparisons in terms of the number of sequences (around 400 chromosome comparisons in the case of mammalian genomes); however, when studying conserved syntenies or evolutionary rearrangements, many sequence comparisons can be skipped for not all will contain significant signals. Subsequently, the scientific community has developed fast heuristics to perform multiple pairwise comparisons between large sequences to determine whether significant sets of conserved similarities exist. The data generation problem is no longer an issue, yet the limitations have shifted toward the analysis of such massive data. Therefore, we present XCout, a Web-based visual analytics application for multiple genome comparisons designed to improve the analysis of large-scale evolutionary studies using novel techniques in Web visualization. XCout enables to work on hundreds of comparisons at once, thus reducing the time of the analysis by identifying significant signals between chromosomes across multiple species. Among others, XCout introduces several techniques to aid in the analysis of large-scale genome rearrangements, particularly (1) an interactive heatmap interface to display comparisons using automatic color scales based on similarity thresholds to ease detection at first sight, (2) an overlay system to detect individual signal contributions between chromosomes, (3) a tracking tool to trace conserved blocks across different species to perform evolutionary studies, and (4) a search engine to search annotations throughout different species.

Tandem repeat interval pattern identifies animal taxa

MOTIVATION

We discover that maximality of information content among intervals of Tandem Repeats (TRs) in animal genome segregates over taxa such that taxa identification becomes swift and accurate. Successive TRs of a motif occur at intervals over the sequence, forming a trail of TRs of the motif across the genome. We present a method, Tandem Repeat Information Mining (TRIM), that mines 4k number of TR trails of all k length motifs from a whole genome sequence and extracts the information content within intervals of the trails. TRIM vector formed from the ordered set of interval entropies becomes instrumental for genome segregation.

RESULTS

Reconstruction of correct phylogeny for animals from whole genome sequences proves precision of TRIM. Identification of animal taxa by TRIM vector upon feature selection is the most significant achievement. These suggest Tandem Repeat Interval Pattern (TRIP) is a taxa-specific constitutional characteristic in animal genome.

AVAILABILITY

Source and executable code of TRIM along with usage manual are made available at https://github.com/BB-BiG/TRIM.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.