Terminal Repeat Research Articles

Cocaine-Induced DNA-Dependent Protein Kinase Relieves RNAP II Pausing by Promoting TRIM28 Phosphorylation and RNAP II Hyperphosphorylation to Enhance HIV Transcription

Drug abuse continues to pose a significant challenge in HIV control efforts. In our investigation, we discovered that cocaine not only upregulates the expression of the DNA-dependent protein kinase (DNA-PK) but also augments DNA-PK activation by enhancing its phosphorylation at S2056. Moreover, DNA-PK phosphorylation triggers the higher localization of the DNA-PK into the nucleus. The finding that cocaine increases the nuclear localization of the DNA-PK provides further support to our observation of enhanced DNA-PK recruitment at the HIV long terminal repeat (LTR) following cocaine exposure. By activating and facilitating the nuclear localization of the DNA-PK, cocaine effectively orchestrates multiple stages of HIV transcription, thereby promoting HIV replication. Additionally, our study demonstrates that the cocaine-induced DNA-PK promotes the hyper-phosphorylation of the RNA polymerase II (RNAP II) carboxyl-terminal domain (CTD) at Ser5 and Ser2 sites, enhancing both the initiation and elongation phases, respectively, of HIV transcription. The cocaine-mediated enhancement of transcriptional initiation is supported by its activation of cyclin-dependent kinase 7 (CDK7). Additionally, the induction of transcriptional elongation is marked by higher LTR recruitment and the increased phosphorylation of CDK9, which indicates the stimulation of positive transcriptional elongation factor b (P-TEFb). We demonstrate for the first time that cocaine, through DNA-PK activation, promotes the specific phosphorylation of TRIM28 at serine 824 (p-TRIM28, S824). This modification converts TRIM28 from a transcriptional inhibitor to a transactivator for HIV transcription. Additionally, we observed that the phosphorylation of TRIM28 (p-TRIM28, S824) promotes the transition from the pausing phase to the elongation phase of HIV transcription, thereby facilitating the production of full-length HIV genomic transcripts. This finding corroborates the previously observed enhanced RNAP II CTD phosphorylation at Ser2, a marker of transcriptional elongation, following cocaine exposure. Accordingly, upon cocaine treatment, we observed the elevated recruitment of p-TRIM28-(S824) at the HIV LTR. Overall, our results unravel the intricate molecular mechanisms underlying cocaine-induced HIV transcription and gene expression. These findings hold promise for the development of highly targeted therapeutics aimed at mitigating the detrimental effects of cocaine in individuals living with HIV.

Genome sequences of five Klebsiella bacteriophages that belong to the genus Jiaodavirus.

We describe the genomes of five lytic Klebsiella pneumoniae myophages, therapeutic candidates, that belong to the family Straboviridae and genus Jiaodavirus. The genomes ranged from 165,574 to 169,768 bp, with ca. 40% GC content, contained 289-300 coding sequences, had 15-16 tRNA genes, and no terminal repeats.

Investigation of Polymorphisms Induced by the Solo Long Terminal Repeats (Solo-LTRs) in Porcine Endogenous Retroviruses (ERVs).

Homologous recombination events take place between the 5' and 3' long terminal repeats (LTRs) of ERVs, resulting in the generation of solo-LTR, which can cause solo-LTR-associated polymorphism across different genomes. In the current study, specific criteria were established for the filtration of solo-LTRs, resulting in an average of 5630 solo-LTRs being identified in 21 genomes. Subsequently, a protocol was developed for detecting solo-LTR polymorphisms in the pig genomes, resulting in the discovery of 927 predicted solo-LTR polymorphic sites. Following verification and filtration processes, 603 highly reliable solo-LTR polymorphic sites were retained, involving 446 solo-LTR presence sites (solo-LTR+) and 157 solo-LTR absence sites (solo-LTR-) relative to the reference genome. Intersection analysis with gene/functional regions revealed that 248 solo-LTR- sites and 23 solo-LTR+ sites overlapped with genes or were in the vicinity of genes or functional regions, impacting a diverse range of gene structures. Moreover, through the utilization of 156 solo-LTR polymorphic sites for population genetic analysis, it was observed that these solo-LTR loci effectively clustered various breeds together, aligning with expectations and underscoring their practical utility. This study successfully established a methodology for detecting solo-LTR polymorphic sites. By applying these methods, a total of 603 high-reliability solo-LTR polymorphic sites were pinpointed, with nearly half of them being linked to genes or functional regions.

Chromosome scale genome assembly and annotation of coconut cultivar Chowghat Green Dwarf

The high-quality genome of coconut (Cocos nucifera L.) is a crucial resource for enhancing agronomic traits and studying genome evolution within the Arecaceae family. We sequenced the Chowghat Green Dwarf cultivar, which is resistant to the root (wilt) disease, utilizing Illumina, PacBio, ONT, and Hi-C technologies to produce a chromosome-level genome of ~ 2.68 Gb with a scaffold N50 of 174 Mb; approximately 97% of the genome could be anchored to 16 pseudo-molecules (2.62 Gb). In total, 34,483 protein-coding genes were annotated; the BUSCO completeness score was 96.80%, while the k-mer completeness was ~ 87%. The assembled genome includes 2.19 Gb (81.64%) of repetitive sequences, with long terminal repeats (LTRs) constituting the most abundant class at 53.76%. Additionally, our analysis confirms two whole-genome duplication (WGD) events in the C. nucifera lineage. A genome-wide analysis of LTR insertion time revealed ancient divergence and proliferation of copia and gypsy elements. In addition, 1368 RGAs were discovered in the CGD genome. We also developed a web server ‘Kalpa Genome Resource’ (http://210.89.54.198:3000/), to manage and store a comprehensive array of genomic data, including genome sequences, genetic markers, structural and functional annotations like metabolic pathways, and transcriptomic profiles. The web server has an embedded genome browser to analyze and visualize the genome, its genomics elements, and transcriptome data. The in-built BLAST server allows sequence homology searches against genome, annotated transcriptome & proteome sequences. The genomic dataset and the database will support comparative genome analysis and can expedite genome-driven breeding and enhancement efforts for tapping genetic gains in coconut.

A Novel Tax-Responsive Reporter T-Cell Line to Analyze Infection of HTLV-1.

Human T-cell leukemia virus type 1 (HTLV-1) infects CD4+ T-cells through close cell-cell contacts. The viral Tax-1 (Tax) protein regulates transcription by transactivating the HTLV-1 U3R promoter in the 5' long terminal repeat of the integrated provirus. Here, we generated a clonal Tax-responsive T-cell line to track HTLV-1 infection at the single-cell level using flow cytometry, bypassing intracellular viral protein staining. Jurkat T-cells stably transduced with the SMPU vector carrying green fluorescent protein (GFP) under control of 18 × 21 bp Tax-responsive element repeats of the U3R were evaluated. Among 40 clones analyzed for Tax responsiveness, the top two were characterized. Upon overexpression of Tax, over 40% of the cells showed GFP positivity, and approximately 90% of the Tax-positive cells were GFP-positive, indicating efficient reporter activity. However, with CREB-deficient Tax mutant M47, both total GFP-positive cell counts and those within the Tax-positive group significantly decreased. Co-culture with chronically HTLV-1-infected MT-2 or C91-PL cells led to an average of 0.9% or 2.4% GFP-positive cells, respectively, confirming the suitability to monitor HTLV-1 transmission and that HTLV-1 infection is very low. Thus, the novel Tax-responsive reporter T-cell line is a suitable tool to monitor infection of HTLV-1 on the single-cell level.

The complete mitochondrial DNA of the carnivorous sponge Lycopodina hypogea is putatively complemented by microDNAs.

Carnivorous sponges (Porifera, Demospongiae, Cladorhizidae), contrary to the usual filter-feeding mechanism of sponges, are specialized in catching larger prey through adhesive surfaces or hook-like spicules. The mitochondrial DNA of sponges overall present several divergences from other metazoans, and while presenting unique features among major transitions, such as in calcarean and glass sponges, poriferan mitogenomes are relatively stable within their groups. Here, we report and discuss the mitogenome of Lycopodina hypogea (Vacelet & Boury-Esnault, 1996), which greatly vary from its subordinal counterparts in both structure and gene order. This mitogenome is seemingly multipartite into three chromosomes, two of them as microDNAs. The main chromosome, chrM1, is unusually large, 31,099 bp in length, has a unique gene order within Poecilosclerida, and presents two rRNA, 13 protein and 19 tRNA coding genes. Intergenic regions comprise approximately 40% of chrM1, bearing several terminal direct and inverted repeats (TDRr and TIRs) but holding no vestiges of former mitochondrial sequences, pseudogenes, or transposable elements. The nd4l and trnI(gau) genes are likely located in microDNAs thus comprising putative mitochondrial chromosomes chrM2, 291 bp, and chrM3, 140 bp, respectively. It is unclear which processes are responsible for the remarkable features of the of L. hypogea mitogenome, including a generalized gene rearrangement, long IGRs, and putative extrachromosomal genes in microDNAs.

Characterization and organization of telomeric-linked helicase (tlh) gene families in Fusarium oxysporum.

Telomere-linked RecQ helicase (tlh) genes have been reported in several fungi and a choanoflagellate in the regions adjacent to the terminal telomere repeats. In this study, we explored the Telomere-linked RecQ helicase (tlh) genes in four strains of Fusarium oxysporum, offering new insights into their genomic structure, functional motifs, and impact on chromosomal ends. We conducted a comprehensive analysis, comparing the tlh genes of F. oxysporum with those previously identified in other organisms and uncovering significant similarities. Through comparative genomics, we identified conserved protein motifs across these genes, including a TLH domain, C2H2, and RecQ helicase motifs. Our phylogenetic analysis positions the F. oxysporum tlh genes in a cluster with other known tlhs, suggesting a shared evolutionary origin. Mutation analysis revealed a relatively low level of deleterious mutations in tlh gene paralogs, with a notable proportion of full-size structures maintained across strains. Analysis of subtelomeric sequences indicates that a region with almost identical sequences flanks the majority of chromosome ends, termed tlh-containing region (TLHcr), across these strains. The presence of TLHcrs at chromosome ends, either as single entities or in arrays, underscores their potential role in telomere function and genome stability. Our findings provide a detailed examination of tlh genes in four strains of F. oxysporum, laying the groundwork for future studies on their biological significance and evolutionary history in fungal genomes.

Stabilization of a single-stranded DNA of adeno-associated virus by inverted terminal repeats

Parvoviruses have evolved to possess a linear single-stranded DNA (ssDNA) genome ranging from 4 to 6.3 kb. Adeno-associated virus (AAV), a member of the Parvoviridae family, contains approximately 5 kb of linear ssDNA within its capsid. This ssDNA features two 145-base inverted terminal repeats (ITRs) positioned at each end. ITRs have a T-shaped hairpin secondary structure, which plays a crucial role in viral replication. To investigate the impact of ITRs on ssDNA stability, we conducted a DNA denaturation-reannealing assay in 10 mM magnesium acetate, 50 mM potassium acetate, and 20 mM Tris-acetate buffer at pH7.9. Conventional double-stranded DNA (dsDNA) fragments retain a reannealing capability of over 50% for sizes under 8.8 kb, gradually losing this capability as sizes increase; however, dsDNA fragments in rAAV ranging from 0.7 to 6.3 kb did not exhibit a reannealing profile. This suggests that the presence of ITRs at both ends hinders annealing between the complementary strands. These results indicate that ITR structures preferentially induce an ssDNA conformation less than 6.3 kb in size, and that the stability of AAV ssDNA contributes to the viral life cycle, including processes such as infection, replication, and packaging. Considering the size of parvovirus genomes, it appears that their genomes require reversible flexibility in complementary DNA strands; simultaneously, this adaptability needs to be regulated by specific palindromic ITRs at both ends.

Natural Diversity of Heat-Induced Transcription of Retrotransposons in Arabidopsis thaliana.

Transposable elements (TEs) are major components of plant genomes, profoundly impacting the fitness of their hosts. However, technical bottlenecks have long hindered our mechanistic understanding of TEs. Using RNA-Seq and long-read sequencing with Oxford Nanopore Technologies' (ONT) direct cDNA sequencing, we analyzed the heat-induced transcription of TEs in three natural accessions of Arabidopsis thaliana (Cvi-0, Col-0, and Ler-1). In addition to the well-studied ONSEN retrotransposon family, we confirmed Copia-35 as a second heat-responsive retrotransposon family with particularly high activity in the relict accession Cvi-0. Our analysis revealed distinct expression patterns of individual TE copies and suggest different mechanisms regulating the GAG protein production in the ONSEN versus Copia-35 families. In addition, analogously to ONSEN, Copia-35 activation led to the upregulation of flanking genes such as APUM9 and potentially to the quantitative modulation of flowering time. ONT data allowed us to test the extent to which read-through formation is important in the regulation of adjacent genes. Unexpectedly, our results indicate that for both families, the upregulation of flanking genes is not predominantly directly initiated by transcription from their 3' long terminal repeats. These findings highlight the intraspecific expressional diversity linked to retrotransposon activation under stress.

Reduction of antisense transcription affects bovine leukemia virus replication and oncogenesis.

In sheep infected with bovine leukemia virus (BLV), transcription of structural, enzymatic, and accessory genes is silenced. However, the BLV provirus transcribes a series of non-coding RNAs that remain undetected by the host immune response. Specifically, three RNAs (AS1-L, AS1-S, and AS2) are consistently expressed from the antisense strand, originating from transcriptional initiation at the 3'-Long Terminal Repeat (LTR). To investigate the role of these non-coding RNAs in viral replication and pathogenesis, a reverse genetics approach was devised, capitalizing on a mechanistic disparity in transcription initiation between the 5' and 3' promoters. A two-nucleotide mutation (GG>TA) in the TFIIB-recognition element (BRE) impaired antisense transcription originating from the 3'-LTR. In the context of the provirus, this 2bp mutation significantly diminished the expression of antisense RNAs, while not notably affecting sense transcription. When inoculated to sheep, the mutated provirus was infectious but exhibited reduced replication levels, shedding light on the role of antisense transcription in vivo. In comparison to lymphoid organs in sheep infected with a wild-type (WT) provirus, the mutant demonstrated alterations in both the spatial distribution and rates of cell proliferation in the lymph nodes and the spleen. Analysis through RNA sequencing and RT-qPCR unveiled an upregulation of the Hmcn1/hemicentin-1 gene in B-lymphocytes from sheep infected with the mutated provirus. Further examination via confocal microscopy and immunohistochemistry revealed an increase in the amount of hemicentin-1 protein encoded by Hmcn1 in peripheral blood mononuclear cells (PBMCs) and lymphoid organs of sheep infected with the mutant. RNA interference targeting Hmcn1 expression impacted the migration of ovine kidney (OVK) cells in vitro. In contrast to the WT, the mutated provirus showed reduced oncogenicity when inoculated into sheep. Collectively, this study underscores the essential role of antisense transcription in BLV replication and pathogenicity. These findings may offer valuable insights into understanding the relevance of antisense transcription in the context of human T-cell leukemia virus (HTLV-1).

Discovery of the first Tn630 member and the closest homolog of IS630 from viruses

IS630/Tc1/mariner (ITm) represents the most widely distributed superfamily of DNA transposons in nature. Currently, bioinformatics research on ITm members primarily involves collecting data of existing and emerging members and organizing them into new groups or families. In the present study, our survey revealed that Tc1 and IS630 members have a broad host range, spanning across all six biological kingdoms (bacteria, fungi, plantae, animalia, archaea and protista) and viruses. The primary discoveries include the first Tn630 member—Tn630-NC1 and the closest homolog of IS630 from viruses—Tc1-C#1. By incorporating our discoveries into existing knowledge, we proposed a model to elucidate the formation of composite transposons. Organization of Tc1 and IS630 members into groups across biological kingdoms facilitates data collection for future research, particularly on their horizontal transfer between different kingdoms. The formation of composite transposons may result from asymmetric of terminal inverted repeats. IS630 should be merged with Tc1 into a single family IS630/Tc1. Furthermore, IS630 and its homologs constitute a valuable resource for studying horizontal gene transfer between gut bacteria and phages, opening up new avenues for research in this field.

The first chicken oocyte nucleus whole transcriptomic profile defines the spectrum of maternal mRNA and non-coding RNA genes transcribed by the lampbrush chromosomes.

Lampbrush chromosomes, with their unusually high rate of nascent RNA synthesis, provide a valuable model for studying mechanisms of global transcriptome up-regulation. Here, we obtained a whole-genomic profile of transcription along the entire length of all lampbrush chromosomes in the chicken karyotype. With nuclear RNA-seq, we obtained information about a wider set of transcripts, including long non-coding RNAs retained in the nucleus and stable intronic sequence RNAs. For a number of protein-coding genes, we visualized their nascent transcripts on the lateral loops of lampbrush chromosomes by RNA-FISH. The set of genes transcribed on the lampbrush chromosomes is required for basic cellular processes and is characterized by a broad expression pattern. We also present the first high-throughput transcriptome characterization of miRNAs and piRNAs in chicken oocytes at the lampbrush chromosome stage. Major targets of predicted piRNAs include CR1 and long terminal repeat (LTR) containing retrotransposable elements. Transcription of tandem repeat arrays was demonstrated by alignment against the whole telomere-to-telomere chromosome assemblies. We show that transcription of telomere-derived RNAs is initiated at adjacent LTR elements. We conclude that hypertranscription on the lateral loops of giant lampbrush chromosomes is required for synthesizing large amounts of transferred to the embryo maternal RNA for thousands of genes.

Chromosome-scale genome assembly of Codonopsis pilosula and comparative genomic analyses shed light on its genome evolution.

Codonopsis pilosula is a significant plant in traditional Chinese medicine, valued for its edible and medicinal properties. However, the lack of available genomic resources has hindered further research. This study presents the first chromosome-scale genome assembly of C. pilosula using PacBio CLR reads and Hi-C scaffolding technology. Additionally, Ks analysis and syntenic depth analysis were performed to elucidate its evolutionary history. The final assembly yielded a high-quality genome of 679.20 Mb, which was anchored to 8 pseudo-chromosomes with an anchoring rate of 96.5% and a scaffold N50 of 80.50 Mb. The genome assembly showed a high completeness of 97.6% based on Benchmarking with Universal Single-Copy Orthologs (BUSCO) analysis. Repetitive elements constituted approximately 76.8% of the genome, with long terminal repeat retrotransposons (LTRs) accounting for about 39.17%. Ks and syntenic depth analyses revealed that the polyploidization history of three platycodonoid clade species involved only the γ-WGT event. Karyotype evolutionary analysis identified an ancestral karyotype with 9 protochromosomes for the three platycodonoid clade species. Moreover, non-WGD genes, particularly those arising from tandem duplications, were found to contribute significantly to gene family expansion. These findings provide essential insights into the genetic diversity and evolutionary biology of C. pilosula, aiding its conservation and sustainable use.

Clinical Significance of HHLA2 as a Novel Therapeutic Target for Colorectal Cancer.

Colorectal cancer (CRC) is a high-indence malignance of the digestive system with a high mortality rate in the world. The results are desired to provide an important theoretical basis for discovering new therapeutic targets for CRC. The expression of human endogenous retrovirus-H-long terminal repeat association protein 2 (HHLA2) in human CRC was detected to explore its correlationship with clinicopathological features and prognosis of patients and its potential in treating CRC. Western blot was employed to detect HHLA2 expression in fresh tissues obtained from 6 CRC patients' excisions, including cancer, paracancer, and normal issues. Immunohistochemical staining was employed to determine HHLA2 expression in paraffin-embedded specimens of 139 patients with colorectal cancer, and its relationship with the clinicopathological profiles and survival was analyzed. Small interfering RNA (siRNA) targeting HHLA2 was used to transfect CRC cells to silent HHLA2. MTT, plate colony formation, cell scratch, and Transwell assay were conducted to observe the proliferation, migration, and invasion of CRC cells. HHLA2 protein was expressed in human colorectal cancer tissues, paracancer tissues and normal tissues, which was significantly upregulated in cancer tissues (P<0.01). HHLA2 expression level in CRC tissues showed a close correlationship with the invasion depth of the tumor (P=0.000), metastasis of regional lymph nodes (P=0.018), clinical stage (P=0.010), and patient survival (P=0.011). Correlation with gender (P=0.873), age (P=0.864), location of the tumor (P=0.768), degree of tumor differentiation (P=0.569) and distant metastasis (P=0.494) exhibited no significance. The survival time of CRC patients with high and low HHLA2 expression groups was 43.231 months and 55.649 months, respectively, with a statistical difference between the two groups (P=0.001). Silencing HHLA2 inhibited proliferation, migration and invasion of CRC cells significantly. HHLA2 is overexpressed in CRC tissues which is associated with poor prognosis of patients. HHLA2 might be recognized as a new candidate for adjuvant diagnosis and prognosis of CRC, as well as a promised new target for immunotherapy of CRC.

First Molecular Characterization of Small Ruminant Lentiviruses in Hungarian Goat Population.

In 2023, a molecular study was conducted on the Hungarian goat population to determine genotypes and subtypes of small ruminant lentiviruses (SRLV) infecting these herds. Ten goat herds seropositive for SRLV infection according to a serosurvey conducted earlier in Hungary were selected, and 135 adult goats (>1 year old) were blood sampled. The two-stage nested real-time PCR (nRT-PCR) was used to detect proviral DNA of SRLV and distinguish between two main viral genotypes (A and B). PCR products were submitted for Sanger dideoxy sequencing, and phylogenetic and molecular evolutionary analyses were conducted on the 200-250 bp-long proviral DNA sequences from the end of long terminal repeat (LTR) region and beginning of gag gene using the MEGA11 software. Reference strains included strains most identical to Hungarian sequences according to the Standard Nucleotide BLAST and prototypic strains for the relevant genotypes and subtypes. Proviral DNA of SRLV was detected in goats from all ten tested herds. A single SRLV genotype was detected in 6 herds-genotype A in three herds and B also in three herds. In four herds, mixed infection with genotypes A and B was confirmed. In total, 110/135 seropositive goats tested positive in the nRT-PCR (81.5%): 49/110 goats (44.5%) for genotype A, 54/110 goats (49.1%) for genotype B, and 7/110 goats (6.4%) for both genotypes. Hungarian sequences belonged to subtypes A1/A18, A2, and subtype B1. This is the first study which shows that Hungarian goats are infected by SRLV belonging to both genotypes A and B.

Inversions encounter relaxed genetic constraints and balance birth and death of TPS genes in Curcuma

Evolutionary dynamics of inversion and its impact on biochemical traits are a puzzling question. Here, we show abundance of inversions in three Curcuma species (turmeric, hidden ginger and Siam tulip). Genes within inversions display higher long terminal repeat content and lower expression level compared with genomic background, suggesting inversions in Curcuma experience relaxed genetic constraints. It is corroborated by depletion of selected SNPs and enrichment of deleterious mutations in inversions detected among 56 Siam tulip cultivars. Functional verification of tandem duplicated terpene synthase (TPS) genes reveals that genes within inversions become pseudogenes, while genes outside retain catalytic function. Our findings suggest that inversions act as a counteracting force against tandem duplication in balancing birth and death of TPS genes and modulating terpenoid contents in Curcuma. This study provides an empirical example that inversions are likely not adaptive but affect biochemical traits.

Clam Genome and Transcriptomes Provide Insights into Molecular Basis of Morphological Novelties and Adaptations in Mollusks.

Bivalve mollusks, comprising animals enclosed in two shell valves, are well-adapted to benthic life in many intertidal zones. Clams have evolved the buried lifestyle, which depends on their unique soft tissue structure and their wedge-shaped muscular foot and long extendible siphons. However, molecular mechanisms of adaptative phenotype evolution remain largely unknown. In the present study, we obtain the high-quality chromosome-level genome of Manila clam R. philippinarum, an economically important marine bivalve in many coastal areas. The genome is constructed by the Hi-C assisted assembly, which yields 19 chromosomes with a total of 1.17 Gb and BUSCO integrity of 92.23%. The de novo assembled genome has a contig N50 length of 307.7 kb and scaffold N50 of 59.5 Mb. Gene family expansion analysis reveals that a total of 24 single-copy gene families have undergone the significant expansion or contraction, including E3 ubiquitin ligase and dynein heavy chain. The significant expansion of transposable elements has been also identified, including long terminal repeats (LTR) and non-LTR retrotransposons. The comparative transcriptomics among different clam tissues reveals that extracellular matrix (ECM) receptors and neuroactive ligand receptors may play the important roles in tissue structural support and neurotransmission during their infaunal life. These findings of gene family expansion and tissue-specific expression may reflect the unique soft tissue structure of clams, suggesting the evolution of lineage-specific morphological novelties. The high-quality genome and transcriptome data of R. philippinarum will not only facilitate the genetic studies on clams but will also provide valuable information on morphological novelties in mollusks.

Multiple functions of the herpesvirus UL14 gene product in viral infection.

Herpesviruses are a family of double-stranded DNA viruses with a tegument structure and a genome composed of a single sequence and terminal repeat (TR) sequences. The herpesvirus UL14 gene encodes the protein UL14 (pUL14), which has various subcellular localizations and plays a vital role in regulating immediate-early (IE) gene transcription and expression, influences the intracellular localization patterns of several proteins belonging to the capsid and the DNA packaging machinery, participates in secondary envelopment, and influences viral particle release. Additionally, pUL14 has roles in maintaining cellular homeostasis and preventing apoptosis. This review discusses how pUL14 engages in the life cycle of herpesviruses and provides new ideas for further research on pUL14's function in viral infection.

Haplotype-resolved genome assembly and resequencing analysis provide insights into genome evolution and allelic imbalance in Pinus densiflora.

Haplotype-level allelic characterization facilitates research on the functional, evolutionary and breeding-related features of extremely large and complex plant genomes. We report a 21.7-Gb chromosome-level haplotype-resolved assembly in Pinus densiflora. We found genome rearrangements involving translocations and inversions between chromosomes 1 and 3 of Pinus species and a proliferation of specific long terminal repeat (LTR) retrotransposons (LTR-RTs) in P. densiflora. Evolutionary analyses illustrated that tandem and LTR-RT-mediated duplications led to an increment of transcription factor (TF) genes in P. densiflora. The haplotype sequence comparison showed allelic imbalances, including presence-absence variations of genes (PAV genes) and their functional contributions to flowering and abiotic stress-related traits in P. densiflora. Allele-aware resequencing analysis revealed PAV gene diversity across P. densiflora accessions. Our study provides insights into key mechanisms underlying the evolution of genome structure, LTR-RTs and TFs within the Pinus lineage as well as allelic imbalances and diversity across P. densiflora.

Template switching during DNA replication is a prevalent source of adaptive gene amplification.

Copy number variants (CNVs)-gains and losses of genomic sequences-are an important source of genetic variation underlying rapid adaptation and genome evolution. However, despite their central role in evolution little is known about the factors that contribute to the structure, size, formation rate, and fitness effects of adaptive CNVs. Local genomic sequences are likely to be an important determinant of these properties. Whereas it is known that point mutation rates vary with genomic location and local DNA sequence features, the role of genome architecture in the formation, selection, and the resulting evolutionary dynamics of CNVs is poorly understood. Previously, we have found that the GAP1 gene in Saccharomyces cerevisiae undergoes frequent and repeated amplification and selection under long-term experimental evolution in glutamine-limiting conditions. The GAP1 gene has a unique genomic architecture consisting of two flanking long terminal repeats (LTRs) and a proximate origin of DNA replication (autonomously replicating sequence, ARS), which are likely to promote rapid GAP1 CNV formation. To test the role of these genomic elements on CNV-mediated adaptive evolution, we performed experimental evolution in glutamine-limited chemostats using engineered strains lacking either the adjacent LTRs, ARS, or all elements. Using a CNV reporter system and neural network simulation-based inference (nnSBI) we quantified the formation rate and fitness effect of CNVs for each strain. We find that although GAP1 CNVs repeatedly form and sweep to high frequency in strains with modified genome architecture, removal of local DNA elements significantly impacts the rate and fitness effect of CNVs and the rate of adaptation. We performed genome sequence analysis to define the molecular mechanisms of CNV formation for 177 CNV lineages. We find that across all four strain backgrounds, between 26% and 80% of all GAP1 CNVs are mediated by Origin Dependent Inverted Repeat Amplification (ODIRA) which results from template switching between the leading and lagging strand during DNA synthesis. In the absence of the local ARS, a distal ARS can mediate CNV formation via ODIRA. In the absence of local LTRs, homologous recombination mechanisms still mediate gene amplification following de novo insertion of retrotransposon elements at the locus. Our study demonstrates the remarkable plasticity of the genome and reveals that template switching during DNA replication is a frequent source of adaptive CNVs.