Isolation and characterization of reverse transcriptase fragments of LTR retrotransposons from the genome of Chenopodium quinoa (Amaranthaceae)

A telomere-to-telomere gap-free assembly integrating multi-omics uncovers the genetic mechanism of fruit quality and important agronomic trait associations in pomegranate

Pomegranate is an important perennial fruit tree distributed worldwide. Reference genomes with gaps and limit gene identification controlling important agronomic traits hinder its functional genomics and genetic improvements. Here, we reported a telomere-to-telomere (T2T) gap-free genome assembly of the distinctive cultivar 'Moshiliu'. The Moshiliu reference genome was assembled into eight chromosomes without gaps, totalling ~366.71 Mb, with 32 158 predicted protein-coding genes. All 16 telomeres and eight centromeres were characterized; combined with FISH analysis, we revealed the atypical telomere units in pomegranate as TTTTAGGG. Furthermore, a total of 16 loci associated with 15 important agronomic traits were identified based on GWAS of 146 accessions. Gene editing and biochemical experiments demonstrated that a 37.2-Kb unique chromosome translocation disrupting the coding domain sequence of PgANS was responsible for anthocyanin-less, knockout of PgANS in pomegranate exhibited a defect in anthocyanin production; a unique repeat expansion in the promoter of PgANR may affected its expression, resulting in black peel; notably, the G → A transversion located at the 166-bp coding domain of PgNST3, which caused a E56K mutation in the PgNST3 protein, closely linked with soft-seed trait. Overexpression of PgNST3A in tomato presented smaller and softer seed coats. The E56K mutation in PgNST3 protein, eliminated the binding ability of PgNST3 to the PgMYB46 promoter, which subsequently affected the thickness of the inner seed coat of soft-seeded pomegranates. Collectively, the validated gap-free genome, the identified genes controlling important traits and the CRISPR-Cas9-mediated gene knockout system all provided invaluable resources for pomegranate precise breeding.

In silico PCR analysis: a comprehensive bioinformatics tool for enhancing nucleic acid amplification assays

Nucleic acid amplification assays represent a pivotal category of methodologies for targeted sequence detection within contemporary biological research, boasting diverse utility in diagnostics, identification, and DNA sequencing. The foundational principles of these assays have been extrapolated to various simple and intricate nucleic acid amplification technologies. Concurrently, a burgeoning trend toward computational or virtual methodologies is exemplified by in silico PCR analysis. In silico PCR analysis is a valuable and productive adjunctive approach for ensuring primer or probe specificity across a broad spectrum of PCR applications encompassing gene discovery through homology analysis, molecular diagnostics, DNA profiling, and repeat sequence identification. The prediction of primer and probe sensitivity and specificity necessitates thorough database searches, accounting for an optimal balance of mismatch tolerance, sequence similarity, and thermal stability. This software facilitates in silico PCR analyses of both linear and circular DNA templates, including bisulfited treatment DNA, enabling multiple primer or probe searches within databases of varying scales alongside advanced search functionalities. This tool is suitable for processing batch files and is essential for automation when working with large amounts of data.

Chromosomal-level genome and multi-omics dataset provides new insights into leaf pigmentation in Acer palmatum

Acer palmatum (A. palmatum), a deciduous shrub or small arbour which belongs to Acer of Aceraceae, is an excellent greening species as well as a beautiful ornamental plant. In this study, a high-quality chromosome-level reference genome for A. palmatum was constructed using Oxford Nanopore sequencing and Hi-C technology. The assembly genome was ∼745.78 Mb long with a contig N50 length of 3.20 Mb, and 95.30 % (710.71 Mb) of the assembly was anchored into 13 pseudochromosomes. A total of 28,559 protein-coding genes were obtained, ∼90.02 % (25,710) of which could be functionally annotated. The genomic evolutionary analysis revealed that A. palmatum is most closely related to A. yangbiense and A. truncatum, and underwent only an ancient gamma whole-genome duplication event. Despite lacking a recent independent WGD, 25,795 (90.32 %) genes of A. palmatum were duplicated, and the unique/expanded gene families were linked with genes involved in plant-pathogen interaction and several metabolic pathways, which might underpin adaptability. A combined genomic, transcriptomic, and metabolomic analysis related to the biosynthesis of anthocyanin in leaves during the different season were characterized. The results indicate that the dark-purple colouration of the leaves in spring was caused by a high amount of anthocyanins, especially delphinidin and its derivatives; and the red colouration of the leaves in autumn by a high amount of cyanidin 3-O-glucoside. In conclusion, these valuable multi-omic resources offer important foundations to explore the molecular regulation mechanism in leaf colouration and also provide a platform for the scientific and efficient utilization of A. palmatum.

Chromosome-Level Genome Assembly of Bupleurum chinense DC Provides Insights Into the Saikosaponin Biosynthesis

Bupleurum chinense DC is a plant widely used in Chinese traditional medicine. Saikosaponins are the major bioactive constituents of B. chinense DC. Saikosaponins biosynthesis in Bupleurum has been more intensively studied than any other metabolic processes or bioactive constituents. However, whole-genome sequencing and chromosome-level assembly for Bupleurum genus have not been reported yet. Here, we report a high-quality chromosome-level genome of B. chinense DC. through the integration of PacBio long-read sequencing, Illumina short-read sequencing, and Hi-C sequencing. The genome was phased into haplotype 0 (621.27 Mb with a contig N50 of 16.86 Mb and a scaffold N50 of 92.25 Mb) and haplotype 1 (600.48 Mb with a contig N50 of 23.90 Mb and a scaffold N50 of 102.68 Mb). A total of 45,909 and 35,805 protein-coding genes were predicted in haplotypes 0 and 1, respectively. The enrichment analyses suggested that the gene families that expanded during the evolution of B. chinense DC are involved in the biosynthesis of isoquinoline alkaloid, tyrosine, and anthocyanin. Furthermore, we analyzed the genes involved in saikosaponin biosynthesis and determined the candidate P450 and UGT genes in the third stage of saikosaponins biosynthetic, which provided new insight into the saikosaponins biosynthetic. The genomic data provide a valuable resource for future investigations of the molecular mechanisms, biological functions, and evolutionary adaptations of B. chinense DC.

Satellite DNA landscapes after allotetraploidization of quinoa (Chenopodium quinoa) reveal unique A and B subgenomes

If two related plant species hybridize, their genomes may be combined and duplicated within a single nucleus, thereby forming an allotetraploid. How the emerging plant balances two co-evolved genomes is still a matter of ongoing research. Here, we focus on satellite DNA (satDNA), the fastest turn-over sequence class in eukaryotes, aiming to trace its emergence, amplification, and loss during plant speciation and allopolyploidization. As a model, we used Chenopodium quinoa Willd. (quinoa), an allopolyploid crop with 2n = 4x = 36 chromosomes. Quinoa originated by hybridization of an unknown female American Chenopodium diploid (AA genome) with an unknown male Old World diploid species (BB genome), dating back 3.3-6.3 million years. Applying short read clustering to quinoa (AABB), C. pallidicaule (AA), and C. suecicum (BB) whole genome shotgun sequences, we classified their repetitive fractions, and identified and characterized seven satDNA families, together with the 5S rDNA model repeat. We show unequal satDNA amplification (two families) and exclusive occurrence (four families) in the AA and BB diploids by read mapping as well as Southern, genomic, and fluorescent in situ hybridization. Whereas the satDNA distributions support C. suecicum as possible parental species, we were able to exclude C. pallidicaule as progenitor due to unique repeat profiles. Using quinoa long reads and scaffolds, we detected only limited evidence of intergenomic homogenization of satDNA after allopolyploidization, but were able to exclude dispersal of 5S rRNA genes between subgenomes. Our results exemplify the complex route of tandem repeat evolution through Chenopodium speciation and allopolyploidization, and may provide sequence targets for the identification of quinoa's progenitors.

The population genetic structure approach adds new insights into the evolution of plant LTR retrotransposon lineages

Long terminal repeat retrotransposons (LTR-RTs) in plant genomes differ in abundance, structure and genomic distribution, reflecting the large number of evolutionary lineages. Elements within lineages can be considered populations, in which each element is an individual in its genomic environment. In this way, it would be reasonable to apply microevolutionary analyses to understand transposable element (TE) evolution, such as those used to study the genetic structure of natural populations. Here, we applied a Bayesian method to infer genetic structure of populations together with classical phylogenetic and dating tools to analyze LTR-RT evolution using the monocot Setaria italica as a model species. In contrast to a phylogeny, the Bayesian clusterization method identifies populations by assigning individuals to one or more clusters according to the most probabilistic scenario of admixture, based on genetic diversity patterns. In this work, each LTR-RT insertion was considered to be one individual and each LTR-RT lineage was considered to be a single species. Nine evolutionary lineages of LTR-RTs were identified in the S. italica genome that had different genetic structures with variable numbers of clusters and levels of admixture. Comprehensive analysis of the phylogenetic, clusterization and time of insertion data allowed us to hypothesize that admixed elements represent sequences that harbor ancestral polymorphic sequence signatures. In conclusion, application of microevolutionary concepts in genome evolution studies is suitable as a complementary approach to phylogenetic analyses to address the evolutionary history and functional features of TEs.

Organization and evolution of two repetitive sequences, 18-24J and 12-13P, in the genome of Chenopodium (Amaranthaceae)

The abundance and chromosomal organization of two repetitive sequences named 12-13P and 18-24J were analyzed in 24 diploid and nine polyploid species of Chenopodium s.l., with special attention to Chenopodium s.s. Both sequences were predominantly present in species of Chenopodium s.s.; however, differences in the amplification levels were observed among the species. The 12-13P repeat was highly amplified in all of the analyzed Eurasian species, whereas the American diploids showed a marked variation in the amplification levels. The 12-13P repeat contains a tandemly arranged 40 bp minisatellite element forming a large proportion of the genome of Chenopodium (up to 3.5%). FISH revealed its localization to the pericentromeric regions of the chromosomes. The chromosomal distribution of 12-13P delivered additional chromosomal marker for B-genome diploids. The 18-24J repeat showed a dispersed organization in all of the chromosomes of the analyzed diploid species and the Eurasian tetraploids. In the American allotetraploids (C. quinoa, C. berlandieri) and Eurasian allohexaploids (e.g., C. album) very intense hybridization signals of 18-24J were observed only on 18 chromosomes that belong to the B subgenome of these polyploids. Combined cytogenetic and molecular analyses suggests that reorganization of these two repeats accompanied the diversification and speciation of diploid (especially A genome) and polyploid species of Chenopodium s.s.

Reverse transcriptase sequences from mulberry LTR retrotransposons: characterization analysis

Copia and Gypsy play important roles in structural, functional and evolutionary dynamics of plant genomes. In this study, a total of 106 and 101, Copia and Gypsy reverse transcriptase (rt) were amplified respectively in the Morus notabilis genome using degenerate primers. All sequences exhibited high levels of heterogeneity, were rich in AT and possessed higher sequence divergence of Copia rt in comparison to Gypsy rt. Two reasons are likely to account for this phenomenon: a) these elements often experience deletions or fragmentation by illegitimate or unequal homologous recombination in the transposition process; b) strong purifying selective pressure drives the evolution of these elements through “selective silencing” with random mutation and eventual deletion from the host genome. Interestingly, mulberry rt clustered with other rt from distantly related taxa according to the phylogenetic analysis. This phenomenon did not result from horizontal transposable element transfer. Results obtained from fluorescence in situ hybridization revealed that most of the hybridization signals were preferentially concentrated in pericentromeric and distal regions of chromosomes, and these elements may play important roles in the regions in which they are found. Results of this study support the continued pursuit of further functional studies of Copia and Gypsy in the mulberry genome.

Characterization, Genomic Organization, Abundance, and Chromosomal Distribution of Ty1-copia Retrotransposons in Erianthus arundinaceus

Erianthus arundinaceus is an important wild species of the genus Saccharum with many valuable traits. However, the composition and structure of its genome are largely unknown, which have hindered its utilization in sugarcane breeding and evolutionary research. Retrotransposons constitute an appreciable fraction of plant genomes and may have played a significant role in the evolution and sequence organization of genomes. In the current study, we investigate the phylogenetic diversity and genomic abundance of Ty1-copia retrotransposons for the first time and inspect their chromosomal distribution patterns in E. arundinaceus. In total, 70 Ty1-copia reverse transcriptase (RT) sequences with significant levels of heterogeneity were obtained. The phylogenetic analysis revealed these Ty1-copia retrotransposons were classified into four distinct evolutionary lineages (Tork/TAR, Tork/Angela, Retrofit/Ale, and Sire/Maximus). Dot-blot analysis showed estimated the total copy number of Ty1-copia retrotransposons to be about 4.5 × 10³ in the E. arundinaceus genome, indicating they were a significant component. Fluorescence in situ hybridization revealed that Ty1-copia retrotransposons from the four lineages had strikingly similar patterns of chromosomal enrichment, being exclusively enriched in the subterminal heterochromatic regions of most E. arundinaceus chromosomes. This is the first clear evidence of the presence of Ty1-copia retrotransposons in the subterminal heterochromatin of E. arundinaceus. Altogether, these results promote the understanding of the diversification of Ty1-copia retrotransposons and shed light on their chromosomal distribution patterns in E. arundinaceus.

Diversification, evolution and methylation of short interspersed nuclear element families in sugar beet and related Amaranthaceae species

Short interspersed nuclear elements (SINEs) are non-autonomous non-long terminal repeat retrotransposons which are widely distributed in eukaryotic organisms. While SINEs have been intensively studied in animals, only limited information is available about plant SINEs. We analysed 22 SINE families from seven genomes of the Amaranthaceae family and identified 34 806 SINEs, including 19 549 full-length copies. With the focus on sugar beet (Beta vulgaris), we performed a comparative analysis of the diversity, genomic and chromosomal organization and the methylation of SINEs to provide a detailed insight into the evolution and age of Amaranthaceae SINEs. The lengths of consensus sequences of SINEs range from 113 nucleotides (nt) up to 224 nt. The SINEs show dispersed distribution on all chromosomes but were found with higher incidence in subterminal euchromatic chromosome regions. The methylation of SINEs is increased compared with their flanking regions, and the strongest effect is visible for cytosines in the CHH context, indicating an involvement of asymmetric methylation in the silencing of SINEs.

Characterization and chromosomal organization of Ty1-copia retrotransposons in wax gourd

Wax gourd (2n=2x=24) is an important vegetable species in Cucurbitaceae. Because it can be stored for a very long period of time, it plays an important role in ensuring the annual supply and regulating off-season supply of the vegetables. However, the availability of genetic information about wax gourd is limited. This study aimed to identify the useful genetic information for wax gourd. The conserved domains of reverse transcriptase (RT) genes of Ty1-copia retrotransposons were isolated from the genome of wax gourd using degenerate oligonucleotide primers. A total of twenty eight RT sequences were obtained, which showed high heterogeneity with the similarity ranging from 47.5% to 94.3%. Sixteen (57.1%) of them were found to be defective, being disrupted by stop codons and/or frameshift mutations. These 28 sequences were divided into five subfamilies. The comparative phylogenetic analysis with other Cucurbitaceae species from GenBank database showed that most retrotransposons derived from the same genus tended to cluster together, although there were a few exceptions. These results indicate that both vertical transmission and horizontal transmission are the sources of Ty1-copia retrotransposons in wax gourd. Fluorescent in situ hybridization (FISH) with Ty1-copia retrotransposon sequences as probes revealed that this kind of retrotransposons had a dispersed genomic organization, physically distributed among all the chromosomes of wax gourd, with clusters in the heterochromatin regions. This is the first report of Ty1-copia retrotransposons in wax gourd, which would be helpful for our understanding about the organization and evolutions of wax gourd genome and also provide valuable information for our utilization of wax gourd retrotransposons.