Terminal Repeat Research Articles

HHLA2 deficiency inhibits pancreatic cancer progression and THP-1 macrophage M2 polarization via EGFR/MAPK/ERK and mTOR/AKT pathway

BackgroundHuman endogenous retrovirus subfamily H long terminal repeat associating protein 2, (HHLA2), a member of B7 family, exhibits heightened expression in various malignant tumors. However, the exact functions of HHLA2 in pancreatic cancer (PC) remain incompletely elucidated.MethodsWe initially conducted an analysis of the B7 family members’ expression pattern in pancreatic tumor samples and adjacent normal tissues using The Cancer Genome Atlas (TCGA) database. Subsequently, immunohistochemistry, RT-qPCR and western blot methods were used to assess HHLA2 expression levels in PC tissues and cell lines. Furthermore, after silencing HHLA2 in PC cell lines, cell migration and proliferation of PC cells were detected by wound healing and CCK-8 assays, and cell invasion of PC cells was detected by transwell assays. We also investigated the regulation of epithelial-mesenchymal transition (EMT) markers and levels of EGFR, MEK, ERK1/2, mTOR and AKT via western blot analysis. Finally, the correlation between HHLA2 expression and immune infiltration was further explored.ResultsSilencing of HHLA2 resulted in the inhibition of PC cell proliferation, migration and invasion, potentially through the suppression of the EGFR/MAPK/ERK and mTOR/AKT signaling pathway. Additionally, silencing HHLA2 led to the inhibition of M2-type polarization of tumor associated macrophages (TAMs).ConclusionThe knockdown of HHLA2 was observed to inhibit the migration and invasion of PC cells through the regulation of the EMT process and EGFR/MAPK/ERK and mTOR/AKT pathway. Furthermore, silencing HHLA2 was found to modulate M2 polarization of TAMs. These finding suggest that HHLA2 could be a promising therapeutic target for Pancreatic cancer.

Fine mapping and candidate gene analysis of CRA8.1.6, which confers clubroot resistance in turnip (Brassica rapa ssp. rapa).

Clubroot disease poses a significant threat to Brassica crops, necessitating ongoing updates on resistance gene sources. In F2 segregants of the clubroot-resistant inbred line BrT18-6-4-3 and susceptible DH line Y510, the genetic analysis identified a single dominant gene responsible for clubroot resistance. Through bulk segregant sequencing analysis and kompetitive allele-specific polymerase chain reaction assays, CRA8.1.6 was mapped within 110 kb (12,255-12,365 Mb) between markers L-CR11 and L-CR12 on chromosome A08. We identified B raA08g015220.3.5C as the candidate gene of CRA8.1.6. Upon comparison with the sequence of disease-resistant material BrT18-6-4-3, we found 249 single-nucleotide polymorphisms, seven insertions, six deletions, and a long terminal repeat (LTR) retrotransposon (5,310 bp) at 909 bp of the first intron. However, the LTR retrotransposon was absent in the coding sequence of the susceptible DH line Y510. Given the presence of a non-functional LTR insertion in other materials, it showed that the LTR insertion might not be associated with susceptibility. Sequence alignment analysis revealed that the fourth exon of the susceptible line harbored two deletions and an insertion, resulting in a frameshift mutation at 8,551 bp, leading to translation termination at the leucine-rich repeat domain's C-terminal in susceptible material. Sequence alignment of the CDS revealed a 99.4% similarity to Crr1a, which indicate that CRA8.1.6 is likely an allele of the Crr1a gene. Two functional markers, CRA08-InDel and CRA08-KASP1, have been developed for marker-assisted selection in CR turnip cultivars. Our findings could facilitate the development of clubroot-resistance turnip cultivars through marker-assisted selection.

Genome-wide identification analysis in wild-type Solanum pinnatisectum reveals some genes defending against Phytophthora infestans.

Solanum pinnatisectum exhibits strong resistance to late blight caused by Phytophthora infestans but only an incomplete genome assembly based on short Illumina reads has been published. In this study, we generated the first chromosome-level draft genome for the wild-type potato species S. pinnatisectum in China using Oxford Nanopore technology sequencing and Hi-C technology. The high-quality assembled genome size is 664Mb with a scaffold N50 value of 49.17 Mb, of which 65.87% was occupied by repetitive sequences, and predominant long terminal repeats (42.51% of the entire genome). The genome of S. pinnatisectum was predicted to contain 34,245 genes, of which 99.34% were functionally annotated. Moreover, 303 NBS-coding disease resistance (R) genes were predicted in the S. pinnatisectum genome to investigate the potential mechanisms of resistance to late blight disease. The high-quality chromosome-level reference genome of S. pinnatisectum is expected to provide potential valuable resources for intensively and effectively investigating molecular breeding and genetic research in the future.

Anaerobic fungi contain abundant, diverse, and transcriptionally active Long Terminal Repeat retrotransposons

Long Terminal Repeat (LTR) retrotransposons are a class of repetitive elements that are widespread in the genomes of plants and many fungi. LTR retrotransposons have been associated with rapidly evolving gene clusters in plants and virulence factor transfer in fungal-plant parasite-host interactions. We report here the abundance and transcriptional activity of LTR retrotransposons across several species of the early-branching Neocallimastigomycota, otherwise known as the anaerobic gut fungi (AGF). The ubiquity of LTR retrotransposons in these genomes suggests key evolutionary roles in these rumen-dwelling biomass degraders, whose genomes also contain many enzymes that are horizontally transferred from other rumen-dwelling prokaryotes. Up to 10% of anaerobic fungal genomes consist of LTR retrotransposons, and the mapping of sequences from LTR retrotransposons to transcriptomes shows that the majority of clusters are transcribed, with some exhibiting expression greater than 104 reads per kilobase million mapped reads (rpkm). Many LTR retrotransposons are strongly differentially expressed upon heat stress during fungal cultivation, with several exhibiting a nearly three-log10 fold increase in expression, whereas growth substrate variation modulated transcription to a lesser extent. We show that some LTR retrotransposons contain carbohydrate-active enzymes (CAZymes), and the expansion of CAZymes within genomes and among anaerobic fungal species may be linked to retrotransposon activity. We further discuss how these widespread sequences may be a source of promoters and other parts towards the bioengineering of anaerobic fungi.

High Diversity of Long Terminal Repeat Retrotransposons in Compact Vertebrate Genomes: Insights from Genomes of Tetraodontiformes.

This study aimed to investigate the evolutionary profile (including diversity, activity, and abundance) of retrotransposons (RTNs) with long terminal repeats (LTRs) in ten species of Tetraodontiformes. These species, Arothron firmamentum, Lagocephalus sceleratus, Pao palembangensis, Takifugu bimaculatus, Takifugu flavidus, Takifugu ocellatus, Takifugu rubripes, Tetraodon nigroviridis, Mola mola, and Thamnaconus septentrionalis, are known for having the smallest genomes among vertebrates. Data mining revealed a high diversity and wide distribution of LTR retrotransposons (LTR-RTNs) in these compact vertebrate genomes, with varying abundances among species. A total of 819 full-length LTR-RTN sequences were identified across these genomes, categorized into nine families belonging to four different superfamilies: ERV (Orthoretrovirinae and Epsilon retrovirus), Copia, BEL-PAO, and Gypsy (Gmr, Mag, V-clade, CsRN1, and Barthez). The Gypsy superfamily exhibited the highest diversity. LTR family distribution varied among species, with Takifugu bimaculatus, Takifugu flavidus, Takifugu ocellatus, and Takifugu rubripes having the highest richness of LTR families and sequences. Additionally, evidence of recent invasions was observed in specific tetraodontiform genomes, suggesting potential transposition activity. This study provides insights into the evolution of LTR retrotransposons in Tetraodontiformes, enhancing our understanding of their impact on the structure and evolution of host genomes.

The first high-quality genome assembly and annotation of Lantana camara, an important ornamental plant and a major invasive species

This study presents the first annotated, haplotype-resolved, chromosome-scale genome of Lantana camara, a flowering shrub native to Central America and known for its dual role as an ornamental plant and an invasive species. Despite its widespread cultivation and ecological impact, the lack of a high-quality genome has hindered the investigation of traits of both ornamental and invasive. This research bridges the gap in genomic resources for L. camara, which is crucial for both ornamental breeding programs and invasive species management. Whole-genome and transcriptome sequencing were utilized to elucidate the genetic complexity of a diploid L. camara breeding line UF-T48. The genome was assembled de novo using HiFi and Hi-C reads, resulting in two phased genome assemblies with high Benchmarking Universal Single-Copy Orthologs (BUSCO) scores of 97.7%, indicating their quality. All 22 chromosomes were assembled with pseudochromosomes averaging 117 Mb. The assemblies revealed 29 telomeres and an extensive presence of repetitive sequences, primarily long terminal repeat transposable elements. The genome annotation identified 83,775 protein-coding genes, with 83% functionally annotated. In particular, the study mapped 42 anthocyanin and carotenoid candidate gene clusters and 12 herbicide target genes to the assembly, identifying 38 genes spread across the genome that are integral to flower color development and 53 genes for herbicide targeting in L. camara. This comprehensive genomic study not only enhances the understanding of L. camara’s genetic makeup but also sets a precedent for genomic research in the Verbenaceae family, offering a foundation for future studies in plant genetics, conservation, and breeding.

NF-κB1 deficiency promotes macrophage-derived adrenal tumors but decreases neurofibromas in HTLV-I LTR-Tax transgenic mice.

Human T-cell leukemia virus type I (HTLV-I) is an oncogenic virus whose infection can cause diverse diseases, most notably adult T-cell leukemia/lymphoma (ATL or ATLL), an aggressive and fatal malignancy of CD4 T cells. The oncogenic ability of HTLV-I is mostly attributed to the viral transcriptional transactivator Tax. Tax alone is sufficient to induce specific tumors in mice depending on the promotor used to drive Tax expression, thereby being used to understand HTLV-I tumorigenesis and model the tumor types developed in Tax transgenic mice. Tax exerts its oncogenic role predominantly by activating the cellular transcription factor NF-κB. Here, we report that genetic deletion of NF-κB1, the prototypic member of the NF-κB family, promotes adrenal medullary tumors but suppresses neurofibromas in mice with transgenic Tax driven by the HTLV-I Long Terminal Repeat (LTR) promoter. The adrenal tumors are derived from macrophages. Neoplastic macrophages also infiltrate the spleen and lymph nodes, causing splenomegaly and lymphadenopathy in mice. Nevertheless, the findings could be human relevant, because macrophages are important target cells of HTLV-I infection and serve as a virus reservoir in vivo. Moreover, the spleen, lymph nodes and adrenal glands are the most common sites of tumor cell infiltration in HTLV-I-infected patients. These data provide new mechanistic insights into the complex interaction between Tax and NF-κB, therefore improving our understanding of HTLV-I oncogenic pathogenesis. They also expand our knowledge and establish a new animal model of macrophage neoplasms and adrenal tumors.

Characterization of the CsCENH3 protein and centromeric DNA profiles reveal the structures of centromeres in cucumber.

Centromeres in eukaryotes mediate the accurate segregation of chromosomes during cell division. They serve as essential functional units of chromosomes and play a core role in the process of genome evolution. Centromeres are composed of satellite repeats and highly repetitive centromeric retrotransposons (CRs), which vary greatly even among closely related species. Cucumber (Cucumis sativus) is a globally cultivated and economically important vegetable and the only species in the Cucumis genus with seven pairs of chromosomes. Therefore, studying the centromeres of the Cucumis subgenus may yield valuable insights into its genome structure and evolution. Using chromatin immunoprecipitation (ChIP) techniques, we isolated centromeric DNA from cucumber reference line 9930. Our investigation into cucumber centromeres uncovered the centromeric satellite sequence, designated as CentCs, and the prevalence of Ty1/Copia long terminal repeat retrotransposons. In addition, active genes were identified in the CsCENH3 nucleosome regions with low transcription levels. To the best of our knowledge, this is the first time that characterization of centromeres has been achieved in cucumber. Meanwhile, our results on the distribution of CentCs and CsCRs in the subgenus Cucumis indicate that the content of centromeric repeats in the wild variants was significantly reduced compared with the cultivated cucumber. The results provide evidence for centromeric DNA amplification that occurred during the domestication process from wild to cultivated cucumber. Furthermore, these findings may offer new information for enhancing our understanding of phylogenetic relationships in the Cucumis genus.

Chromosome-scale genome assembly clarifies the mechanism of flooding tolerance and evolutionary history of Myricaria laxiflora

Myricaria laxiflora, belonging to the Tamaricaceae family, is a valuable plant resource naturally occurring along the riverbanks of the Yangtze River Basin. In this study, we constructed a chromosome-level M. laxiflora genome which was 1.34 Gb in size, with 12 chromosomes and a contig N50 of 11.69 Mb. In total, 23,719 protein-coding genes were discovered across the M. laxiflora genome. Phylogenetic investigation indicated that Tamaricaceae and Polygonaceae diverged between 65 and 81 million years ago (Mya). M. laxiflora underwent one whole genome duplication (WGD) between 19.74 and 24.60 Mya, and experienced a long terminal repeat (LTR) insertion around 0.15 Mya. These genomic changes may be related to environmental upheavals (such as the uplift of the Qinghai-Tibet Plateau) experienced by M. laxiflora during its evolution, and also may have provided a foundation for its evolution and adaptation. To analyze the flooding adaption mechanism of M. laxiflora, we compared morphological, physiological, and transcriptomic changes under control and submerged conditions. Overall, M. laxiflora was found to adopt a "quiescence" strategy when completely submerged. Specifically, in response to submersion, photosynthesis weakens, growth stagnates, glycolysis is activated, and the antioxidant enzyme system is enhanced. The “ethylene response pathway” was ascertained as a likely regulator of this process, and a MlERF-VII gene (Mla11G026380) may be the key gene regulating this process. Our findings shed light on the underlying mechanisms behind the flooding resistance displayed by M. laxiflora, while also establishing a basis for the conservation of its genetic resources.

A unified classification system for HIV-1 5' long terminal repeats.

The HIV-1 provirus mainly consists of internal coding region flanked by 1 long terminal repeats (LTRs) at each terminus. The LTRs play important roles in HIV-1 reverse transcription, integration, and transcription. However, despite of the significant study advances of the internal coding regions of HIV-1 by using definite reference classification, there are no systematic and phylogenetic classifications for HIV-1 5' LTRs, which hinders our elaboration on 5' LTR and a better understanding of the viral origin, spread and therapy. Here, by analyzing all available resources of 5' LTR sequences in public databases following 4 recognized principles for the reference classification, 83 representatives and 14 consensus sequences were identified as representatives of 2 groups, 6 subtypes, 6 sub-subtypes, and 9 CRFs. To test the reliability of the supplemented classification system, the constructed references were applied to identify the 5' LTR assignment of the 22 clinical isolates in China. The results revealed that 16 out of 22 tested strains showed a consistent subtype classification with the previous LTR-independent classification system. However, 6 strains, for which recombination events within 5' LTR were demonstrated, unexpectedly showed a different subtype classification, leading a significant change of binding sites for important transcription factors including SP1, p53, and NF-κB. The binding change of these transcriptional factors would probably affect the transcriptional activity of 5' LTR. This study supplemented a unified classification system for HIV-1 5' LTRs, which will facilitate HIV-1 characterization and be helpful for both basic and clinical research fields.

KSHV genome harbors both constitutive and lytically induced enhancers.

Kaposi's sarcoma-associated herpesvirus (KSHV) belongs to the gamma-herpesvirus family and is a well-known human oncogenic virus. In infected cells, the viral genome of 165 kbp is circular DNA wrapped in chromatin. The tight control of gene expression is critical for latency, the transition into the lytic phase, and the development of viral-associated malignancies. Distal cis-regulatory elements, such as enhancers and silencers, can regulate gene expression in a position- and orientation-independent manner. Open chromatin is another characteristic feature of enhancers. To systematically search for enhancers, we cloned all the open chromatin regions in the KSHV genome downstream of the luciferase gene and tested their enhancer activity in infected and uninfected cells. A silencer was detected upstream of the latency-associated nuclear antigen promoter. Two constitutive enhancers were identified in the K12p-OriLyt-R and ORF29 Intron regions, where ORF29 Intron is a tissue-specific enhancer. The following promoters: OriLyt-L, PANp, ALTp, and the terminal repeats (TRs) acted as lytically induced enhancers. The expression of the replication and transcription activator (RTA), the master regulator of the lytic cycle, was sufficient to induce the activity of lytic enhancers in uninfected cells. We propose that the TRs that span about 24 kbp region serve as a "viral super-enhancer" that integrates the repressive effect of the latency-associated nuclear antigen (LANA) with the activating effect of RTA. Utilizing CRISPR activation and interference techniques, we determined the connections between these enhancers and their regulated genes. The silencer and enhancers described here provide an additional layer to the complex gene regulation of herpesviruses.IMPORTANCEIn this study, we performed a systematic functional assay to identify cis-regulatory elements within the genome of the oncogenic herpesvirus, Kaposi's sarcoma-associated herpesvirus (KSHV). Similar to other herpesviruses, KSHV presents both latent and lytic phases. Therefore, our assays were performed in uninfected cells, during latent infection, and under lytic conditions. We identified two constitutive enhancers, one of which seems to be a tissue-specific enhancer. In addition, four lytically induced enhancers, which are all responsive to the replication and transcription activator (RTA), were identified. Furthermore, a silencer was identified between the major latency promoter and the lytic gene locus. Utilizing CRISPR activation and interference techniques, we determined the connections between these enhancers and their regulated genes. The terminal repeats, spanning a region of about 24 kbp, seem like a "viral super-enhancer" that integrates the repressive effect of the latency-associated nuclear antigen (LANA) with the activating effect of RTA to regulate latency to lytic transition.

Chromosome-level genome assembly of Pedicularis kansuensis illuminates genome evolution of facultative parasitic plant.

Parasitic plants have a heterotrophic lifestyle, in which they withdraw all or part of their nutrients from their host through the haustorium. Despite the release of many draft genomes of parasitic plants, the genome evolution related to the parasitism feature of facultative parasites remains largely unknown. In this study, we present a high-quality chromosomal-level genome assembly for the facultative parasite Pedicularis kansuensis (Orobanchaceae), which invades both legume and grass host species in degraded grasslands on the Qinghai-Tibet Plateau. This species has the largest genome size compared with other parasitic species, and expansions of long terminal repeat retrotransposons accounting for 62.37% of the assembly greatly contributed to the genome size expansion of this species. A total of 42,782 genes were annotated, and the patterns of gene loss in P. kansuensis differed from other parasitic species. We also found many mobile mRNAs between P. kansuensis and one of its host species, but these mobile mRNAs could not compensate for the functional losses of missing genes in P. kansuensis. In addition, we identified nine horizontal gene transfer (HGT) events from rosids and monocots, as well as one single-gene duplication events from HGT genes, which differ distinctly from that of other parasitic species. Furthermore, we found evidence for HGT through transferring genomic fragments from phylogenetically remote host species. Taken together, these findings provide genomic insights into the evolution of facultative parasites and broaden our understanding of the diversified genome evolution in parasitic plants and the molecular mechanisms of plant parasitism.

Identification of transcription factors contributing to vitamin C synthesis during Rosa roxburghii fruit development by integrating transcriptomics and metabolomics

Ascorbic acid, also referred to as vitamin C (Vc), is an important nutrient found in fruits and vegetables that promotes produce quality and human health. Rosa roxburghii is an underutilized natural fruit that contains very high levels of Vc. However, the Vc content of R. roxburghii varies considerably during plant development and ripening. To better understand the molecular mechanisms that underlie fluctuations in Vc content of R. roxburghii fruit at different developmental stages, we performed transcriptomic and metabolomic analyses and identified two significant gene networks/modules and 168 transcription factors directly involved in Vc synthesis. Promoter analysis of two core genes involved in Vc synthesis, RrGGP and RrGalUR, revealed the presence of a retroviral long terminal repeat (LTR) insert in the RrGalUR promoter. Using yeast one-hybrid and dual-luciferase assays, we demonstrated that the transcription factors RrHY5H and RrZIP9 bind to the promoter of RrGGP to promote its expression. RrZIP6 and RrWRKY4 bind to the LTR in the RrGalUR promoter to promote its expression. Our results reveal a molecular mechanism that controls Vc synthesis and accumulation in R. roxburghii fruit.

Elucidating the multichromosomal structure within the Brasenia schreberi mitochondrial genome through assembly and analysis

Brasenia schreberi, a plant species traditionally utilized in Chinese medicine and cuisine, represents an early evolutionary stage among flowering plants (angiosperms). While the plastid genome of this species has been published, its mitochondrial genome (mitogenome) has not been extensively explored, with a notable absence of thorough comparative analyses of its organellar genomes. In our study, we had assembled the entire mitogenome of B. schreberi utilizing the sequencing data derived from both Illumina platform and Oxford Nanopore. The B. schreberi mitogenome mostly exists as six circular DNA molecules, with the largest being 628,257 base pairs (bp) and the smallest 110,220 bp, amounting to 1.49 megabases (Mb). Then we annotated the mitogenome of B. schreberi. The mitogenome encompasses a total of 71 genes: 40 of these are coding proteins genes (PCGs), 28 are genes for transfer RNA (tRNA), and the remaining 3 are genes for ribosomal RNA (rRNA). In the analysis of codon usage, we noted a unique codon preference specific to each amino acid. The most commonly used codons exhibited an average RSCU of 1.36, indicating a noticeable bias in codon selection. In the repeat sequence analysis, a total of 553 simple sequence repeats (SSRs) were identified, 1,822 dispersed repeats (comprising 1,015 forward and 807 palindromic repeats), and 608 long terminal repeats (LTRs). Additionally, in the analysis of homologous sequences between organelle genomes, we detected 38 homologous sequences derived from the plastid genome, each exceeding 500 bp, within the B. schreberi mitochondrial genome. Notably, ten tRNA genes (trnC-GCA, trnM-CAU, trnI-CAU, trnQ-UUG, trnN-GUU, trnT-GGU, trnW-CCA, trnA-UGC, trnI-GAU, and trnV-GAC) appear to have been completely transferred from the chloroplast to the mitogenome. Utilizing the Deepred-mt to predict the RNA editing sites in the mitogenome, we have identified 675 high-quality RNA editing sites in the 40 mitochondrial PCGs. In the final stage of our study, we performed an analysis of colinearity and inferred the phylogenetic relationship of B. schreberi with other angiosperms, utilizing the mitochondrial PCGs as a basis. The results showed that the non-coding regions of the B. schreberi mitogenome are characterized by an abundance of repetitive sequences and exogenous sequences, and B. schreberi is more closely related with Euryale ferox.

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.