Genome Evolution Research Articles

Establishing CRISPRi for Programmable Gene Repression and Genome Evolution in Cupriavidus necator.

Cupriavidus necator H16 is a "Knallgas" bacterium with the ability to utilize various carbon sources and has been employed as a versatile microbial cell factory to produce a wide range of value-added compounds. However, limited genome engineering, especially gene regulation methods, has constrained its full potential as a microbial production platform. The advent of CRISPR/Cas9 technology has shown promise in addressing this limitation. Here, we developed an optimized CRISPR interference (CRISPRi) system for gene repression in C. necator by expressing a codon-optimized deactivated Cas9 (dCas9) and appropriate single guide RNAs (sgRNAs). CRISPRi was proven to be a programmable and controllable tool and could successfully repress both exogenous and endogenous genes. As a case study, we decreased the accumulation of polyhydroxyalkanoate (PHB) via CRISPRi and rewired the carbon fluxes to the synthesis of lycopene. Additionally, by disturbing the expression of DNA mismatch repair gene mutS with CRISPRi, we established CRISPRi-Mutator for genome evolution, rapidly generating mutant strains with enhanced hydrogen peroxide tolerance and robustness in microbial electrosynthesis (MES) system. Our work provides an efficient CRISPRi toolkit for advanced genetic manipulation and optimization of C. necator cell factories for diverse biotechnology applications.

Comparative Analysis of Chloroplast Genome of Meconopsis (Papaveraceae) Provides Insights into Their Genomic Evolution and Adaptation to High Elevation.

The Meconopsis species are widely distributed in the Qinghai-Tibet Plateau, Himalayas, and Hengduan Mountains in China, and have high medicinal and ornamental value. The high diversity of plant morphology in this genus poses significant challenges for species identification, given their propensity for highland dwelling, which makes it a question worth exploring how they cope with the harsh surroundings. In this study, we recently generated chloroplast (cp) genomes of two Meconopsis species, Meconopsis paniculata (M. paniculata) and M. pinnatifolia, and compared them with those of ten Meconopsis cp genomes to comprehend cp genomic features, their phylogenetic relationships, and what part they might play in plateau adaptation. These cp genomes shared a great deal of similarities in terms of genome size, structure, gene content, GC content, and codon usage patterns. The cp genomes were between 151,864 bp and 154,997 bp in length, and contain 133 predictive genes. Through sequence divergence analysis, we identified three highly variable regions (trnD-psbD, ccsA-ndhD, and ycf1 genes), which could be used as potential markers or DNA barcodes for phylogenetic analysis. Between 22 and 38 SSRs and some long repeat sequences were identified from 12 Meconopsis species. Our phylogenetic analysis confirmed that 12 species of Meconopsis clustered into a monophyletic clade in Papaveraceae, which corroborated their intrageneric relationships. The results indicated that M. pinnatifolia and M. paniculata are sister species in the phylogenetic tree. In addition, the atpA and ycf2 genes were positively selected in high-altitude species. The functions of these two genes might be involved in adaptation to the extreme environment in the cold and low CO2 concentration conditions at the plateau.

Intracellular signaling in proto-eukaryotes evolves to alleviate regulatory conflicts of endosymbiosis.

The complex eukaryotic cell resulted from a merger between simpler prokaryotic cells, yet the role of the mitochondrial endosymbiosis with respect to other eukaryotic innovations has remained under dispute. To investigate how the regulatory challenges associated with the endosymbiotic state impacted genome and network evolution during eukaryogenesis, we study a constructive computational model where two simple cells are forced into an obligate endosymbiosis. Across multiple in silico evolutionary replicates, we observe the emergence of different mechanisms for the coordination of host and symbiont cell cycles, stabilizing the endosymbiotic relationship. In most cases, coordination is implicit, without signaling between host and symbiont. Signaling only evolves when there is leakage of regulatory products between host and symbiont. In the fittest evolutionary replicate, the host has taken full control of the symbiont cell cycle through signaling, mimicking the regulatory dominance of the nucleus over the mitochondrion that evolved during eukaryogenesis.

Extraction and selection of high-molecular-weight DNA for long-read sequencing from Chlamydomonas reinhardtii.

Recent advances in long-read sequencing technologies have enabled the complete assembly of eukaryotic genomes from telomere to telomere by allowing repeated regions to be fully sequenced and assembled, thus filling the gaps left by previous short-read sequencing methods. Furthermore, long-read sequencing can also help characterizing structural variants, with applications in the fields of genome evolution or cancer genomics. For many organisms, the main bottleneck to sequence long reads remains the lack of robust methods to obtain high-molecular-weight (HMW) DNA. For this purpose, we developed an optimized protocol to extract DNA suitable for long-read sequencing from the unicellular green alga Chlamydomonas reinhardtii, based on CTAB/phenol extraction followed by a size selection step for long DNA molecules. We provide validation results for the extraction protocol, as well as statistics obtained with Oxford Nanopore Technologies sequencing.

Oligo-barcode illuminates holocentric karyotype evolution in Rhynchospora (Cyperaceae).

Holocentric karyotypes are assumed to rapidly evolve through chromosome fusions and fissions due to the diffuse nature of their centromeres. Here, we took advantage of the recent availability of a chromosome-scale reference genome for Rhynchospora breviuscula, a model species of this holocentric genus, and developed the first set of oligo-based barcode probes for a holocentric plant. These probes were applied to 13 additional species of the genus, aiming to investigate the evolutionary dynamics driving the karyotype evolution in Rhynchospora. The two sets of probes were composed of 27,392 (green) and 23,968 (magenta) oligonucleotides (45-nt long), and generated 15 distinct FISH signals as a unique barcode pattern for the identification of all five chromosome pairs of the R. breviuscula karyotype. Oligo-FISH comparative analyzes revealed different types of rearrangements, such as fusions, fissions, putative inversions and translocations, as well as genomic duplications among the analyzed species. Two rounds of whole genome duplication (WGD) were demonstrated in R. pubera, but both analyzed accessions differed in the complex chain of events that gave rise to its large, structurally diploidized karyotypes with 2n = 10 or 12. Considering the phylogenetic relationships and divergence time of the species, the specificity and synteny of the probes were maintained up to species with a divergence time of ~25 My. However, karyotype divergence in more distant species hindered chromosome mapping and the inference of specific events. This barcoding system is a powerful tool to study chromosomal variations and genomic evolution in holocentric chromosomes of Rhynchospora species.

Telomere-to-telomere gapless genome assembly of the Chinese sea bass (Lateolabrax maculatus)

Chinese sea bass (Lateolabrax maculatus) is a highly sought-after commercial seafood species in Asian regions due to its excellent nutritional value. With the rapid advancement of bioinformatics, higher standards for genome analysis compared to previously published reference genomes are now necessary. This study presents a gapless assembly of the Chinese sea bass genome, which has a length of 632.75 Mb. The sequences were assembled onto 24 chromosomes with a coverage of over 99% (626.61 Mb), and telomeres were detected on 34 chromosome ends. Analysis using Merqury indicated a high level of accuracy, with an average consensus quality value of 54.25. The ONT ultralong and PacBio HiFi data were aligned with the assembly using minimap2, resulting in a mapping rate of 99.9%. The study also identified repeating elements in 20.90% (132.25 Mb) of the genome and inferred 22,014 protein-coding genes. These results establish meaningful groundwork for exploring the evolution of the Chinese sea bass genome and advancing molecular breeding techniques.

The chromosome-scale genome and population genomics reveal the adaptative evolution of Populus pruinosa to desertification environment.

The Populus pruinosa is a relic plant that has managed to survive in extremely harsh desert environments. Owing to intensifying global warming and desertification, research into ecological adaptation and speciation of P. pruinosa has attracted considerable interest, but the lack of a chromosome-scale genome has limited adaptive evolution research. Here, a 521.09Mb chromosome-level reference genome of P. pruinosa was reported. Genome evolution and comparative genomic analysis revealed that tandemly duplicated genes and expanded gene families in P. pruinosa contributed to adaptability to extreme desert environments (especially high salinity and drought). The long terminal repeat retrotransposons (LTR-RTs) inserted genes in the gene body region might drive the adaptive evolution of P. pruinosa and species differentiation in saline-alkali desert environments. We recovered genetic differentiation in the populations of the northern Tianshan Mountain and southern Tianshan Mountain through whole-genome resequencing of 156 P. pruinosa individuals from 25 populations in China. Further analyses revealed that precipitation drove the local adaptation of P. pruinosa populations via some genetic sites, such as MAG2-interacting protein 2 (MIP2) and SET domain protein 25 (SDG25). This study will provide broad implications for adaptative evolution and population studies by integrating internal genetic and external environmental factors in P. pruinosa.

Rare subclonal sequencing of breast cancers indicates putative metastatic driver mutations are predominately acquired after dissemination

BackgroundEvolutionary models of breast cancer progression differ on the extent to which metastatic potential is pre-encoded within primary tumors. Although metastatic recurrences often harbor putative driver mutations that are not detected in their antecedent primary tumor using standard sequencing technologies, whether these mutations were acquired before or after dissemination remains unclear.MethodsTo ascertain whether putative metastatic driver mutations initially deemed specific to the metastasis by whole exome sequencing were, in actuality, present within rare ancestral subclones of the primary tumors from which they arose, we employed error-controlled ultra-deep sequencing (UDS-UMI) coupled with FFPE artifact mitigation by uracil-DNA glycosylase (UDG) to assess the presence of 132 “metastasis-specific” mutations within antecedent primary tumors from 21 patients. Maximum mutation detection sensitivity was ~1% of primary tumor cells. A conceptual framework was developed to estimate relative likelihoods of alternative models of mutation acquisition.ResultsThe ancestral primary tumor subclone responsible for seeding the metastasis was identified in 29% of patients, implicating several putative drivers in metastatic seeding including LRP5 A65V and PEAK1 K140Q. Despite this, 93% of metastasis-specific mutations in putative metastatic driver genes remained undetected within primary tumors, as did 96% of metastasis-specific mutations in known breast cancer drivers, including ERRB2 V777L, ESR1 D538G, and AKT1 D323H. Strikingly, even in those cases in which the rare ancestral subclone was identified, 87% of metastasis-specific putative driver mutations remained undetected. Modeling indicated that the sequential acquisition of multiple metastasis-specific driver or passenger mutations within the same rare subclonal lineage of the primary tumor was highly improbable.ConclusionsOur results strongly suggest that metastatic driver mutations are sequentially acquired and selected within the same clonal lineage both before, but more commonly after, dissemination from the primary tumor, and that these mutations are biologically consequential. Despite inherent limitations in sampling archival primary tumors, our findings indicate that tumor cells in most patients continue to undergo clinically relevant genomic evolution after their dissemination from the primary tumor. This provides further evidence that metastatic recurrence is a multi-step, mutation-driven process that extends beyond primary tumor dissemination and underscores the importance of longitudinal tumor assessment to help guide clinical decisions.

Temporal Dynamics and Genomic Landscape of SARS-CoV-2 After Four Years of Evolution.

Introduction Since its emergence, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has undergone extensive genomic evolution, impacting public health policies, diagnosis, medication, and vaccine development. This study leverages advanced bioinformatics to assess the virus's temporal and regional genomic evolution from December 2019 to October 2023. Methods Our analysis incorporates 16,575 complete SARS-CoV-2 sequences collected from 214 countries. These samples were comparatively analyzed, with a detailed characterization of nucleic mutations, lineages, distribution, and evolutionary patterns during each year, using the Wuhan-Hu-1 strain as the reference. Results Our analysis has identified a total of 21,580 mutations that we classified into transient mutations, which diminished over time, and persistent mutations with steadily increasing frequencies. This mutation landscape led to a notable surge in the evolutionary rate, rising from 13 mutations per sample in 2020 to 96 by 2023, with minor geographic variations. The phylogenetic analysis unveiled three distinct evolutionary branches, each representing unique viral evolution pathways. These lineages exhibited a tendency for a reduced duration of dominance with a shortening prevalence period over time, as dominant strains were consistently replaced by more fit variants. Notably, the emergence of the Alpha and Delta variants in 2021 was followed by the subsequent dominance of Omicron clade variants that have branched into several recombinant variants in 2022, marking a significant shift in the viral landscape. Conclusion This study sheds light on the dynamic nature of SARS-CoV-2 evolution, emphasizing the importance of continuous and vigilant genomic surveillance. The dominance of recombinant lineages, coupled with the disappearance of local variants, underscores the virus's adaptability.

Genome-Wide Identification and Characterization of the TIFY Gene Family and Their Expression Patterns in Response to MeJA and Aluminum Stress in Centipedegrass (Eremochloa ophiuroides)

The TIFY family is a group of novel plant-specific transcription factors involved in plant development, signal transduction, and responses to stress and hormones. TIFY genes have been found and functionally characterized in a number of plant species. However, there is no information about this family in warm-season grass plants. The current study identified 24 TIFY genes in Eremochloa ophiuroides, a well-known perennial warm-season grass species with a high tolerance to aluminum toxicity and good adaptability to the barren acidic soils. All of the 24 EoTIFYs were unevenly located on six out of nine chromosomes and could be classified into two subfamilies (ZIM/ZML and JAZ), consisting of 3 and 21 genes, respectively, with the JAZ subfamily being further divided into five subgroups (JAZ I to JAZ V). The amino acids of 24 EoTIFYs showed apparent differences between the two subfamilies based on the analysis of gene structures and conserved motifs. MCScanX analysis revealed the tandem duplication and segmental duplication of several EoTIFY genes occurred during E. ophiuroides genome evolution. Syntenic analyses of TIFY genes between E. ophiuroides and other five plant species (including A. thaliana, O. sativa, B. distachyon, S. biocolor, and S. italica) provided valuable clues for understanding the potential evolution of the EoTIFY family. qRT-PCR analysis revealed that EoTIFY genes exhibited different spatial expression patterns in different tissues. In addition, the expressions of EoTIFY genes were highly induced by MeJA and all of the EoTIFY family members except for EoJAZ2 displayed upregulated expression by MeJA. Ten EoTIFY genes (EoZML1, EoZML1, EoJAZ1, EoJAZ3, EoJAZ5, EoJAZ6, EoJAZ8, EoJAZ9, EoJAZ10, and EoJAZ21) were observed to be highly expressed under both exogenous MeJA treatment and aluminum stress, respectively. These results suggest that EoTIFY genes play a role in the JA-regulated pathway of plant growth and aluminum resistance as well. The results of this study laid a foundation for further understanding the function of TIFY genes in E. ophiuroides, and provided useful information for future aluminum tolerance related breeding and gene function research in warm-season grass plants.

Transmission restriction and genomic evolution co-shape the genetic diversity patterns of influenza A virus

Influenza A virus (IAV) shows an extensive host range and rapid genomic variations, leading to continuous emergence of novel viruses with significant antigenic variations and the potential for cross-species transmission. This causes global pandemics and seasonal flu outbreaks, posing sustained threats worldwide. Thus, studying all IAVs' evolutionary patterns and underlying mechanisms is crucial for effective prevention and control. We developed FluTyping to identify IAV genotypes, to explore overall genetic diversity patterns and their restriction factors. FluTyping groups isolates based on genetic distance and phylogenetic relationships using entire genomes, enabling identification of each isolate's genotype. Three distinct genetic diversity patterns were observed: one genotype domination pattern comprising only H1N1 and H3N2 seasonal influenza subtypes, multi-genotypes co-circulation pattern including majority avian influenza subtypes and swine influenza H1N2, and hybrid-circulation pattern involving H7N9 and three H5 subtypes of influenza viruses. Furthermore, the IAVs in multi-genotypes co-circulation pattern showed region-specific dominant genotypes, implying the restriction of virus transmission is a key factor contributing to distinct genetic diversity patterns, and the genomic evolution underlying different patterns showed more influenced by host-specific factors. In summary, a comprehensive picture of the evolutionary patterns of overall IAVs is provided by the FluTyping's identified genotypes, offering important theoretical foundations for future prevention and control of these viruses.

Abstract PR012: Evolutionary dynamics of whole-genome doubling in ovarian cancer

Abstract Chromosomal instability (CIN) is a major driver of tumor progression and treatment resistance in many cancers. CIN is characterized by ongoing chromosome missegregation, segmental aneuploidy and whole-genome doubling (WGD), generating copy number heterogeneity that provides a substrate for natural selection. Although well characterized in model systems, the evolutionary dynamics and genomic impact of WGD is underexplored in primary tumors and metastases. Newly developed single-cell genomics assays now allow for the measurement of the whole genomes of single cells obtained from patient tissues, thus providing a direct view into the evolution of tumor genomes undergoing WGD. We performed multi-region single-cell whole genome sequencing (scWGS) in a cohort of patients with high-grade serous ovarian cancer (HGSOC), an archetype of high CIN cancers. We profiled a total of 70,644 single cells from 63 pre-treatment samples obtained for 36 HGSOC patients. For each patient, we computed total and allele-specific copy number per cell, and inferred the evolutionary history of each patient’s tumor, including clonal and subclonal WGD events and cell-specific whole-chromosome or arm-level copy number changes indicative of missegregations. To determine the role of CIN-associated micronuclei (MN) formation, we used high-resolution cGAS immunofluorescence (IF) microscopy to profile FFPE tissues site-matched to scWGS, and quantify spatially-resolved rates of MN formation and rupture across 120 samples. A majority of tumors were characterized by early clonal WGD (61%, 22/36 cases). The remaining patients exhibited subclonal WGD subpopulations averaging 10% of cancer cells in total. To understand the interplay between recent WGD and CIN across the cohort, we subset all patients cells into non-WGD, subclonal WGD, and clonal WGD and computed mis-segregation rates for these subsets. Non-WGD cells showed the lowest levels of mis-segregation. Subclonal WGD exhibited the highest CIN rates, including mis-segregation and multipolar cells indicative of CIN following a recent WGD event. Clonal WGD exhibited intermediate levels of mis-segregation indicative of stabilization following a historical WGD event. CIN following WGD often resulted in a reduction in ploidy relative to the ancestral WGD cell. Micronucleus rupture was ubiquitously detected in IF across HGSOC tumors, and MN rates varied widely between patients and across primary and metastatic sites. Micronucleus rupture commonly resulted from the impact of WGD, which increased sequestration of missegregated chromosomes into MN. Through multimodal single cell measurements of HGSOC tumors, we were able to characterize genomic complexity from individual cell division errors in HGSOC, and shed light into how the transient processes governing ovarian cancer cell evolution impact the genomes of cell populations. Citation Format: Ignacio Vázquez-García, Matthew Myers, Marc J. Williams, Samuel S. Freeman, Jun Li, Hongyu Shi, Gryte Satas, Nicholas Ceglia, Florian Uhlitz, MSK SPECTRUM Consortium, Britta Weigelt, Dmitriy Zamarin, Andrew McPherson, Samuel F. Bakhoum, Sohrab P. Shah. Evolutionary dynamics of whole-genome doubling in ovarian cancer [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Translating Cancer Evolution and Data Science: The Next Frontier; 2023 Dec 3-6; Boston, Massachusetts. Philadelphia (PA): AACR; Cancer Res 2024;84(3 Suppl_2):Abstract nr PR012.

WITHDRAWN: Chloroplast genome evolution of Berberis (Berberidaceae): Implications for phylogeny and metabarcoding

Berberidis Radix (Sankezhen), a typical multi-origin Chinese medicinal material, originates from the dried roots of the Berberis genus plants and is used to treat a range of ailments. These species have similar morphology potentially leading to misidentification that can impact the medicine's efficacy. Therefore, it is imperative to develop suitable molecular markers to identify medicinal species. Furthermore, discrepancies exist in the taxonomy of the Berberis genus. In the current study, we de novo assembled the chloroplast genomes of six Berberis species (Berberis julianae, Berberis tsienii, Berberis pruinosa, Berberis thunbergii, Berberis poiretii, and Berberis wilsoniae) that commonly constitute Berberidis Radix and compared them with previously reported genomes. Our comparative analysis revealed similarities in genome structure, relative synonymous codon usage, amino acid frequency, repeats, and substitutions. Higher synonymous substitutions, indicative of predominant purifying selection on protein-coding genes, were observed compared to non-synonymous substitutions. However, positive selection was identified for six genes across 29 Berberis species—accD, matK, ndhD, rbcL, ycf1 and ycf2—highlighting their potential roles in adaptive responses to specific environmental conditions within the genus. Inverted repeats expansion and contraction affected the rate of mutations and is associated with the phylogenetic classification of Berberis. Our phylogenetic analysis supported the division of the Berberis complex into four genera, which corroborates the findings of previous studies involving extensive sampling. We identified the ndhD-ccsA region as the most polymorphic region and applied this region to Chinese patent medicines containing Berberidis Radix through metabarcoding. Metabarcoding analysis confirmed that five Berberis species commonly constitute Berberidis Radix in Chinese patent medicines. In conclusion, this study provides insight into the molecular evolution of the chloroplast genome and the phylogeny of the Berberis genus. In addition, through metabarcoding, provides insight into the species composition of Berberidis Radix in Chinese patent medicines.

Abstract B029: Centuries of genome instability and evolution in soft-shell clam (Mya arenaria) transmissible cancer

Abstract Transmissible cancers are infectious parasitic clones that metastasize to new hosts, living past the death of the founder animal in which the cancer initiated. Using genomic and transcriptomic analyses, we investigated the evolutionary history of a recently identified transmissible cancer lineage that has spread though the soft-shell clam (Mya arenaria) population. We first assembled a chromosome-scale soft-shell clam reference genome and characterized somatic mutations in the cancer lineage, discovering a novel mutational signature that estimates the lineage to be >200 years old and observing a wide variety of mutation types indicative of an extremely unstable cancer genome. We next quantified gene expression, observing differential expression of stereotypical cancer hallmarks in addition to key pathways that may have facilitated the cancer’s ability to survive seawater transfer to new hosts and evade immune rejection. Finally, we developed genetic assays to track bivalve transmissible neoplasia in the soft-shell clam and identify a novel lineage in the Baltic clam (Macoma balthica). Taken together, this study reveals the long-term survival of an invertebrate cancer lineage and identifies adaptive mechanisms to evade physical and immune barriers to cancer transmission, which may have been facilitated in part by the adaptive potential of its unstable genome. Citation Format: Samuel F.M. Hart, Michael J. Metzger. Centuries of genome instability and evolution in soft-shell clam (Mya arenaria) transmissible cancer [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Translating Cancer Evolution and Data Science: The Next Frontier; 2023 Dec 3-6; Boston, Massachusetts. Philadelphia (PA): AACR; Cancer Res 2024;84(3 Suppl_2):Abstract nr B029.

Whole-genome sequencing and evolutionary analysis of the wild edible mushroom, Morchella eohespera.

Morels (Morchella, Ascomycota) are an extremely desired group of edible mushrooms with worldwide distribution. Morchella eohespera is a typical black morel species, belonging to the Elata clade of Morchella species. The biological and genetic studies of this mushroom are rare, largely hindering the studies of molecular breeding and evolutionary aspects. In this study, we performed de novo sequencing and assembly of the M. eohespera strain m200 genome using the third-generation nanopore sequencing platform. The whole-genome size of M. eohespera was 53.81 Mb with a contig N50 of 1.93 Mb, and the GC content was 47.70%. A total of 9,189 protein-coding genes were annotated. Molecular dating showed that M. eohespera differentiated from its relative M. conica at ~19.03 Mya (million years ago) in Burdigalian. Evolutionary analysis showed that 657 gene families were contracted and 244 gene families expanded in M. eohespera versus the related morel species. The non-coding RNA prediction results showed that there were 336 tRNAs, 76 rRNAs, and 45 snRNAs in the M. eohespera genome. Interestingly, there was a high degree of repetition (20.93%) in the M. eohespera genome, and the sizes of long interspersed nuclear elements, short interspersed nuclear elements, and long terminal repeats were 0.83 Mb, 0.009 Mb, and 4.56 Mb, respectively. Additionally, selection pressure analysis identified that a total of 492 genes in the M. eohespera genome have undergone signatures of positive selection. The results of this study provide new insights into the genome evolution of M. eohespera and lay the foundation for in-depth research into the molecular biology of the genus Morchella in the future.

Transposon dynamics in the emerging oilseed crop Thlaspi arvense.

Genome evolution is partly driven by the mobility of transposable elements (TEs) which often leads to deleterious effects, but their activity can also facilitate genetic novelty and catalyze local adaptation. We explored how the intraspecific diversity of TE polymorphisms might contribute to the broad geographic success and adaptive capacity of the emerging oil crop Thlaspi arvense (field pennycress). We classified the TE inventory based on a high-quality genome assembly, estimated the age of retrotransposon TE families and comprehensively assessed their mobilization potential. A survey of 280 accessions from 12 regions across the Northern hemisphere allowed us to quantify over 90,000 TE insertion polymorphisms (TIPs). Their distribution mirrored the genetic differentiation as measured by single nucleotide polymorphisms (SNPs). The number and types of mobile TE families vary substantially across populations, but there are also shared patterns common to all accessions. Ty3/Athila elements are the main drivers of TE diversity in T. arvense populations, while a single Ty1/Alesia lineage might be particularly important for transcriptome divergence. The number of retrotransposon TIPs is associated with variation at genes related to epigenetic regulation, including an apparent knockout mutation in BROMODOMAIN AND ATPase DOMAIN-CONTAINING PROTEIN 1 (BRAT1), while DNA transposons are associated with variation at the HSP19 heat shock protein gene. We propose that the high rate of mobilization activity can be harnessed for targeted gene expression diversification, which may ultimately present a toolbox for the potential use of transposition in breeding and domestication of T. arvense.

Comparative genomics of Japanese encephalitis virus shows low rates of recombination and a small subset of codon positions under episodic diversifying selection.

Orthoflavivirus japonicum (JEV) is the dominant cause of viral encephalitis in the Asian region with 100,000 cases and 25,000 deaths reported annually. The genome is comprised of a single polyprotein that encodes three structural and seven non-structural proteins. We collated a dataset of 349 complete genomes from a number of public databases, and analysed the data for recombination, evolutionary selection and phylogenetic structure. There are low rates of recombination in JEV, subsequently recombination is not a major evolutionary force shaping JEV. We found a strong overall signal of purifying selection in the genome, which is the main force affecting the evolutionary dynamics in JEV. There are also a small number of genomic sites under episodic diversifying selection, especially in the envelope protein and non-structural proteins 3 and 5. Overall, these results support previous analyses of JEV evolutionary genomics and provide additional insight into the evolutionary processes shaping the distribution and adaptation of this important pathogenic arbovirus.

ContScout: sensitive detection and removal of contamination from annotated genomes

Contamination of genomes is an increasingly recognized problem affecting several downstream applications, from comparative evolutionary genomics to metagenomics. Here we introduce ContScout, a precise tool for eliminating foreign sequences from annotated genomes. It achieves high specificity and sensitivity on synthetic benchmark data even when the contaminant is a closely related species, outperforms competing tools, and can distinguish horizontal gene transfer from contamination. A screen of 844 eukaryotic genomes for contamination identified bacteria as the most common source, followed by fungi and plants. Furthermore, we show that contaminants in ancestral genome reconstructions lead to erroneous early origins of genes and inflate gene loss rates, leading to a false notion of complex ancestral genomes. Taken together, we offer here a tool for sensitive removal of foreign proteins, identify and remove contaminants from diverse eukaryotic genomes and evaluate their impact on phylogenomic analyses.

Comparative genomic analyses provide new insights into evolutionary history and conservation genomics of gorillas

Genome sequencing is a powerful tool to understand species evolutionary history, uncover genes under selection, which could be informative of local adaptation, and infer measures of genetic diversity, inbreeding and mutational load that could be used to inform conservation efforts. Gorillas, critically endangered primates, have received considerable attention and with the recently sequenced Bwindi mountain gorilla population, genomic data is now available from all gorilla subspecies and both mountain gorilla populations. Here, we reanalysed this rich dataset with a focus on evolutionary history, local adaptation and genomic parameters relevant for conservation. We estimate a recent split between western and eastern gorillas of 150,000–180,000 years ago, with gene flow around 20,000 years ago, primarily between the Cross River and Grauer’s gorilla subspecies. This gene flow event likely obscures evolutionary relationships within eastern gorillas: after excluding putatively introgressed genomic regions, we uncover a sister relationship between Virunga mountain gorillas and Grauer’s gorillas to the exclusion of Bwindi mountain gorillas. This makes mountain gorillas paraphyletic. Eastern gorillas are less genetically diverse and more inbred than western gorillas, yet we detected lower genetic load in the eastern species. Analyses of indels fit remarkably well with differences in genetic diversity across gorilla taxa as recovered with nucleotide diversity measures. We also identified genes under selection and unique gene variants specific for each gorilla subspecies, encoding, among others, traits involved in immunity, diet, muscular development, hair morphology and behavior. The presence of this functional variation suggests that the subspecies may be locally adapted. In conclusion, using extensive genomic resources we provide a comprehensive overview of gorilla genomic diversity, including a so-far understudied Bwindi mountain gorilla population, identify putative genes involved in local adaptation, and detect population-specific gene flow across gorilla species.

Genomic and cytogenetic analyses reveal satellite repeat signature in allotetraploid okra (Abelmoschus esculentus)

BackgroundSatellite repeats are one of the most rapidly evolving components in eukaryotic genomes and play vital roles in genome regulation, genome evolution, and speciation. As a consequence, the composition, abundance and chromosome distribution of satellite repeats often exhibit variability across various species, genome, and even individual chromosomes. However, we know little about the satellite repeat evolution in allopolyploid genomes.ResultsIn this study, we investigated the satellite repeat signature in five okra (Abelmoschus esculentus) accessions using genomic and cytogenetic methods. In each of the five accessions, we identified eight satellite repeats, which exhibited a significant level of intraspecific conservation. Through fluorescence in situ hybridization (FISH) experiments, we observed that the satellite repeats generated multiple signals and exhibited variations in copy number across chromosomes. Intriguingly, we found that five satellite repeats were interspersed with centromeric retrotransposons, signifying their involvement in centromeric satellite repeat identity. We confirmed subgenome-biased amplification patterns of these satellite repeats through existing genome assemblies or dual-color FISH, indicating their distinct dynamic evolution in the allotetraploid okra subgenome. Moreover, we observed the presence of multiple chromosomes harboring the 35 S rDNA loci, alongside another chromosomal pair carrying the 5 S rDNA loci in okra using FISH assay. Remarkably, the intensity of 35 S rDNA hybridization signals varied among chromosomes, with the signals predominantly localized within regions of relatively weak DAPI staining, associated with GC-rich heterochromatin regions. Finally, we observed a similar localization pattern between 35 S rDNA and three satellite repeats with high GC content and confirmed their origin in the intergenic spacer region of the 35 S rDNA.ConclusionsOur findings uncover a unique satellite repeat signature in the allotetraploid okra, contributing to our understanding of the composition, abundance, and chromosomal distribution of satellite repeats in allopolyploid genomes, further enriching our understanding of their evolutionary dynamics in complex allopolyploid genomes.