A Predictive Approach to Infer the Activity and Natural Variation of Retrotransposon Families in Plants

Shotgun sequencing of airborne eDNA achieves rapid assessment of whole biomes, population genetics and genomic variation

Biodiversity and its associated genetic diversity are being lost at an unprecedented rate. Simultaneously, the distributions of flora, fauna, fungi, microbes and pathogens are rapidly changing. Novel technology can help to capture and record genetic diversity before it is lost and to measure population shifts and pathogen distributions. Here we report the rapid application of shotgun long-read environmental DNA (eDNA) analysis for non-invasive biodiversity, genetic diversity and pathogen assessments from air. We also compared air eDNA with water and soil eDNA. Coupling long-read sequencing with established cloud-based biodiversity pipelines enabled a 2-day turnaround from airborne sample collection to completed analysis by a single investigator. To determine the full utility of airborne eDNA, we also conducted a local bioinformatic analysis and deep short-read shotgun sequencing. From outdoor air eDNA alone, comprehensive genetic analysis was performed, including population genetics (phylogenetic placement) of a charismatic mammal (bobcat, Lynx rufus) and a venomous spider (golden silk orb weaver, Trichonephila clavipes), and haplotyping humans (Homo sapiens) from natural complex community settings, such as subtropical forests and temperate locations. The rich datasets also enabled deeper analysis of specific species and genomic regions of interest, including viral variant calling, human variant analysis and antimicrobial resistance gene surveillance from airborne DNA. Our results highlight the speed, versatility and specificity of pan-biodiversity monitoring via non-invasive eDNA sampling using current benchtop/portable and cloud-based approaches. Furthermore, they reveal the future feasibility of scaling down (equipment and temporally) these approaches for near real-time analysis. Together these approaches can enable rapid simultaneous detection of all life and its genetic diversity from air, water and sediment samples for unbiased non-targeted information-rich genomics-empowered (1) biodiversity monitoring, (2) population genetics, (3) pathogen and disease-vector genomic surveillance, (4) allergen and narcotic surveillance, (5) antimicrobial resistance surveillance and (6) bioprospecting.

Unraveling the genetic basis of Rhizobium rhizogenes-mediated transformation and hairy root formation in rose using a genome-wide association study

KEY MESSAGE

Multiple QTLs reveal the polygenic nature of R. rhizogenes-mediated transformation and hairy root formation in roses, with five key regions explaining 12.0-26.9% of trait variability and transformation-related candidate genes identified. Understanding genetic mechanisms of plant transformation remains crucial for biotechnology. This is particularly relevant for roses and other woody ornamentals that exhibit recalcitrant behavior in transformation procedures. Rhizobium rhizogenes-mediated transformation leading to hairy root (HR) formation provides an excellent model system to study transformation processes and host-pathogen interactions. Therefore, this study aimed to identify quantitative trait loci (QTLs) associated with HR formation and explore their relationship with adventitious root (AR) formation in rose as a model for woody ornamentals. A diversity panel of 104 in vitro grown rose genotypes was transformed with R. rhizogenes strain ATCC 15834 carrying a green fluorescent protein reporter gene. Phenotypic data on callus and root formation were collected for laminae and petioles. A genome-wide association study using 23,419 single-nucleotide polymorphism markers revealed significant QTLs on chromosomes one and two for root formation traits. Five key genomic regions explained 12.0-26.9% of trait variability, with some peaks overlapping previously reported QTLs for AR formation. This genetic overlap was supported by weak to moderate correlations between HR and AR formation traits, particularly in petioles. Candidate gene identification through literature review and transcriptomic data analysis revealed ten candidate genes involved in bacterial response, hormone signaling, and stress responses. Our findings provide new insights into the genetic control of HR formation in roses and highlight potential targets for improving transformation efficiency in ornamental crops, thereby facilitating future research and breeding applications.

Long-read genome assembly of the insect model organism Tribolium castaneum reveals spread of satellite DNA in gene-rich regions by recurrent burst events

Eukaryotic genomes are replete with satellite DNAs (satDNAs), large stretches of tandemly repeated sequences that are mostly underrepresented in genome assemblies. Here we combined nanopore long-read sequencing with a reference-guided assembly approach to generate an improved, high-quality genome assembly, TcasONT, of the model beetle Tribolium castaneum Enriched by 45 Mb in repetitive regions, the new assembly comprises almost the entire genome sequence. We use the enhanced assembly to conduct global and in-depth analyses of abundant euchromatic satDNAs. Unexpectedly, we show the extensive spread of satDNAs in gene-rich regions, including long arrays. The sequence similarity relationships between satDNA monomers and arrays indicate a recent exchange of satDNA arrays between different chromosomes. We propose a scenario of their genome dynamics characterized by repeated bursts of satDNAs spreading through euchromatin, followed by a process of elongation and homogenization of arrays. We find that suppressed recombination on the X Chromosome has no significant effect on the spread of satDNAs but the X rather tolerates the amplification of satDNAs into longer arrays. Analyses of arrays' neighboring regions show a tendency of one satDNA to be associated with transposable-like elements. Using 2D electrophoresis followed by Southern blotting, we prove Cast satDNAs' presence in the fraction of extrachromosomal circular DNA (eccDNA). We point to two mechanisms that enable this satDNA spread to occur: transposition by transposable elements and insertion mediated by eccDNA. The presence of such a large proportion of satDNA in gene-rich regions inevitably gives rise to speculation about their possible influence on gene expression.

PTMoreR-enabled cross-species PTM mapping and comparative phosphoproteomics across mammals

To support PTM proteomic analysis and annotation in different species, we developed PTMoreR, a user-friendly tool that considers the surrounding amino acid sequences of PTM sites during BLAST, enabling a motif-centric analysis across species. By controlling sequence window similarity, PTMoreR can map phosphoproteomic results between any two species, perform site-level functional enrichment analysis, and generate kinase-substrate networks. We demonstrate that the majority of real P-sites in mice can be inferred from experimentally derived human P-sites with PTMoreR mapping. Furthermore, the compositions of 129 mammalian phosphoproteomes can also be predicted using PTMoreR. The method also identifies cross-species phosphorylation events that occur on proteins with an increased tendency to respond to the environmental factors. Moreover, the classic kinase motifs can be extracted across mammalian species, offering an evolutionary angle for refining current motifs. PTMoreR supports PTM proteomics in non-human species and facilitates quantitative phosphoproteomic analysis.

The Genome Organization of 5S rRNA Genes in the Model Organism Tribolium castaneum and Its Sibling Species Tribolium freemani

5S ribosomal DNAs (rDNAs) are arranged in tandem and are often under-represented in genome assemblies. In the present study, we performed a global and in-depth analysis of the 5S rDNAs in the model insect Tribolium castaneum and its closely related species Tribolium freemani. To accomplish this goal, we used our recently published genome assemblies based on Nanopore and PacBio long-read sequencing. Although these closely related species share the 5S rRNA gene sequence with high homology, they show a different organization of the 5S rDNA locus. Analysis of 5S rDNA arrays in T. castaneum revealed a typical tandemly repeated organization characterized by repeat units consisting of the 121 bp long 5S rRNA gene and the 71 bp long nontranscribed spacer (NTS). In contrast, T. freemani showed a much more complex organization of 5S rDNA arrays characterized by two patterns. The first is based on the association of 5S rRNA gene with arrays of a satellite DNA, representing the NTS sequence of the 5S rDNA genes in T. freemani. The second, more complex type is characterized by a somewhat less frequent occurrence of the 5S rRNA gene and its association with longer satellite DNA arrays that are regularly interrupted by Jockey-like retrotransposons. This organization, in which the ribosomal gene is associated with two completely different repetitive elements such as satellite DNAs and retrotransposons, suggests that the 5S rRNA gene, regardless of its crucial function in the genome, could be a subject of extremely dynamic genomic rearrangements.

Populating a Continent: Phylogenomics Reveal the Timing of Australian Frog Diversification

The Australian continent's size and isolation make it an ideal place for studying the accumulation and evolution of biodiversity. Long separated from the ancient supercontinent Gondwana, most of Australia's plants and animals are unique and endemic, including the continent's frogs. Australian frogs comprise a remarkable ecological and morphological diversity categorized into a small number of distantly related radiations. We present a phylogenomic hypothesis based on an exon-capture dataset that spans the main clades of Australian myobatrachoid, pelodryadid hyloid, and microhylid frogs. Our time-calibrated phylogenomic-scale phylogeny identifies great disparity in the relative ages of these groups that vary from Gondwanan relics to recent immigrants from Asia and include arguably the continent's oldest living vertebrate radiation. This age stratification provides insight into the colonization of, and diversification on, the Australian continent through deep time, during periods of dramatic climatic and community changes. Contemporary Australian frog diversity highlights the adaptive capacity of anurans, particularly in response to heat and aridity, and explains why they are one of the continent's most visible faunas. [Anuran; adaptive radiation; Gondwana; phylogenetics].

Portable BLAST-like algorithm library and its implementations for command line, Python, and R

Basic local-alignment search tool (BLAST) is a versatile and commonly used sequence analysis tool in bioinformatics. BLAST permits fast and flexible sequence similarity searches across nucleotide and amino acid sequences, leading to diverse applications such as protein domain identification, orthology searches, and phylogenetic annotation. Most BLAST implementations are command line tools which produce output as comma-separated values files. However, a portable, modular and embeddable implementation of a BLAST-like algorithm, is still missing from our toolbox. Here we present nsearch, a command line tool and C++11 library which provides BLAST-like functionality that can easily be embedded in any application. As an example of this portability we present Blaster which leverages nsearch to provide native BLAST-like functionality for the R programming language, as well as npysearch which provides similar functionality for Python. These packages permit embedding BLAST-like functionality into larger frameworks such as Shiny or Django applications. Benchmarks show that nsearch, npysearch, and Blaster are comparable in speed and accuracy to other commonly used modern BLAST implementations such as VSEARCH and BLAST+. We envision similar implementations of nsearch for other languages commonly used in data science such as Julia to facilitate sequence similarity comparisons. Nsearch, Blaster and npysearch are free to use under the BSD 3.0 license and available on Github Conda, CRAN (Blaster) and PyPi (npysearch).

The Low-Copy-Number Satellite DNAs of the Model Beetle Tribolium castaneum

The red flour beetle Tribolium castaneum is an important pest of stored agricultural products and the first beetle whose genome was sequenced. So far, one high-copy-number and ten moderate-copy-number satellite DNAs (satDNAs) have been described in the assembled part of its genome. In this work, we aimed to catalog the entire collection of T. castaneum satDNAs. We resequenced the genome using Illumina technology and predicted potential satDNAs via graph-based sequence clustering. In this way, we discovered 46 novel satDNAs that occupied a total of 2.1% of the genome and were, therefore, considered low-copy-number satellites. Their repeat units, preferentially 140-180 bp and 300-340 bp long, showed a high A + T composition ranging from 59.2 to 80.1%. In the current assembly, we annotated the majority of the low-copy-number satDNAs on one or a few chromosomes, discovering mainly transposable elements in their vicinity. The current assembly also revealed that many of the in silico predicted satDNAs were organized into short arrays not much longer than five consecutive repeats, and some of them also had numerous repeat units scattered throughout the genome. Although 20% of the unassembled genome sequence masked the genuine state, the predominance of scattered repeats for some low-copy satDNAs raises the question of whether these are essentially interspersed repeats that occur in tandem only sporadically, with the potential to be satDNA "seeds".

Reference-Guided De Novo Genome Assembly of the Flour Beetle Tribolium freemani

The flour beetle Tribolium freemani is a sibling species of the model organism and important pest Tribolium castaneum. The two species are so closely related that they can produce hybrid progeny, but the genetic basis of their differences has not been revealed. In this work, we sequenced the T. freemani genome by applying PacBio HiFi technology. Using the well-assembled T. castaneum genome as a reference, we assembled 262 Mb of the T. freemani genomic sequence and anchored it in 10 linkage groups corresponding to nine autosomes and sex chromosome X. The assembly showed 99.8% completeness of conserved insect genes, indicating a high-quality reference genome. Comparison with the T. castaneum assembly revealed that the main differences in genomic sequence between the two sibling species come from repetitive DNA, including interspersed and tandem repeats. In this work, we also provided the complete assembled mitochondrial genome of T. freemani. Although the genome assembly needs to be ameliorated in tandemly repeated regions, the first version of the T. freemani reference genome and the complete mitogenome presented here represent useful resources for comparative evolutionary studies of related species and for further basic and applied research on different biological aspects of economically important pests.