Repeat Elements Research Articles

Genetic variation of the HIV-1 subtype C transmitted/founder viruses long terminal repeat elements and the impact on transcription activation potential and clinical disease outcomes.

A genetic bottleneck is a hallmark of HIV-1 transmission such that only very few viral strains, termed transmitted/founder (T/F) variants establish infection in a newly infected host. Phenotypic characteristics of these variants may determine the subsequent course of disease. The HIV-1 5' long terminal repeat (LTR) promoter drives viral gene transcription and is genetically identical to the 3' LTR. We hypothesized that HIV-1 subtype C (HIV-1C) T/F virus LTR genetic variation is a determinant of transcriptional activation potential and clinical disease outcome. The 3'LTR was amplified from plasma samples of 41 study participants acutely infected with HIV-1C (Fiebig stages I and V/VI). Paired longitudinal samples were also available at one year post-infection for 31 of the 41 participants. 3' LTR amplicons were cloned into a pGL3-basic luciferase expression vector, and transfected alone or together with Transactivator of transcription (tat) into Jurkat cells in the absence or presence of cell activators (TNF-α, PMA, Prostratin and SAHA). Inter-patient T/F LTR sequence diversity was 5.7% (Renge: 2-12) with subsequent intrahost viral evolution observed in 48.4% of the participants analyzed at 12 months post-infection. T/F LTR variants exhibited differential basal transcriptional activity, with significantly higher Tat-mediated transcriptional activity compared to basal (p<0.001). Basal and Tat-mediated T/F LTR transcriptional activity showed significant positive correlation with contemporaneous viral loads and negative correlation with CD4 T cell counts (p<0.05) during acute infection respectively. Furthermore, Tat-mediated T/F LTR transcriptional activity significanly correlated positively with viral load set point and viral load; and negatively with CD4 T cell counts at one year post infection (all p<0.05). Lastly, PMA, Prostratin, TNF-α and SAHA cell stimulation resulted in enhanced yet heterologous transcriptional activation of different T/F LTR variants. Our data suggest that T/F LTR variants may influence viral transcriptional activity, disease outcomes and sensitivity to cell activation, with potential implications for therapeutic interventions.

Contribution to the Knowledge of Dicranoptychini (Diptera, Tipuloidea, Limoniidae) in China, with the First Mitochondrial Genome of the Tribe and Its Phylogenetic Implications.

Dicranoptychini is a tribe in the subfamily Limoniinae (Diptera, Tipuloidea, and Limoniidae) and includes only the genus Dicranoptycha Osten Sacken, 1860. However, the species diversity of the tribe in China was seriously underestimated, and the taxonomic status of Dicranoptycha has long been controversial. In this study, types of Chinese Dicranoptycha species and specimens collected from several localities in China were examined, and the first mitochondrial (mt) genome of the tribe Dicranoptychini is presented. Two Dicranoptycha species, D. jiufengshana sp. nov. and D. shandongensis sp. nov., from China, are described and illustrated as new to science. A Palaearctic species, D. prolongata Alexander, 1938, is recorded in China for the first time. In addition, the complete mt genome of D. shandongensis sp. nov. is sequenced and annotated, indicating that it is a typical circular DNA molecule with a length of 16,157 bp and shows a similar gene order, nucleotide composition, and codon usage to mt genomes of other Tipuloidea species. The two pairs of repeat elements are found in its control region. Phylogenetic results confirm the sister-group relationship between Cylindrotomidae and Tipulidae, question the position of the genus Epiphragma Osten Sacken, 1860 in Limoniidae, and indicate that Dicranoptychini may be a basal lineage within Limoniinae.

Genome-wide landscapes of genes and repeatome reveal the genomic differences between the two subspecies of peanut (Arachis hypogaea)

Distribution and structural features of genes, repeat elements and transposable elements (TEs) were studied to identify the genomic differences between the two subspecies (ssp. hypogaea and ssp. fastigiata) of peanut (Arachis hypogaea L.). A total of 67128 predicted genes, 2738666 copies of TEs and 162361 tandem repeats from the reference genome of Tifrunner were employed for this study. Of the 67128 genes, 33622 were reading on the plus strand, while 33506 were traced on the minus strand. Though B03 had the highest number of genes (4524), A08 recorded the highest density (53 genes/Mb) of genes in the genome. Telomeric regions had the highest density of genes. The average length of the genes was 3971 bp with majority of the genes (39228) containing one to five exons. The gene Arahy.11.0DU9MH had the insertion of 28 different types of TEs, and was the longest gene in the peanut genome. A total of 15731 genes were monomorphic in terms of SNPs across 179 accessions, while 7401 genes showed polymorphism at one nucleotide, indicating very low allelic variation at these genes. Remaining 66% of the genes had two or more SNPs, and therefore showed relatively high allelic variation. Among the 101 unique types of TEs, the Retro elements (869279) followed by CACTA (272596) and Mu (250248) TEs were most predominant. Telomeric regions showed less density of TEs than the regions. On an average, each gene contained 1.8 copies of TEs, and 35706 genes did not have the insertion of any TEs. A QTL-Seq approach could identify 186 SNPs and 26 gene differences between the two subspecies of A. hypogaea. Two of the 26 genes showed allelic variation in terms of SNPs and TEs.

Unveiling the hidden players: The crucial role of transposable elements in the placenta and their potential contribution to pre-eclampsia

The human placenta is a vital connection between maternal and fetal tissues, allowing for the exchange of molecules and modulation of immune interactions during pregnancy. Interestingly, some of the placenta's unique functions can be attributed to transposable elements (TEs), which are DNA sequences that have mobilised into the genome. Co-option throughout mammalian evolution has led to the generation of TE-derived regulators and TE-derived genes, some of which are expressed in the placenta but silenced in somatic tissues. TE genes encompass both TE-derived genes with a repeat element in the coding region and TE-derived regulatory regions such as alternative promoters and enhancers. Placental-specific TE genes are known to contribute to the placenta's unique functions, and interestingly, they are also expressed in some cancers and share similar functions. There is evidence to support that aberrant activity of TE genes may contribute to placental pathologies, cancer and autoimmunity. In this review, we highlight the crucial roles of TE genes in placental function, and how their dysregulation may lead to pre-eclampsia, a common and dangerous placental condition. We provide a summary of the functional TE genes in the placenta to offer insight into their significance in normal and abnormal human development. Ultimately, this review highlights an opportunity for future research to investigate the potential dysregulation of TE genes in the development of placental pathologies such as pre-eclampsia. Further understanding of TE genes and their role in the placenta could lead to significant improvements in maternal and fetal health.

Characterization of repetitive DNA on the genome of the marsh rat Holochilus nanus (Cricetidae: Sigmodontinae).

Repetitive DNA are sequences repeated hundreds or thousands of times and an abundant part of eukaryotic genomes. SatDNA represents the majority of the repetitive sequences, followed by transposable elements. The species Holochilus nanus (HNA) belongs to the rodent tribe Oryzomyini, the most taxonomically diverse of Sigmodontinae subfamily. Cytogenetic studies on Oryzomyini reflect such diversity by revealing an exceptional range of karyotype variability. However, little is known about the repetitive DNA content and its involvement in chromosomal diversification of these species. In the search for a more detailed understanding about the composition of repetitive DNA on the genome of HNA and other species of Oryzomyini, we employed a combination of bioinformatic, cytogenetic and molecular techniques to characterize the repetitive DNA content of these species. RepeatExplorer analysis showed that almost half of repetitive content of HNA genome are composed by Long Terminal Repeats and a less significant portion are composed by Short Interspersed Nuclear Elements and Long Interspersed Nuclear Elements. RepeatMasker showed that more than 30% of HNA genome are composed by repetitive sequences, with two main waves of repetitive element insertion. It was also possible to identify a satellite DNA sequence present in the centromeric region of Oryzomyini species, and a repetitive sequence enriched on the long arm of HNA X chromosome. Also, comparative analysis between HNA genome with and without B chromosome did not evidence any repeat element enriched on the supernumerary, suggesting that B chromosome of HNA is composed by a fraction of repeats from all the genome.

Considering admixture when producing draft genomes: an example in North American ratsnakes (Pantherophis alleghaniensis/Pantherophis obsoletus).

The number of reference genomes of snakes lags behind several other vertebrate groups (e.g. birds and mammals). However, in the last two years, a concerted effort by researchers from around the world has produced new genomes of snakes representing members from several new families. Here, we present a high-quality, annotated genome of the central ratsnake (Pantherophis alleghaniensis), a member of the most diverse snake lineage, Colubroidea. Pantherophis alleghaniensis is found in the central part of the Nearctic, east of the Mississippi River. This genome was sequenced using 10X Chromium synthetic long reads and polished using Illumina short reads. The final genome assembly had an N50 of 21.82 Mb and an L50 of 22 scaffolds with a maximum scaffold length of 82.078 Mb. The genome is composed of 49.24% repeat elements dominated by long interspersed elements. We annotated this genome using transcriptome assemblies from 14 tissue types and recovered 28,368 predicted proteins. Finally, we estimated admixture proportions between two species of ratsnakes and discovered that this specimen is an admixed individual containing genomes from the western (Pantherophis obsoletus) and central ratsnakes (P. alleghaniensis). We discuss the importance of considering interspecific admixture in downstream approaches for inferring demography and phylogeny.

Low coverage sequencing provides insights into the key features of the nuclear and mitochondrial genomes of the Alligator Snapping Turtle Macrochelys temminckii

The alligator snapping turtle Macrochelys temminckii is a culturally, ecologically, and evolutionary relevant species of conservation concern. In this study, we conducted a genome survey of M. temminckii. Using a low-coverage short read sequencing strategy, this study estimated the genome size, repetitive genome content, annotated and quantified repetitive elements, assembled the 45S rRNA DNA operon, and characterized in detail the mitochondrial genome of M. temminckii. Using a k-mer strategy, the estimated haploid genome size varied between 3.77 and 3.19 Gbp, which is within the range previously reported for other representatives of the family Chelydridae. Repetitive genome content estimates using different k-mers (21 to 51) indicated that more than 75 % of the genome of M. temminckii comprised repetitive elements. Taking into account only annotated repetitive elements, the most common repetitive elements were classified as Class I - Long Interspersed Nuclear Element (LINE) which were more abundant than Class I – Penelope and Class I - Long Terminal Repeat (LTR) Ty3 mobile elements. Less abundant repeat element families in the nuclear genome of M. temminckii included Class I - DIRS mobile elements and Satellite DNA. The nuclear ribosomal operon was partially assembled into three contigs, one encoding the complete ssrRNA gene, a second encoding the complete 5.8S rRNA gene, and a third comprising the full lsrRNA gene. The AT-rich complete mitochondrial genome was 16,570 bp long. These new genomic resources are of utmost importance to aid in the development of conservation plans for this iconic freshwater turtle.

Differential effect of histone H3.3 depletion on retroviral repression in embryonic stem cells

BackgroundIntegration of retroviruses into the host genome can impair the genomic and epigenomic integrity of the cell. As a defense mechanism, epigenetic modifications on the proviral DNA repress retroviral sequences in mouse embryonic stem cells (ESC). Here, we focus on the histone 3 variant H3.3, which is abundant in active transcription zones, as well as centromeres and heterochromatinized repeat elements, e.g., endogenous retroviruses (ERV).ResultsTo understand the involvement of H3.3 in the epigenetic silencing of retroviral sequences in ESC, we depleted the H3.3 genes in ESC and transduced the cells with GFP-labeled MLV pseudovirus. This led to altered retroviral repression and reduced Trim28 recruitment, which consequently led to a loss of heterochromatinization in proviral sequences. Interestingly, we show that H3.3 depletion has a differential effect depending on which of the two genes coding for H3.3, H3f3a or H3f3b, are knocked out. Depletion of H3f3a resulted in a transient upregulation of incoming retroviral expression and ERVs, while the depletion of H3f3b did not have the same effect and repression was maintained. However, the depletion of both genes resulted in a stable activation of the retroviral promoter. These findings suggest that H3.3 is important for regulating retroviral gene expression in mouse ESC and provide evidence for a distinct function of the two H3.3 genes in this regulation. Furthermore, we show that Trim28 is needed for depositing H3.3 in retroviral sequences, suggesting a functional interaction between Trim28 recruitment and H3.3 loading.ConclusionsIdentifying the molecular mechanisms by which H3.3 and Trim28 interact and regulate retroviral gene expression could provide a deeper understanding of the fundamental processes involved in retroviral silencing and the general regulation of gene expression, thus providing new answers to a central question of stem cell biology.

How doctors build community and socialize into a clinical department through morning reports. A positioning theory study.

The morning report is one of the longest surviving hospital practices. Most studies of the morning report focus on the effectiveness of formal medical training, while focus on social and communicative aspects is rarer. This study explores the social interactions and communication in morning reports, examining the ways in which they contribute to the construction of professional identity and socialization into the community of the clinical department. We used a qualitative explorative design with video observations of morning reports. Our data consisted of 43 video-recorded observations (in all, 15.5 hours) from four different hospital departments in Denmark. These were analyzed using the theoretical framework of positioning theory. A key finding was that each department followed its own individual structure. This order was not articulated as such but played out implictly. Two alternative storylines unfolded in the elements of the morning report: 1) being equal members of the specialty and department, and 2) preserving the hierarchical community and its inherent positions. The morning report can be seen as playing an important role in community making. It unfolds as a "dance" of repeated elements in a complex collegial space. Within this complexity, the morning report is a space for positioning oneself and others as a collegial "we", i.e., equal members of a department and specialty, at the same time as "having a place" in a hierarchal community. Thus, morning reports contribute to developing professional identity and socialization into the medical community.

Telomere-to-Telomere genome assemblies of human-infecting Encephalitozoon species

BackgroundMicrosporidia are diverse spore forming, fungal-related obligate intracellular pathogens infecting a wide range of hosts. This diversity is reflected at the genome level with sizes varying by an order of magnitude, ranging from less than 3 Mb in Encephalitozoon species (the smallest known in eukaryotes) to more than 50 Mb in Edhazardia spp. As a paradigm of genome reduction in eukaryotes, the small Encephalitozoon genomes have attracted much attention with investigations revealing gene dense, repeat- and intron-poor genomes characterized by a thorough pruning of molecular functions no longer relevant to their obligate intracellular lifestyle. However, because no Encephalitozoon genome has been sequenced from telomere-to-telomere and since no methylation data is available for these species, our understanding of their overall genetic and epigenetic architectures is incomplete.MethodsIn this study, we sequenced the complete genomes from telomere-to-telomere of three human-infecting Encephalitozoon spp. —E. intestinalis ATCC 50506, E. hellem ATCC 50604 and E. cuniculi ATCC 50602— using short and long read platforms and leveraged the data generated as part of the sequencing process to investigate the presence of epigenetic markers in these genomes. We also used a mixture of sequence- and structure-based computational approaches, including protein structure prediction, to help identify which Encephalitozoon proteins are involved in telomere maintenance, epigenetic regulation, and heterochromatin formation.ResultsThe Encephalitozoon chromosomes were found capped by TTAGG 5-mer telomeric repeats followed by telomere associated repeat elements (TAREs) flanking hypermethylated ribosomal RNA (rRNA) gene loci featuring 5-methylcytosines (5mC) and 5-hemimethylcytosines (5hmC), themselves followed by lesser methylated subtelomeres and hypomethylated chromosome cores. Strong nucleotide biases were identified between the telomeres/subtelomeres and chromosome cores with significant changes in GC/AT, GT/AC and GA/CT contents. The presence of several genes coding for proteins essential to telomere maintenance, epigenetic regulation, and heterochromatin formation was further confirmed in the Encephalitozoon genomes.ConclusionAltogether, our results strongly support the subtelomeres as sites of heterochromatin formation in Encephalitozoon genomes and further suggest that these species might shutdown their energy-consuming ribosomal machinery while dormant as spores by silencing of the rRNA genes using both 5mC/5hmC methylation and facultative heterochromatin formation at these loci.

The mitochondrial genome of the bioluminescent fish Malacosteusniger Ayres, 1848 (Stomiidae, Actinopterygii) is large and complex, and contains an inverted-repeat structure.

We determined the complete mitogenome sequence of the bioluminescent fish Malacosteusniger using long-read sequencing technologies. The 21,263 bp mitogenome features a complex structure with two copies of a 1198-bp inverted-repeat and a region of 2616-bp containing alternating copies of 16 and 26 bp repeat elements. Whole mitogenome phylogenies inferred from both nucleotide and amino-acid datasets place M.niger among Melanostomiinae. The need for additional complete mitogenome sequences from the subfamily Malacosteinae is discussed.

ICOR: improving codon optimization with recurrent neural networks

BackgroundIn protein sequences—as there are 61 sense codons but only 20 standard amino acids—most amino acids are encoded by more than one codon. Although such synonymous codons do not alter the encoded amino acid sequence, their selection can dramatically affect the expression of the resulting protein. Codon optimization of synthetic DNA sequences is important for heterologous expression. However, existing solutions are primarily based on choosing high-frequency codons only, neglecting the important effects of rare codons. In this paper, we propose a novel recurrent-neural-network based codon optimization tool, ICOR, that aims to learn codon usage bias on a genomic dataset of Escherichia coli. We compile a dataset of over 7,000 non-redundant, high-expression, robust genes which are used for deep learning. The model uses a bidirectional long short-term memory-based architecture, allowing for the sequential context of codon usage in genes to be learned. Our tool can predict synonymous codons for synthetic genes toward optimal expression in Escherichia coli.ResultsWe demonstrate that sequential context achieved via RNN may yield codon selection that is more similar to the host genome. Based on computational metrics that predict protein expression, ICOR theoretically optimizes protein expression more than frequency-based approaches. ICOR is evaluated on 1,481 Escherichia coli genes as well as a benchmark set of 40 select DNA sequences whose heterologous expression has been previously characterized. ICOR’s performance is measured across five metrics: the Codon Adaptation Index, GC-content, negative repeat elements, negative cis-regulatory elements, and codon frequency distribution.ConclusionsThe results, based on in silico metrics, indicate that ICOR codon optimization is theoretically more effective in enhancing recombinant expression of proteins over other established codon optimization techniques. Our tool is provided as an open-source software package that includes the benchmark set of sequences used in this study.

Abstract 137: Epigenetic regulation of APOBEC3A mutagenesis and tumor evolution during targeted therapy in non-small cell lung cancer

Abstract Acquired drug resistance to even the most effective anti-cancer targeted therapies remains an unsolved clinical problem. Although many drivers of acquired drug resistance have been identified, the underlying molecular mechanisms shaping tumor evolution during treatment are incompletely understood. We recently demonstrated that lung cancer targeted therapies commonly used in the clinic induce the expression of cytidine deaminase APOBEC3A (A3A), leading to sustained mutagenesis in drug tolerant cancer cells persisting during therapy. Preventing therapy-induced A3A mutagenesis by gene deletion decreased the accumulation of chromosomal aberrations and delayed the emergence of drug resistance. Thus, we hypothesize that inhibition of A3A may represent a potential therapeutic strategy to halt the evolution of resistant clones and prevent acquired resistance to lung cancer targeted therapies. Understanding the molecular mechanism of A3A induction during targeted therapy is a crucial step to develop targeting A3A therapies.Using a targeted drug combination screen, we found that DNA methyltransferase (DNMT) inhibitors induce expression of APOBEC3A in non-small cell lung cancer cells, phenocopying the effects of targeted therapies. RNA-seq profiling revealed that both targeted therapies and DNMT inhibitors increase expression of non-coding repeat RNAs. Repeatome analysis identified distinct classes of non-codding RNAs including endogenous retrovirus elements (ERV) in drug tolerant persister cells. Activation of endogenous intracellular viral sensing pathways by exogenous nucleic acids mimicking viral infection induced APOBEC3A, and genetic perturbation of RIG-I or MAVS reduced APOBEC3A induced by targeted therapy. These findings suggest that epigenetic derepression of retroviral repeat elements may underly APOBEC3A mutagenic activity and tumor evolution during lung cancer targeted therapy. Citation Format: Hideko Isozaki, Ramin Sakhtemani, Naveed Nikpour, Susanna Monroe, Jessica Lin, Lecia Sequist, Zofia Piotrowska, Justin Gainor, Rémi Buisson, Michael Lawrence, Aaron Hata. Epigenetic regulation of APOBEC3A mutagenesis and tumor evolution during targeted therapy in non-small cell lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 137.

Deconstructed beetles: Bilayered composite materials produce green coloration with remarkably high near-infrared reflectance

Beetle elytra (hardened forewings) are a promising source of inspiration to develop or enhance the performance of human-fabricated composite materials. The structures responsible for optical properties in the ultra-violet to visible spectrum (300–700 nm) have been extensively characterised, but we have limited knowledge of optical properties and their physical origin in the near-infrared (NIR; 700–1700 nm). We examined the elytra of three species of green scarab beetles (Xylonichus eucalypti, Anoplognathus prasinus and Paraschizognathus olivaceus) with very high NIR reflectance. We manually separated layers in the elytra to disambiguate their contributions to the overall optical response. We show that unlike other scarabs, nanostructures within the cuticle layer do not produce notable reflectance. Instead, the cuticle resembles a pigment-based filter with 50% transmittance in the NIR and absorption in the visible spectrum contributing to the green appearance. Each species has a layer below the cuticle that appears white to the naked eye and produces broadband reflectance, particularly in the near-infrared; however, the structure of the white underlay differs markedly between the three species. In A. prasinus and P. olivaceus, the structure is disordered (no regular, repeated elements at optical length scales); whereas in Xylonichus eucalypti, the white underlay was notably thinner and comprised quasi-ordered hollow cylindrical structures embedded in a chitin matrix. We modelled the coherent scattering produced by this structure to demonstrate that it is responsible for broadband visible and NIR reflectance. We discuss biological implications and technological applications of the composite structure of beetle elytra.

Comparative analysis of the organelle genomes of three Rhodiola species provide insights into their structural dynamics and sequence divergences

BackgroundPlant organelle genomes are a valuable resource for evolutionary biology research, yet their genome architectures, evolutionary patterns and environmental adaptations are poorly understood in many lineages. Rhodiola species is a type of flora mainly distributed in highland habitats, with high medicinal value. Here, we assembled the organelle genomes of three Rhodiola species (R. wallichiana, R. crenulata and R. sacra) collected from the Qinghai-Tibet plateau (QTP), and compared their genome structure, gene content, structural rearrangements, sequence transfer and sequence evolution rates.ResultsThe results demonstrated the contrasting evolutionary pattern between plastomes and mitogenomes in three Rhodiola species, with the former possessing more conserved genome structure but faster evolutionary rates of sequence, while the latter exhibiting structural diversity but slower rates of sequence evolution. Some lineage-specific features were observed in Rhodiola mitogenomes, including chromosome fission, gene loss and structural rearrangement. Repeat element analysis shows that the repeats occurring between the two chromosomes may mediate the formation of multichromosomal structure in the mitogenomes of Rhodiola, and this multichromosomal structure may have recently formed. The identification of homologous sequences between plastomes and mitogenomes reveals several unidirectional protein-coding gene transfer events from chloroplasts to mitochondria. Moreover, we found that their organelle genomes contained multiple fragments of nuclear transposable elements (TEs) and exhibited different preferences for TEs insertion type. Genome-wide scans of positive selection identified one gene matR from the mitogenome. Since the matR is crucial for plant growth and development, as well as for respiration and stress responses, our findings suggest that matR may participate in the adaptive response of Rhodiola species to environmental stress of QTP.ConclusionThe study analyzed the organelle genomes of three Rhodiola species and demonstrated the contrasting evolutionary pattern between plastomes and mitogenomes. Signals of positive selection were detected in the matR gene of Rhodiola mitogenomes, suggesting the potential role of this gene in Rhodiola adaptation to QTP. Together, the study is expected to enrich the genomic resources and provide valuable insights into the structural dynamics and sequence divergences of Rhodiola species.

Transposable Element Interactions Shape the Ecology of the Deer Mouse Genome.

The genomic landscape of transposable elements (TEs) varies dramatically across species, with some TEs demonstrating greater success in colonizing particular lineages than others. In mammals, long interspersed nuclear element (LINE) retrotransposons are typically more common than any other TE. Here, we report an unusual genomic landscape of TEs in the deer mouse, Peromyscus maniculatus. In contrast to other previously examined mammals, long terminal repeat elements occupy more of the deer mouse genome than LINEs (11% and 10%, respectively). This pattern reflects a combination of relatively low LINE activity and a massive invasion of lineage-specific endogenous retroviruses (ERVs). Deer mouse ERVs exhibit diverse origins spanning the retroviral phylogeny suggesting they have been host to a wide range of exogenous retroviruses. Notably, we trace the origin of one ERV lineage, which arose ∼5-18 million years ago, to a close relative of feline leukemia virus, revealing inter-ordinal horizontal transmission. Several lineage-specific ERV subfamilies have very high copy numbers, with the top five most abundant accounting for ∼2% of the genome. We also observe a massive amplification of Kruppel-associated box domain-containing zinc finger genes, which likely control ERV activity and whose expansion may have been facilitated by ectopic recombination between ERVs. Finally, we find evidence that ERVs directly impacted the evolutionary trajectory of LINEs by outcompeting them for genomic sites and frequently disrupting autonomous LINE copies. Together, our results illuminate the genomic ecology that shaped the unique deer mouse TE landscape, shedding light on the evolutionary processes that give rise to variation in mammalian genome structure.

Complete sequences of six major histocompatibility complex haplotypes, including all the major MHC class II structures.

Accurate and comprehensive immunogenetic reference panels are key to the successful implementation of population-scale immunogenomics. The 5Mbp Major Histocompatibility Complex (MHC) is the most polymorphic region of the human genome and associated with multiple immune-mediated diseases, transplant matching and therapy responses. Analysis of MHC genetic variation is severely complicated by complex patterns of sequence variation, linkage disequilibrium and a lack of fully resolved MHC reference haplotypes, increasing the risk of spurious findings on analyzing this medically important region. Integrating Illumina, ultra-long Nanopore, and PacBio HiFi sequencing as well as bespoke bioinformatics, we completed five of the alternative MHC reference haplotypes of the current (GRCh38/hg38) build of the human reference genome and added one other. The six assembled MHC haplotypes encompass the DR1 and DR4 haplotype structures in addition to the previously completed DR2 and DR3, as well as six distinct classes of the structurally variable C4 region. Analysis of the assembled haplotypes showed that MHC class II sequence structures, including repeat element positions, are generally conserved within the DR haplotype supergroups, and that sequence diversity peaks in three regions around HLA-A, HLA-B+C, and the HLA class II genes. Demonstrating the potential for improved short-read analysis, the number of proper read pairs recruited to the MHC was found to be increased by 0.06%-0.49% in a 1000 Genomes Project read remapping experiment with seven diverse samples. Furthermore, the assembled haplotypes can serve as references for the community and provide the basis of a structurally accurate genotyping graph of the complete MHC region.

Identification and Characterization of the HERV-K (HML-8) Group of Human Endogenous Retroviruses in the Genome.

Human endogenous retroviruses (HERVs) can be vertically transmitted in a Mendelian fashion, are stably maintained in the human genome, and are estimated to constitute ∼8% of the genome. HERVs affect human physiology and pathology through their provirus-encoded protein or long terminal repeat (LTR) element effect. Characterization of the genomic distribution is an essential step to understanding the relationships between endogenous retrovirus expression and diseases. However, the poor characterization of human MMTV-like (HML)-8 prevents a detailed understanding of the regulation of the expression of this family in humans and its impact on the host genome. In light of this, the definition of an accurate and updated HERV-K HML-8 genomic map is urgently needed. In this study, we report the results of a comprehensive analysis of HERV-K HML-8 sequence presence and distribution within the human genome and hominoids, with a detailed description of the different structural and phylogenetic aspects characterizing the group. A total of 40 proviruses and 5 solo LTR elements for human were characterized, which included a detailed description of provirus structure, integration time, potentially regulated genes, transcription factor-binding sites, and primer-binding site features. Besides, 9 chimpanzee sequences, 8 gorilla sequences, and 10 orangutan sequences belonging to the HML-8 subgroup were identified. The integration time results showed that the HML-8 elements were integrated into the primate lineage around 35 and 42 million years ago (mya), during primates evolutionary speciation. Overall, the results clarified the composition of the HML-8 groups, providing an exhaustive background for subsequent functional studies.

Structures of Ric-8B in complex with Gα protein folding clients reveal isoform specificity mechanisms

Mammalian Ric-8 proteins act as chaperones to regulate the cellular abundance of heterotrimeric G protein α subunits. The Ric-8A isoform chaperones Gαi/o, Gα12/13, and Gαq/11 subunits, while Ric-8B acts on Gαs/olf subunits. Here, we determined cryoelectron microscopy (cryo-EM) structures of Ric-8B in complex with Gαs and Gαolf, revealing isoform differences in the relative positioning and contacts between the C-terminal α5 helix of Gα within the concave pocket formed by Ric-8 α-helical repeat elements. Despite the overall architectural similarity with our earlier structures of Ric-8A complexed to Gαq and Gαi1, Ric-8B distinctly accommodates an extended loop found only in Gαs/olf proteins. The structures, along with results from Ric-8 protein thermal stability assays and cell-based Gαolf folding assays, support a requirement for the Gα C-terminal region for binding specificity, and highlight that multiple structural elements impart specificity for Ric-8/G protein binding.

Towards routine chromosome-scale haplotype-resolved reconstruction in cancer genomics

Cancer genomes are highly complex and heterogeneous. The standard short-read sequencing and analytical methods are unable to provide the complete and precise base-level structural variant landscape of cancer genomes. In this work, we apply high-resolution long accurate HiFi and long-range Hi-C sequencing to the melanoma COLO829 cancer line. Also, we develop an efficient graph-based approach that processes these data types for chromosome-scale haplotype-resolved reconstruction to characterise the cancer precise structural variant landscape. Our method produces high-quality phased scaffolds on the chromosome level on three healthy samples and the COLO829 cancer line in less than half a day even in the absence of trio information, outperforming existing state-of-the-art methods. In the COLO829 cancer cell line, here we show that our method identifies and characterises precise somatic structural variant calls in important repeat elements that were missed in short-read-based call sets. Our method also finds the precise chromosome-level structural variant (germline and somatic) landscape with 19,956 insertions, 14,846 deletions, 421 duplications, 52 inversions and 498 translocations at the base resolution. Our simple pstools approach should facilitate better personalised diagnosis and disease management, including predicting therapeutic responses.