Terminal Repeat Research Articles

Transposable elements in Rosaceae: insights into genome evolution, expression dynamics, and syntenic gene regulation.

Transposable elements (TEs) exert significant influence on plant genomic structure and gene expression. Here, we explored TE-related aspects across 14 Rosaceae genomes, investigating genomic distribution, transposition activity, expression patterns, and nearby differentially expressed genes (DEGs). Analyses unveiled distinct long terminal repeat retrotransposon (LTR-RT) evolutionary patterns, reflecting varied genome size changes among nine species over the past million years. In the past 2.5 million years, Rubus idaeus showed a transposition rate twice as fast as Fragaria vesca, while Pyrus bretschneideri displayed significantly faster transposition compared with Crataegus pinnatifida. Genes adjacent to recent TE insertions were linked to adversity resistance, while those near previous insertions were functionally enriched in morphogenesis, enzyme activity, and metabolic processes. Expression analysis revealed diverse responses of LTR-RTs to internal or external conditions. Furthermore, we identified 3695 pairs of syntenic DEGs proximal to TEs in Malus domestica cv. 'Gala' and M. domestica (GDDH13), suggesting TE insertions may contribute to varietal trait differences in these apple varieties. Our study across representative Rosaceae species underscores the pivotal role of TEs in plant genome evolution within this diverse family. It elucidates how these elements regulate syntenic DEGs on a genome-wide scale, offering insights into Rosaceae-specific genomic evolution.

Pangenome analysis reveals transposon-driven genome evolution in cotton

BackgroundTransposable elements (TEs) have a profound influence on the trajectory of plant evolution, driving genome expansion and catalyzing phenotypic diversification. The pangenome, a comprehensive genetic pool encompassing all variations within a species, serves as an invaluable tool, unaffected by the confounding factors of intraspecific diversity. This allows for a more nuanced exploration of plant TE evolution.ResultsHere, we constructed a pangenome for diploid A-genome cotton using 344 accessions from representative geographical regions, including 223 from China as the main component. We found 511 Mb of non-reference sequences (NRSs) and revealed the presence of 5479 previously undiscovered protein-coding genes. Our comprehensive approach enabled us to decipher the genetic underpinnings of the distinct geographic distributions of cotton. Notably, we identified 3301 presence-absence variations (PAVs) that are closely tied to gene expression patterns within the pangenome, among which 2342 novel expression quantitative trait loci (eQTLs) were found residing in NRSs. Our investigation also unveiled contrasting patterns of transposon proliferation between diploid and tetraploid cotton, with long terminal repeat (LTR) retrotransposons exhibiting a synchronized surge in polyploids. Furthermore, the invasion of LTR retrotransposons from the A subgenome to the D subgenome triggered a substantial expansion of the latter following polyploidization. In addition, we found that TE insertions were responsible for the loss of 36.2% of species-specific genes, as well as the generation of entirely new species-specific genes.ConclusionsOur pangenome analyses provide new insights into cotton genomics and subgenome dynamics after polyploidization and demonstrate the power of pangenome approaches for elucidating transposon impacts and genome evolution.

Celine, a long interspersed nuclear element retrotransposon, colonizes in the centromeres of poplar chromosomes.

Centromeres in most multicellular eukaryotes are composed of long arrays of repetitive DNA sequences. Interestingly, several transposable elements, including the well-known long terminal repeat centromeric retrotransposon of maize (CRM), were found to be enriched in functional centromeres marked by the centromeric histone H3 (CENH3). Here, we report a centromeric long interspersed nuclear element (LINE), Celine, in Populus species. Celine has colonized preferentially in the CENH3-associated chromatin of every poplar chromosome, with 84% of the Celine elements localized in the CENH3-binding domains. In contrast, only 51% of the CRM elements were bound to CENH3 domains in Populus trichocarpa. These results suggest different centromere targeting mechanisms employed by Celine and CRM elements. Nevertheless, the high target specificity seems to be detrimental to further amplification of the Celine elements, leading to a shorter life span and patchy distribution among plant species compared with the CRM elements. Using a phylogenetically guided approach, we were able to identify Celine-like LINE elements in tea plant (Camellia sinensis) and green ash tree (Fraxinus pennsylvanica). The centromeric localization of these Celine-like LINEs was confirmed in both species. We demonstrate that the centromere targeting property of Celine-like LINEs is of primitive origin and has been conserved among distantly related plant species.

Enhanced transgene expression from single-stranded AAV vectors in human cells in vitro and in murine hepatocytes in vivo

We identified that distal 10 nucleotides in the D-sequence in AAV2 inverted terminal repeat (ITR) share partial sequence homology to ½ binding site of glucocorticoid receptor-binding element (GRE). Here, we describe that (i) purified GR binds to AAV2 D-sequence, and the D-sequence competes with GR binding to its cognate binding site; (ii) dexamethasone-mediated activation of GR pathway significantly increases the transduction efficiency of AAV2 vectors in human cells; (iii) human osteosarcoma cells, U2OS, which lack expression of GR, are poorly transduced by AAV2 vectors, but stable transfection with a GR expression plasmid restores vector-mediated transgene expression; (iv) replacement of the distal 10 nucleotides in the D-sequence of the AAV2 ITR with a full-length GRE consensus sequence significantly enhances transgene expression in human cells in vitro and in murine hepatocytes in vivo, and (v) none of the ITRs in AAV1, AAV3, AAV4, AAV5, and AAV6 genomes contains the GRE ½ binding site, and insertion of with a full-length GRE consensus sequence in the AAV6-ITR also significantly enhances transgene expression from AAV6 vectors, both in vitro and in vivo. These novel vectors, termed generation Y (GenY) AAV vectors, which are serotype-, transgene-, or promoter-agnostic, should be useful in human gene therapy.

A PAX3 insertion in the Celestial breed and certain feline breeding lines with dominant blue eyes.

During the last 60 years many inherited traits in domestic outbred cats were selected and retained giving birth to new breeds characterised by singular coat or morphological phenotypes. Among them, minimal white spotting associated with blue eyes was selected by feline breeders to create the Altai, Topaz, and Celestial breeds. Various established breeds also introduced this trait in their lineages. The trait, that was confirmed as autosomal dominant by breeding data, was first described in domestic cats from Kazakhstan and Russia, in British shorthair and British longhair from Russia, and in Maine Coon cats from the Netherlands, suggesting different founding effects. Using a genome-wide association study we identified a single region on chromosome C1 that was associated with the minimal white spotting and blue eyes phenotype (also called DBE by breeders for dominant blue eyes) in the French Celestial breed. Within that region we identified Paired Box3 (PAX3) as the strongest candidate gene, since PAX3 is a key regulator of MITF (Melanocyte-Inducing Transcription Factor) and PAX3 variants have been previously identified in various species showing white spotting with or without blue eyes including the mouse and the horse. Whole genome sequencing of a Celestial cat revealed an endogenous retrovirus LTR (long terminal repeat) insertion within PAX3 intron 4 known to contain regulatory sequences (conserved non-coding element [CNE]) involved in PAX3 expression. The insertion is in the vicinity of CNE2 and CNE3. All 52 Celestial and Celestial-mixed cats with a DBE phenotype presented the insertion, that was absent in their 22 non-DBE littermates and in 87 non-DBE cats from various breeds. The outbred Celestial founder was also heterozygous for the insertion. Additionally, the variant was found in nine DBE Maine Coon cats related to the Celestial founder and four DBE Siberian cats with an uncertain origin. Segregation of the variant in the Celestial breed is consistent with dominant inheritance and does not appear to be associated with deafness. We propose that this NC_018730.3:g.206974029_206974030insN[395] variant represents the DBECEL (Celestial Dominant Blue Eyes) allele in the domestic cat.

Look4LTRs: a Long terminal repeat retrotransposon detection tool capable of cross species studies and discovering recently nested repeats

Plant genomes include large numbers of transposable elements. One particular type of these elements is flanked by two Long Terminal Repeats (LTRs) and can translocate using RNA. Such elements are known as LTR-retrotransposons; they are the most abundant type of transposons in plant genomes. They have many important functions involving gene regulation and the rise of new genes and pseudo genes in response to severe stress. Additionally, LTR-retrotransposons have several applications in biotechnology. Due to the abundance and the importance of LTR-retrotransposons, multiple computational tools have been developed for their detection. However, none of these tools take advantages of the availability of related genomes; they process one chromosome at a time. Further, recently nested LTR-retrotransposons (multiple elements of the same family are inserted into each other) cannot be annotated accurately — or cannot be annotated at all — by the currently available tools. Motivated to overcome these two limitations, we built Look4LTRs, which can annotate LTR-retrotransposons in multiple related genomes simultaneously and discover recently nested elements. The methodology of Look4LTRs depends on techniques imported from the signal-processing field, graph algorithms, and machine learning with a minimal use of alignment algorithms. Four plant genomes were used in developing Look4LTRs and eight plant genomes for evaluating it in contrast to three related tools. Look4LTRs is the fastest while maintaining better or comparable F1 scores (the harmonic average of recall and precision) to those obtained by the other tools. Our results demonstrate the added benefit of annotating LTR-retrotransposons in multiple related genomes simultaneously and the ability to discover recently nested elements. Expert human manual examination of six elements — not included in the ground truth — revealed that three elements belong to known families and two elements are likely from new families. With respect to examining recently nested LTR-retrotransposons, three out of five were confirmed to be valid elements. Look4LTRs — with its speed, accuracy, and novel features — represents a true advancement in the annotation of LTR-retrotransposons, opening the door to many studies focused on understanding their functions in plants.

PCBP1 interacts with the HTLV-1 Tax oncoprotein to potentiate NF-κB activation.

Human T-cell leukemia virus type 1 (HTLV-1) is the etiological agent of adult T-cell leukemia/lymphoma. The HTLV-1 Tax constitutively activates nuclear factor-κB (NF-κB) to promote the survival and transformation of HTLV-1-infected T cells. Despite extensive study of Tax, how Tax interacts with host factors to regulate NF-κB activation and HTLV-1-driven cell proliferation is not entirely clear. Here, we showed that overexpression of Poly (rC)-binding protein 1 (PCBP1) promoted Tax-mediated IκB kinase (IKK)-NF-κB signaling activation, whereas knockdown of PCBP1 attenuated Tax-dependent IKK-NF-κB activation. However, Tax activation of HTLV-1 long terminal repeat was unaffected by PCBP1. Furthermore, depletion of PCBP1 led to apoptosis and reduced proliferation of HTLV-1-transformed cells. Mechanistically, PCBP1 interacted and co-localized with Tax in the cytoplasm, and PCBP1 KH3 domain was indispensable for the interaction between PCBP1 and Tax. Moreover, PCBP1 facilitated the assembly of Tax/IKK complex. Collectively, our results demonstrated that PCBP1 may exert an essential effect in Tax/IKK complex combination and subsequent NF-κB activation, which provides a novel insight into the pathogenetic mechanisms of HTLV-1.

The identification of retro-DNAs in primate genomes as DNA transposons mobilizing via retrotransposition

Background Mobile elements (MEs) constitute a major portion of the genome in primates and other higher eukaryotes, and they play important role in genome evolution and gene function. MEs can be divided into two fundamentally different classes: DNA transposons which transpose in the genome in a “cut-and-paste” style, and retrotransposons which propagate in a “copy-and-paste” fashion via a process involving transcription and reverse-transcription. In primate genomes, DNA transposons are mostly dead, while many retrotransposons are still highly active. We report here the identification of a new type of MEs, which we call “retro-DNAs”, for their combined characteristics of these two fundamentally different ME classes. Methods A comparative computational genomic approach was used to analyze the reference genome sequences of 10 primate species consisting of five apes, four monkeys, and marmoset. Results From our analysis, we identified a total of 1,750 retro-DNAs, representing 748 unique insertion events in the genomes of ten primate species including human. These retro-DNAs contain sequences of DNA transposons but lack the terminal inverted repeats (TIRs), the hallmark of DNA transposons. Instead, they show characteristics of retrotransposons, such as polyA tails, longer target-site duplications (TSDs), and the “TT/AAAA” insertion site motif, suggesting the use of the L1-based target-primed reverse transcription (TPRT) mechanism. At least 40% of these retro-DNAs locate into genic regions, presenting potentials for impacting gene function. More interestingly, some retro-DNAs, as well as their parent sites, show certain levels of expression, suggesting that they have the potential to create more retro-DNA copies in the present primate genomes. Conclusions Although small in number, the identification of these retro-DNAs reveals a new mechanism for propagating DNA transposons in primate genomes without active canonical DNA transposon activity. Our data also suggest that the TPRT machinery may transpose a wider variety of DNA sequences in the genomes.

Fine Mapping and Candidate Gene Validation of Tomato Gene Carpelloid Stamen and Parthenocarpy (CSP)

Parthenocarpy and male sterility are highly desirable traits in tomato breeding and molecular study. The stamen carpelloid mutant generally displays male sterility. A natural mutant displaying carpelloid stamen and parthenocarpy (csp) was identified in our research group. In this study, the csp locus was finely mapped to a 65 kb interval, which contained six putative genes. One of them, Solyc04g081000, encodes the tomato class B MADS box gene TAP3 (syn. SlDEF). Sequencing data revealed that a copia long terminal repeat retrotransposon was inserted in the first intron of the TAP3 gene of the csp mutant. qRT-PCR showed that the expression of TAP3 was significantly down-regulated in the petals and stamens of the csp mutant. A phenotypic analysis of the TAP3 gene-edited mutants and allelism tests indicated that TAP3 was the gene underlying csp, and csp was a novel allelic mutation of TAP3. The results of this study will lay the foundation for a further analysis of the function of TAP3 and provide materials and a basis for a further study of the functional differentiation of tomato B-class genes.

Adaptive expansion of ERVK solo-LTRs is associated with Passeriformes speciation events

Endogenous retroviruses (ERVs) are ancient retroviral remnants integrated in host genomes, and commonly deleted through unequal homologous recombination, leaving solitary long terminal repeats (solo-LTRs). This study, analysing the genomes of 362 bird species and their reptilian and mammalian outgroups, reveals an unusually higher level of solo-LTRs formation in birds, indicating evolutionary forces might have purged ERVs during evolution. Strikingly in the order Passeriformes, and especially the parvorder Passerida, endogenous retrovirus K (ERVK) solo-LTRs showed bursts of formation and recurrent accumulations coinciding with speciation events over past 22 million years. Moreover, our results indicate that the ongoing expansion of ERVK solo-LTRs in these bird species, marked by high transcriptional activity of ERVK retroviral genes in reproductive organs, caused variation of solo-LTRs between individual zebra finches. We experimentally demonstrated that cis-regulatory activity of recently evolved ERVK solo-LTRs may significantly increase the expression level of ITGA2 in the brain of zebra finches compared to chickens. These findings suggest that ERVK solo-LTRs expansion may introduce novel genomic sequences acting as cis-regulatory elements and contribute to adaptive evolution. Overall, our results underscore that the residual sequences of ancient retroviruses could influence the adaptive diversification of species by regulating host gene expression.

A chromosome-level genome reveals genome evolution and molecular basis of anthraquinone biosynthesis in Rheum palmatum

BackgroundRhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism.ResultsWe assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening.ConclusionOverall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies.

Genomic, transcriptomic, and metabolomic analyses provide insights into the evolution and development of a medicinal plant Saposhnikovia divaricata (Apiaceae).

Saposhnikovia divaricata, 2n = 2x = 16, as a perennial species, is widely distributed in China, Mongolia, Russia, etc. It is a traditional Chinese herb used to treat tetanus, rubella pruritus, rheumatic arthralgia, and other diseases. Here, we assembled a 2.07Gb and N50 scaffold length of 227.67Mb high-quality chromosome-level genome of S. divaricata based on the PacBio Sequel II sequencing platform. The total number of genes identified was 42 948, and 42 456 of them were functionally annotated. A total of 85.07% of the genome was composed of repeat sequences, comprised mainly of long terminal repeats (LTRs) which represented 73.7% of the genome sequence. The genome size may have been affected by a recent whole-genome duplication event. Transcriptional and metabolic analyses revealed bolting and non-bolting S. divaricata differed in flavonoids, plant hormones, and some pharmacologically active components. The analysis of its genome, transcriptome, and metabolome helped to provide insights into the evolution of bolting and non-bolting phenotypes in wild and cultivated S. divaricata and lays the basis for genetic improvement of the species.

Chromosome-level genome assembly and characterization of the Calophaca sinica genome.

Calophaca sinica is a rare plant endemic to northern China which belongs to the Fabaceae family and possesses rich nutritional value. To support the preservation of the genetic resources of this plant, we have successfully generated a high-quality genome of C. sinica (1.06 Gb). Notably, transposable elements (TEs) constituted ~73% of the genome, with long terminal repeat retrotransposons (LTR-RTs) dominating this group of elements (~54% of the genome). The average intron length of the C. sinica genome was noticeably longer than what has been observed for closely related species. The expansion of LTR-RTs and elongated introns emerged had the largest influence on the enlarged genome size of C. sinica in comparison to other Fabaceae species. The proliferation of TEs could be explained by certain modes of gene duplication, namely, whole genome duplication (WGD) and dispersed duplication (DSD). Gene family expansion, which was found to enhance genes associated with metabolism, genetic maintenance, and environmental stress resistance, was a result of transposed duplicated genes (TRD) and WGD. The presented genomic analysis sheds light on the genetic architecture of C. sinica, as well as provides a starting point for future evolutionary biology, ecology, and functional genomics studies centred around C. sinica and closely related species.

Molecular characterization of Marek's Disease virus reveals reticuloendotheliosis virus-long terminal repeat integration in the genome of the field isolates in Egypt

The highly contagious, immunosuppressive, and cancer-causing Marek's disease virus (MDV) infects chickens. The financial costs of Marek's disease (MD) are significant for the chicken industry. In this study, a total of 180 samples from chicken farms suspected to be MDV-infected were collected. The chickens were sampled during the period between the months of October 2016 and February 2018 at Dakahlia and Damietta Governorates, Egypt. A total of 36 pooled samples were created. The prepared samples were inoculated into embryonated chicken eggs (ECEs). Indirect fluorescent antibody technique (IFAT) and ICP4 gene-based polymerase chain reaction (PCR) were used for MDV identification. For the genetic characterization of the identified virus, The ICP4 gene sequence was identified and compared with the sequences available from various regions of the world. Furthermore, the genomes of all detected MDVs were screened for the long terminal repeat (LTR) region of reticuloendotheliosis (REV) in their genomes. The results showed that 31 out of 36 pooled samples (86.1%) inoculated into ECEs displayed the characteristic pock lesions. By using IFAT and PCR to identify MDV in ECEs, positive results were found in 27 samples (75%). The Egyptian virus is thought to be genetically closely related to MDVs circulating in Ethiopia, China, and India. REV-LTR was amplified from 6 out of 27 field isolates genomes (22.2 %) while MDV vaccine strains were free from REV-LTR insertion. The integrated REV-LTRs depicted a close genetic relationship with those integrated in fowl poxvirus (FWPV) circulating in Egypt as well as those integrated in FWPVs and MDVs from China, USA, South Africa, and Australia. To the best of our knowledge, this investigation represents the first identification and characterization of REV-LTR insertions in Egyptian MDV field isolates. Given the findings above, additional research in the future seems crucial to determine how the REV-LTR insertions affect MDV pathogenesis, virulence, and insufficient vaccination protection.

Identification and characterization of endogenous retroviruses upon SARS-CoV-2 infection.

Endogenous retroviruses (ERVs) derived from the long terminal repeat (LTR) family of transposons constitute a significant portion of the mammalian genome, with origins tracing back to ancient viral infections. Despite comprising approximately 8% of the human genome, the specific role of ERVs in the pathogenesis of COVID-19 remains unclear. In this study, we conducted a genome-wide identification of ERVs in human peripheral blood mononuclear cells (hPBMCs) and primary lung epithelial cells from monkeys and mice, both infected and uninfected with SARS-CoV-2. We identified 405, 283, and 206 significantly up-regulated transposable elements (TEs) in hPBMCs, monkeys, and mice, respectively. This included 254, 119, 68, and 28 ERVs found in hPBMCs from severe and mild COVID-19 patients, monkeys, and transgenic mice expressing the human ACE2 receptor (hACE2) and infected with SARS-CoV-2. Furthermore, analysis using the Genomic Regions Enrichment of Annotations Tool (GREAT) revealed certain parental genomic sequences of these up-regulated ERVs in COVID-19 patients may be involved in various biological processes, including histone modification and viral replication. Of particular interest, we identified 210 ERVs specifically up-regulated in the severe COVID-19 group. The genes associated with these differentially expressed ERVs were enriched in processes such as immune response activation and histone modification. HERV1_I-int: ERV1:LTR and LTR7Y: ERV1:LTR were highlighted as potential biomarkers for evaluating the severity of COVID-19. Additionally, validation of our findings using RT-qPCR in Bone Marrow-Derived Macrophages (BMDMs) from mice infected by HSV-1 and VSV provided further support to our results. This study offers insights into the expression patterns and potential roles of ERVs following viral infection, providing a valuable resource for future studies on ERVs and their interaction with SARS-CoV-2.

A chromosome-level genome assembly of the spider mite Tetranychus piercei McGregor

Despite the rapid advances in sequencing technology, limited genomic resources are currently available for phytophagous spider mites, which include many important agricultural pests. One of these pests is Tetranychus piercei (McGregor), a serious banana pest in East Asia exhibiting remarkable tolerance to high temperature. In this study, we assembled a high-quality genome of T. piercei using a combination of PacBio long reads and Illumina short reads sequencing. With the assistance of chromatin conformation capture technology, 99.9% of the contigs were anchored into three pseudochromosomes with a total size of 86.02 Mb. Repetitive elements, accounting for 14.16% of this genome (12.20 Mb), are predominantly composed of long-terminal repeats (30.7%). By combining evidence of ab initio prediction, transcripts, and homologous proteins, we annotated 11,881 protein-coding genes. Both the genome and proteins have high BUSCO completeness scores (>94%). This high-quality genome, along with reliable annotation, provides a valuable resource for investigating the high-temperature tolerance of this species and exploring the genomic basis that underlies the host range evolution of spider mites.

Base pair dynamics, electrostatics, and thermodynamics at the LTR-III quadruplex:duplex junction

G-quadruplexes (GQs) play key regulatory roles within the human genome and have also been identified to play similar roles in other eukaryotes, bacteria, archaea, and viruses. Human immunodeficiency virus 1, the etiological agent of acquired immunodeficiency syndrome, can form two GQs in its long terminal repeat (LTR) promoter region, each of which act to regulate viral gene expression in opposing manners. The major LTR GQ, called LTR-III, is a distinct hybrid GQ containing a 12-nucleotide duplex loop attached to the quadruplex motif. The resulting quadruplex:duplex junction (QDJ) has been hypothesized to serve as a selective drug targeting site. To better understand the dynamics of this QDJ, we performed conventional and enhanced-sampling molecular dynamics simulations using the Drude-2017 force field. We observed unbiased and reversible formation of additional base pairs in the QDJ, between Ade4:Thy14 and Gua3:Thy14. Both base pairs were electrostatically favored, but geometric constraints within the junction may drive the formation of, and preference for, the Ade4:Thy14 base pair. Finally, we demonstrated that the base pairs are separated only by small energy barriers that may enable transitions between both base-paired states. Together, these simulations provide new insights into the dynamics, electrostatics, and thermodynamics of the LTR-III QDJ.

Karyotype and LTR-RTs analysis provide insights into oak genomic evolution

BackgroundWhole-genome duplication and long terminal repeat retrotransposons (LTR-RTs) amplification in organisms are essential factors that affect speciation, local adaptation, and diversification of organisms. Understanding the karyotype projection and LTR-RTs amplification could contribute to untangling evolutionary history. This study compared the karyotype and LTR-RTs evolution in the genomes of eight oaks, a dominant lineage in Northern Hemisphere forests.ResultsKaryotype projections showed that chromosomal evolution was relatively conservative in oaks, especially on chromosomes 1 and 7. Modern oak chromosomes formed through multiple fusions, fissions, and rearrangements after an ancestral triplication event. Species-specific chromosomal rearrangements revealed fragments preserved through natural selection and adaptive evolution. A total of 441,449 full-length LTR-RTs were identified from eight oak genomes, and the number of LTR-RTs for oaks from section Cyclobalanopsis was larger than in other sections. Recent amplification of the species-specific LTR-RTs lineages resulted in significant variation in the abundance and composition of LTR-RTs among oaks. The LTR-RTs insertion suppresses gene expression, and the suppressed intensity in gene regions was larger than in promoter regions. Some centromere and rearrangement regions indicated high-density peaks of LTR/Copia and LTR/Gypsy. Different centromeric regional repeat units (32, 78, 79 bp) were detected on different Q. glauca chromosomes.ConclusionChromosome fusions and arm exchanges contribute to the formation of oak karyotypes. The composition and abundance of LTR-RTs are affected by its recent amplification. LTR-RTs random retrotransposition suppresses gene expression and is enriched in centromere and chromosomal rearrangement regions. This study provides novel insights into the evolutionary history of oak karyotypes and the organization, amplification, and function of LTR-RTs.

Isolation, characterisation and description of the roseoflavin producer Streptomyces berlinensis sp. nov.

The Gram-positive bacteria Streptomyces davaonensis and Streptomyces cinnabarinus have been the only organisms known to produce roseoflavin, a riboflavin (vitamin B2) derived red antibiotic. Using a selective growth medium and a phenotypic screening, we were able to isolate a novel roseoflavin producer from a German soil sample. The isolation procedure was repeated twice, that is, the same strain could be isolated from the same location in Berlin 6 months and 12 months after its first isolation. Whole genome sequencing of the novel roseoflavin producer revealed an unusual chromosomal arrangement and the deposited genome sequence of the new isolate (G + C content of 71.47%) contains 897 genes per inverted terminal repeat, 6190 genes in the core and 107 genes located on an illegitimate terminal end. We identified the roseoflavin biosynthetic genes rosA, rosB and rosC and an unusually high number of riboflavin biosynthetic genes. Overexpression of rosA, rosB and rosC in Escherichia coli and enzyme assays confirmed their predicted functions in roseoflavin biosynthesis. A full taxonomic analysis revealed that the isolate represents a previously unknown Streptomyces species and we propose the name Streptomyces berlinensis sp. nov. for this roseoflavin producer.

The near‐complete genome assembly of Reynoutria multiflora reveals the genetic basis of stilbene and anthraquinone biosynthesis

AbstractReynoutria multiflora is a widely used medicinal plant in China. Its medicinal compounds are mainly stilbenes and anthraquinones which possess important pharmacological activities in anti‐aging, anti‐inflammatory and anti‐oxidation, but their biosynthetic pathways are still largely unresolved. Here, we reported a near‐complete genome assembly of R. multiflora consisting of 1.39 Gb with a contig N50 of 122.91 Mb and only one gap left. Genome evolution analysis revealed that two recent bursts of long terminal repeats (LTRs) contributed significantly to the increased genome size of R. multiflora, and numerous large chromosome rearrangements were observed between R. multiflora and Fagopyrum tataricum genomes. Comparative genomics analysis revealed that a recent whole‐genome duplication specific to Polygonaceae led to a significant expansion of gene families associated with disease tolerance and the biosynthesis of stilbenes and anthraquinones in R. multiflora. Combining transcriptomic and metabolomic analyses, we elucidated the molecular mechanisms underlying the dynamic changes in content of medicinal ingredients in R. multiflora roots across different growth years. Additionally, we identified several putative key genes responsible for anthraquinone and stilbene biosynthesis. We identified a stilbene synthase gene PM0G05131 highly expressed in roost, which may exhibit an important role in the accumulation of stilbenes in R. multiflora. These genomic data will expedite the discovery of anthraquinone and stilbenes biosynthesis pathways in medicinal plants.