Aligning multiple genomic sequences with the threaded blockset aligner

FishTEDB 2.0: an update fish transposable element (TE) database with new functions to facilitate TE research

We launched the initial version of FishTEDB in 2018, which aimed to establish an open-source, user-friendly, data-rich transposable element (TE) database. Over the past 5 years, FishTEDB 1.0 has gained approximately 10 000 users, accumulating more than 450 000 interactions. With the unveiling of extensive fish genome data and the increasing emphasis on TE research, FishTEDB needs to extend the richness of data and functions. To achieve the above goals, we introduced 33 new fish species to FishTEDB 2.0, encompassing a wide array of fish belonging to 48 orders. To make the updated database more functional, we added a genome browser to visualize the positional relationship between TEs and genes and the estimated TE insertion time in different species. In conclusion, we released a new version of the fish TE database, FishTEDB 2.0, designed to assist researchers in the future study of TE functions and promote the progress of biological theories related to TEs. Database URL: https://www.fishtedb.com/.

ATP1A3 Disease Spectrum Includes Paroxysmal Weakness and Encephalopathy Not Triggered by Fever

Background and Objectives

Heterozygous pathogenic variants in ATP1A3, which encodes the catalytic alpha subunit of neuronal Na+/K+-ATPase, cause primarily neurologic disorders with widely variable features that can include episodic movement deficits. One distinctive presentation of ATP1A3-related disease is recurrent fever-triggered encephalopathy. This can occur with generalized weakness and/or ataxia and is described in the literature as relapsing encephalopathy with cerebellar ataxia. This syndrome displays genotype-phenotype correlation with variants at p.R756 causing temperature sensitivity of ATP1A3. We report clinical and in vitro functional evidence for a similar phenotype not triggered by fever but associated with protein loss-of-function.

Methods

We describe the phenotype of an individual with de novo occurrence of a novel heterozygous ATP1A3 variant, NM_152296.5:c.388_390delGTG; p.(V130del). We confirmed the pathogenicity of p.V130del by cell survival complementation assay in HEK293 cells and then characterized its functional impact on enzymatic ion transport and extracellular sodium binding by two-electrode voltage clamp electrophysiology in Xenopus oocytes. To determine whether variant enzymes reach the cell surface, we surface-biotinylated oocytes expressing N-tagged ATP1A3.

Results

The proband is a 7-year-old boy who has had 2 lifetime episodes of paroxysmal weakness, encephalopathy, and ataxia not triggered by fever. He had speech regression and intermittent hand tremors after the second episode but otherwise spontaneously recovered after episodes and is at present developmentally appropriate. The p.V130del variant was identified on clinical trio exome sequencing, which did not reveal any other variants possibly associated with the phenotype. p.V130del eliminated ATP1A3 function in cell survival complementation assay. In Xenopus oocytes, p.V130del variant Na+/K+-ATPases showed complete loss of ion transport activity and marked abnormalities of extracellular Na+ binding at room temperature. Despite this clear loss-of-function effect, surface biotinylation under the same conditions revealed that p.V130del variant enzymes were still present at the oocyte's cell membrane.

Discussion

This individual's phenotype expands the clinical spectrum of ATP1A3-related recurrent encephalopathy to include presentations without fever-triggered events. The total loss of ion transport function with p.V130del, despite enzyme presence at the cell membrane, indicates that haploinsufficiency can cause relatively mild phenotypes in ATP1A3-related disease.

Genomic signatures associated with recurrent scale loss in cyprinid fish

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

Identification and functional characterization of conserved cis-regulatory elements responsible for early fruit development in cucurbit crops

A number of cis-regulatory elements (CREs) conserved during evolution have been found to be responsible for phenotypic novelty and variation. Cucurbit crops such as cucumber (Cucumis sativus), watermelon (Citrullus lanatus), melon (Cucumis melo), and squash (Cucurbita maxima) develop fruits from an inferior ovary and share some similar biological processes during fruit development. Whether conserved regulatory sequences play critical roles in fruit development of cucurbit crops remains to be explored. In six well-studied cucurbit species, we identified 392,438 conserved noncoding sequences (CNSs), including 82,756 that are specific to cucurbits, by comparative genomics. Genome-wide profiling of accessible chromatin regions (ACRs) and gene expression patterns mapped 20,865 to 43,204 ACRs and their potential target genes for two fruit tissues at two key developmental stages in six cucurbits. Integrated analysis of CNSs and ACRs revealed 4,431 syntenic orthologous CNSs, including 1,687 cucurbit-specific CNSs that overlap with ACRs that are present in all six cucurbit crops and that may regulate the expression of 757 adjacent orthologous genes. CRISPR mutations targeting two CNSs present in the 1,687 cucurbit-specific sequences resulted in substantially altered fruit shape and gene expression patterns of adjacent NAC1 (NAM, ATAF1/2, and CUC2) and EXT-like (EXTENSIN-like) genes, validating the regulatory roles of these CNSs in fruit development. These results not only provide a number of target CREs for cucurbit crop improvement, but also provide insight into the roles of CREs in plant biology and during evolution.

ACMGA: a reference-free multiple-genome alignment pipeline for plant species

BACKGROUND

The short-read whole-genome sequencing (WGS) approach has been widely applied to investigate the genomic variation in the natural populations of many plant species. With the rapid advancements in long-read sequencing and genome assembly technologies, high-quality genome sequences are available for a group of varieties for many plant species. These genome sequences are expected to help researchers comprehensively investigate any type of genomic variants that are missed by the WGS technology. However, multiple genome alignment (MGA) tools designed by the human genome research community might be unsuitable for plant genomes.

RESULTS

To fill this gap, we developed the AnchorWave-Cactus Multiple Genome Alignment (ACMGA) pipeline, which improved the alignment of repeat elements and could identify long (> 50 bp) deletions or insertions (INDELs). We conducted MGA using ACMGA and Cactus for 8 Arabidopsis (Arabidopsis thaliana) and 26 Maize (Zea mays) de novo assembled genome sequences and compared them with the previously published short-read variant calling results. MGA identified more single nucleotide variants (SNVs) and long INDELs than did previously published WGS variant callings. Additionally, ACMGA detected significantly more SNVs and long INDELs in repetitive regions and the whole genome than did Cactus. Compared with the results of Cactus, the results of ACMGA were more similar to the previously published variants called using short-read. These two MGA pipelines identified numerous multi-allelic variants that were missed by the WGS variant calling pipeline.

CONCLUSIONS

Aligning de novo assembled genome sequences could identify more SNVs and INDELs than mapping short-read. ACMGA combines the advantages of AnchorWave and Cactus and offers a practical solution for plant MGA by integrating global alignment, a 2-piece-affine-gap cost strategy, and the progressive MGA algorithm.

Evolution of diapause in the African turquoise killifish by remodeling the ancient gene regulatory landscape

Suspended animation states allow organisms to survive extreme environments. The African turquoise killifish has evolved diapause as a form of suspended development to survive a complete drought. However, the mechanisms underlying the evolution of extreme survival states are unknown. To understand diapause evolution, we performed integrative multi-omics (gene expression, chromatin accessibility, and lipidomics) in the embryos of multiple killifish species. We find that diapause evolved by a recent remodeling of regulatory elements at very ancient gene duplicates (paralogs) present in all vertebrates. CRISPR-Cas9-based perturbations identify the transcription factors REST/NRSF and FOXOs as critical for the diapause gene expression program, including genes involved in lipid metabolism. Indeed, diapause shows a distinct lipid profile, with an increase in triglycerides with very-long-chain fatty acids. Our work suggests a mechanism for the evolution of complex adaptations and offers strategies to promote long-term survival by activating suspended animation programs in other species.

CONSERVATION ASSESSMENT OF HUMAN SPLICE SITE ANNOTATION BASED ON A 470-GENOME ALIGNMENT

Despite many improvements over the years, the annotation of the human genome remains imperfect, and different annotations of the human reference genome sometimes contradict one another. The use of evolutionarily conserved sequences provides a strategy for selecting a high-confidence subset of the annotation that is more likely to be related to biological functions, and the rapidly growing number of genomes from other species increases its power. Using the latest whole genome alignment, we found that splice sites from protein-coding genes in the high-quality MANE annotation are consistently conserved across more than 400 species. We also studied splice sites from the RefSeq, GENCODE, and CHESS databases that are not present in MANE. We trained a logistic regression classifier to distinguish between the conservation exhibited by sites from MANE versus sites chosen randomly from neutrally evolving sequence. We found that splice sites classified by our model as conserved have lower SNP rates and better transcriptomic support. We then computed a subset of transcripts only using either "conserved" splice sites or ones from MANE. This subset is enriched in high-confidence transcripts of the major gene catalogs that appear to be under purifying selection and are more likely to be correct and functionally relevant.

The rise of baobab trees in Madagascar

The baobab trees (genus Adansonia) have attracted tremendous attention because of their striking shape and distinctive relationships with fauna1. These spectacular trees have also influenced human culture, inspiring innumerable arts, folklore and traditions. Here we sequenced genomes of all eight extant baobab species and argue that Madagascar should be considered the centre of origin for the extant lineages, a key issue in their evolutionary history2,3. Integrated genomic and ecological analyses revealed the reticulate evolution of baobabs, which eventually led to the species diversity seen today. Past population dynamics of Malagasy baobabs may have been influenced by both interspecific competition and the geological history of the island, especially changes in local sea levels. We propose that further attention should be paid to the conservation status of Malagasy baobabs, especially of Adansonia suarezensis and Adansonia grandidieri, and that intensive monitoring of populations of Adansonia za is required, given its propensity for negatively impacting the critically endangered Adansonia perrieri.

The NCBI Comparative Genome Viewer (CGV) is an interactive visualization tool for the analysis of whole-genome eukaryotic alignments

We report a new visualization tool for analysis of whole-genome assembly-assembly alignments, the Comparative Genome Viewer (CGV) (https://ncbi.nlm.nih.gov/genome/cgv/). CGV visualizes pairwise same-species and cross-species alignments provided by National Center for Biotechnology Information (NCBI) using assembly alignment algorithms developed by us and others. Researchers can examine large structural differences spanning chromosomes, such as inversions or translocations. Users can also navigate to regions of interest, where they can detect and analyze smaller-scale deletions and rearrangements within specific chromosome or gene regions. RefSeq or user-provided gene annotation is displayed where available. CGV currently provides approximately 800 alignments from over 350 animal, plant, and fungal species. CGV and related NCBI viewers are undergoing active development to further meet needs of the research community in comparative genome visualization.

The rate of adaptive molecular evolution in wild and domesticated Saccharomyces cerevisiae populations

Through its fermentative capacities, Saccharomyces cerevisiae was central in the development of civilisation during the Neolithic period, and the yeast remains of importance in industry and biotechnology, giving rise to bona fide domesticated populations. Here, we conduct a population genomic study of domesticated and wild populations of S. cerevisiae. Using coalescent analyses, we report that the effective population size of yeast populations decreased since the divergence with S. paradoxus. We fitted models of distributions of fitness effects to infer the rate of adaptive (ω a ) and non-adaptive (ω na ) non-synonymous substitutions in protein-coding genes. We report an overall limited contribution of positive selection to S. cerevisiae protein evolution, albeit with higher rates of adaptive evolution in wild compared to domesticated populations. Our analyses revealed the signature of background selection and possibly Hill-Robertson interference, as recombination was found to be negatively correlated withω na and positively correlated withω a . However, the effect of recombination onω a was found to be labile, as it is only apparent after removing the impact of codon usage bias on the synonymous site frequency spectrum and disappears if we control for the correlation withω na , suggesting that it could be an artefact of the decreasing population size. Furthermore, the rate of adaptive non-synonymous substitutions is significantly correlated with the residue solvent exposure, a relation that cannot be explained by the population's demography. Together, our results provide a detailed characterisation of adaptive mutations in protein-coding genes across S. cerevisiae populations.

A mammalian tripartite enhancer cluster controls hypothalamic Pomc expression, food intake, and body weight

Food intake and energy balance are tightly regulated by a group of hypothalamic arcuate neurons expressing the proopiomelanocortin (POMC) gene. In mammals, arcuate-specific POMC expression is driven by two cis-acting transcriptional enhancers known as nPE1 and nPE2. Because mutant mice lacking these two enhancers still showed hypothalamic Pomc mRNA, we searched for additional elements contributing to arcuate Pomc expression. By combining molecular evolution with reporter gene expression in transgenic zebrafish and mice, here, we identified a mammalian arcuate-specific Pomc enhancer that we named nPE3, carrying several binding sites also present in nPE1 and nPE2 for transcription factors known to activate neuronal Pomc expression, such as ISL1, NKX2.1, and ERα. We found that nPE3 originated in the lineage leading to placental mammals and remained under purifying selection in all mammalian orders, although it was lost in Simiiformes (monkeys, apes, and humans) following a unique segmental deletion event. Interestingly, ablation of nPE3 from the mouse genome led to a drastic reduction (>70%) in hypothalamic Pomc mRNA during development and only moderate (<33%) in adult mice. Comparison between double (nPE1 and nPE2) and triple (nPE1, nPE2, and nPE3) enhancer mutants revealed the relative contribution of nPE3 to hypothalamic Pomc expression and its importance in the control of food intake and adiposity in male and female mice. Altogether, these results demonstrate that nPE3 integrates a tripartite cluster of partially redundant enhancers that originated upon a triple convergent evolutionary process in mammals and that is critical for hypothalamic Pomc expression and body weight homeostasis.

Genomic insights into the seawater adaptation in Cyprinidae

BACKGROUND

Cyprinidae, the largest fish family, encompasses approximately 367 genera and 3006 species. While they exhibit remarkable adaptability to diverse aquatic environments, it is exceptionally rare to find them in seawater, with the Far Eastern daces being of few exceptions. Therefore, the Far Eastern daces serve as a valuable model for studying the genetic mechanisms underlying seawater adaptation in Cyprinidae.

RESULTS

Here, we sequenced the chromosome-level genomes of two Far Eastern daces (Pseudaspius brandtii and P. hakonensis), the two known cyprinid fishes found in seawater, and performed comparative genomic analyses to investigate their genetic mechanism of seawater adaptation. Demographic history reconstruction of the two species reveals that their population dynamics are correlated with the glacial-interglacial cycles and sea level changes. Genomic analyses identified Pseudaspius-specific genetic innovations related to seawater adaptation, including positively selected genes, rapidly evolving genes, and conserved non-coding elements (CNEs). Functional assays of Pseudaspius-specific variants of the prolactin (prl) gene showed enhanced cell adaptation to greater osmolarity. Functional assays of Pseudaspius specific CNEs near atg7 and usp45 genes suggest that they exhibit higher promoter activity and significantly induced at high osmolarity.

CONCLUSIONS

Our results reveal the genome-wide evidence for the evolutionary adaptation of cyprinid fishes to seawater, offering valuable insights into the molecular mechanisms supporting the survival of migratory fish in marine environments. These findings are significant as they contribute to our understanding of how cyprinid fishes navigate and thrive in diverse aquatic habitats, providing useful implications for the conservation and management of marine ecosystems.

The rapidly evolving X-linked MIR-506 family fine-tunes spermatogenesis to enhance sperm competition

Despite rapid evolution across eutherian mammals, the X-linked MIR-506 family miRNAs are located in a region flanked by two highly conserved protein-coding genes (SLITRK2 and FMR1) on the X chromosome. Intriguingly, these miRNAs are predominantly expressed in the testis, suggesting a potential role in spermatogenesis and male fertility. Here, we report that the X-linked MIR-506 family miRNAs were derived from the MER91C DNA transposons. Selective inactivation of individual miRNAs or clusters caused no discernible defects, but simultaneous ablation of five clusters containing 19 members of the MIR-506 family led to reduced male fertility in mice. Despite normal sperm counts, motility, and morphology, the KO sperm were less competitive than wild-type sperm when subjected to a polyandrous mating scheme. Transcriptomic and bioinformatic analyses revealed that these X-linked MIR-506 family miRNAs, in addition to targeting a set of conserved genes, have more targets that are critical for spermatogenesis and embryonic development during evolution. Our data suggest that the MIR-506 family miRNAs function to enhance sperm competitiveness and reproductive fitness of the male by finetuning gene expression during spermatogenesis.

GPN-MSA: an alignment-based DNA language model for genome-wide variant effect prediction

Whereas protein language models have demonstrated remarkable efficacy in predicting the effects of missense variants, DNA counterparts have not yet achieved a similar competitive edge for genome-wide variant effect predictions, especially in complex genomes such as that of humans. To address this challenge, we here introduce GPN-MSA, a novel framework for DNA language models that leverages whole-genome sequence alignments across multiple species and takes only a few hours to train. Across several benchmarks on clinical databases (ClinVar, COSMIC, OMIM), experimental functional assays (DMS, DepMap), and population genomic data (gnomAD), our model for the human genome achieves outstanding performance on deleteriousness prediction for both coding and non-coding variants.

Pathogenic variants in human DNA damage repair genes mostly arose in recent human history

BACKGROUND

Genome stability is maintained by the DNA damage repair (DDR) system composed of multiple DNA repair pathways of hundreds of genes. Germline pathogenic variation (PV) in DDR genes damages function of the affected DDR genes, leading to genome instability and high risk of diseases, in particular, cancer. Knowing evolutionary origin of the PVs in human DDR genes is essential to understand the etiology of human diseases. However, answer to the issue remains largely elusive. In this study, we analyzed evolutionary origin for the PVs in human DDR genes.

METHODS

We identified 169 DDR genes by referring to various databases and identified PVs in the DDR genes of modern humans from ClinVar database. We performed a phylogenetic analysis to analyze the conservation of human DDR PVs in 100 vertebrates through cross-species genomic data comparison using the phyloFit program of the PHAST package and visualized the results using the GraphPad Prism software and the ggplot module. We identified DDR PVs from over 5000 ancient humans developed a database to host the DDR PVs ( https://genemutation.fhs.um.edu.mo/dbDDR-AncientHumans ). Using the PV data, we performed a molecular archeological analysis to compare the DDR PVs between modern humans and ancient humans. We analyzed evolution selection of DDR genes across 20 vertebrates using the CodeML in PAML for phylogenetic analysis.

RESULTS

Our phylogenic analysis ruled out cross-species conservation as the origin of human DDR PVs. Our archeological approach identified rich DDR PVs shared between modern and ancient humans, which were mostly dated within the last 5000 years. We also observed similar pattern of quantitative PV distribution between modern and ancient humans. We further detected a set of ATM, BRCA2 and CHEK2 PVs shared between human and Neanderthals.

CONCLUSIONS

Our study reveals that human DDR PVs mostly arose in recent human history. We propose that human high cancer risk caused by DDR PVs can be a by-product of human evolution.

Unveiling hidden structural patterns in the SARS-CoV-2 genome: Computational insights and comparative analysis

SARS-CoV-2, the causative agent of COVID-19, is known to exhibit secondary structures in its 5' and 3' untranslated regions, along with the frameshifting stimulatory element situated between ORF1a and 1b. To identify additional regions containing conserved structures, we utilized a multiple sequence alignment with related coronaviruses as a starting point. We applied a computational pipeline developed for identifying non-coding RNA elements. Our pipeline employed three different RNA structural prediction approaches. We identified forty genomic regions likely to harbor structures, with ten of them showing three-way consensus substructure predictions among our predictive utilities. We conducted intracomparisons of the predictive utilities within the pipeline and intercomparisons with four previously published SARS-CoV-2 structural datasets. While there was limited agreement on the precise structure, different approaches seemed to converge on regions likely to contain structures in the viral genome. By comparing and combining various computational approaches, we can predict regions most likely to form structures, as well as a probable structure or ensemble of structures. These predictions can be used to guide surveillance, prophylactic measures, or therapeutic efforts. Data and scripts employed in this study may be found at https://doi.org/10.5281/zenodo.8298680.

miR-317-3p and miR-283-5p Play a Crucial Role in Regulating the Resistance to Indoxacarb in Spodoptera frugiperda by Targeting GSTs4

Spodoptera frugiperda is primarily controlled through chemical insecticides. Our RNA-seq data highlight the overexpression of GSTs4 in indoxacarb-resistant S. frugiperda. However, the exact role of GSTs4 in indoxacarb resistance and its regulatory mechanisms remains elusive. Therefore, we investigated the functional role of GSTs4 in S. frugiperda and explored the underlying post-transcriptional regulatory mechanisms. GSTs4 was highly overexpressed (27.6-fold) in the indoxacarb-resistant strain, and GSTs4 silencing significantly increases the susceptibility of S. frugiperda to indoxacarb, increasing mortality by 27.3%. miR-317-3p and miR-283-5p can bind to the 3'UTR of GSTs4, and the targeting relationship was confirmed by dual-luciferase reporter assays. Injecting miR-317-3p and miR-283-5p agomirs reduces GSTs4 levels by 64.8 and 42.3%, respectively, resulting in an increased susceptibility of S. frugiperda to indoxacarb. Conversely, the administration of miR-317-3p and miR-283-5pantagomirs increases GSTs4 expression and reduces larval susceptibility to indoxacarb. These findings demonstrate that miR-317-3p and miR-283-5p contribute to indoxacarb resistance in S. frugiperda by regulating the overexpression of GSTs4.

Pangenome graph construction from genome alignments with Minigraph-Cactus

Pangenome references address biases of reference genomes by storing a representative set of diverse haplotypes and their alignment, usually as a graph. Alternate alleles determined by variant callers can be used to construct pangenome graphs, but advances in long-read sequencing are leading to widely available, high-quality phased assemblies. Constructing a pangenome graph directly from assemblies, as opposed to variant calls, leverages the graph's ability to represent variation at different scales. Here we present the Minigraph-Cactus pangenome pipeline, which creates pangenomes directly from whole-genome alignments, and demonstrate its ability to scale to 90 human haplotypes from the Human Pangenome Reference Consortium. The method builds graphs containing all forms of genetic variation while allowing use of current mapping and genotyping tools. We measure the effect of the quality and completeness of reference genomes used for analysis within the pangenomes and show that using the CHM13 reference from the Telomere-to-Telomere Consortium improves the accuracy of our methods. We also demonstrate construction of a Drosophila melanogaster pangenome.

A phylogenetic method linking nucleotide substitution rates to rates of continuous trait evolution

Genomes contain conserved non-coding sequences that perform important biological functions, such as gene regulation. We present a phylogenetic method, PhyloAcc-C, that associates nucleotide substitution rates with changes in a continuous trait of interest. The method takes as input a multiple sequence alignment of conserved elements, continuous trait data observed in extant species, and a background phylogeny and substitution process. Gibbs sampling is used to assign rate categories (background, conserved, accelerated) to lineages and explore whether the assigned rate categories are associated with increases or decreases in the rate of trait evolution. We test our method using simulations and then illustrate its application using mammalian body size and lifespan data previously analyzed with respect to protein coding genes. Like other studies, we find processes such as tumor suppression, telomere maintenance, and p53 regulation to be related to changes in longevity and body size. In addition, we also find that skeletal genes, and developmental processes, such as sprouting angiogenesis, are relevant.

The growth factor EPIREGULIN promotes basal progenitor cell proliferation in the developing neocortex

Neocortex expansion during evolution is linked to higher numbers of neurons, which are thought to result from increased proliferative capacity and neurogenic potential of basal progenitor cells during development. Here, we show that EREG, encoding the growth factor EPIREGULIN, is expressed in the human developing neocortex and in gorilla cerebral organoids, but not in the mouse neocortex. Addition of EPIREGULIN to the mouse neocortex increases proliferation of basal progenitor cells, whereas EREG ablation in human cortical organoids reduces proliferation in the subventricular zone. Treatment of cortical organoids with EPIREGULIN promotes a further increase in proliferation of gorilla but not of human basal progenitor cells. EPIREGULIN competes with the epidermal growth factor (EGF) to promote proliferation, and inhibition of the EGF receptor abrogates the EPIREGULIN-mediated increase in basal progenitor cells. Finally, we identify putative cis-regulatory elements that may contribute to the observed inter-species differences in EREG expression. Our findings suggest that species-specific regulation of EPIREGULIN expression may contribute to the increased neocortex size of primates by providing a tunable pro-proliferative signal to basal progenitor cells in the subventricular zone.

BRD2 and BRD3 genes independently evolved RNA structures to control unproductive splicing

The mammalian BRD2 and BRD3 genes encode structurally related proteins from the bromodomain and extraterminal domain protein family. The expression of BRD2 is regulated by unproductive splicing upon inclusion of exon 3b, which is located in the region encoding a bromodomain. Bioinformatic analysis indicated that BRD2 exon 3b inclusion is controlled by a pair of conserved complementary regions (PCCR) located in the flanking introns. Furthermore, we identified a highly conserved element encoding a cryptic poison exon 5b and a previously unknown PCCR in the intron between exons 5 and 6 of BRD3, however, outside of the homologous bromodomain. Minigene mutagenesis and blockage of RNA structure by antisense oligonucleotides demonstrated that RNA structure controls the rate of inclusion of poison exons. The patterns of BRD2 and BRD3 expression and splicing show downregulation upon inclusion of poison exons, which become skipped in response to transcription elongation slowdown, further confirming a role of PCCRs in unproductive splicing regulation. We conclude that BRD2 and BRD3 independently acquired poison exons and RNA structures to dynamically control unproductive splicing. This study describes a convergent evolution of regulatory unproductive splicing mechanisms in these genes, providing implications for selective modulation of their expression in therapeutic applications.

Detection of evolutionary conserved and accelerated genomic regions related to adaptation to thermal niches in Anolis lizards

Understanding the genetic basis for adapting to thermal environments is important due to serious effects of global warming on ectothermic species. Various genes associated with thermal adaptation in lizards have been identified mainly focusing on changes in gene expression or the detection of positively selected genes using coding regions. Only a few comprehensive genome-wide analyses have included noncoding regions. This study aimed to identify evolutionarily conserved and accelerated genomic regions using whole genomes of eight Anolis lizard species that have repeatedly adapted to similar thermal environments in multiple lineages. Evolutionarily conserved genomic regions were extracted as regions with overall sequence conservation (regions with fewer base substitutions) across all lineages compared with the neutral model. Genomic regions that underwent accelerated evolution in the lineage of interest were identified as those with more base substitutions in the target branch than in the entire background branch. Conserved elements across all branches were relatively abundant in "intergenic" genomic regions among noncoding regions. Accelerated regions (ARs) of each lineage contained a significantly greater proportion of noncoding RNA genes than the entire multiple alignment. Common genes containing ARs within 5 kb of their vicinity in lineages with similar thermal habitats were identified. Many genes associated with circadian rhythms and behavior were found in hot-open and cool-shaded habitat lineages. These genes might play a role in contributing to thermal adaptation and assist future studies examining the function of genes involved in thermal adaptation via genome editing.

Identification of a distal enhancer regulating hedgehog interacting protein gene in human lung epithelial cells

BACKGROUND

An intergenic region at chromosome 4q31 is one of the most significant regions associated with COPD susceptibility and lung function in GWAS. In this region, the implicated causal gene HHIP has a unique epithelial expression pattern in adult human lungs, in contrast to dominant expression in fibroblasts in murine lungs. However, the mechanism underlying the species-dependent cell type-specific regulation of HHIP remains largely unknown.

METHODS

We employed snATAC-seq analysis to identify open chromatin regions within the COPD GWAS region in various human lung cell types. ChIP-quantitative PCR, reporter assays, chromatin conformation capture assays and Hi-C assays were conducted to characterize the regulatory element in this region. CRISPR/Cas9-editing was performed in BEAS-2B cells to generate single colonies with stable knockout of the regulatory element. RT-PCR and Western blot assays were used to evaluate expression of HHIP and epithelial-mesenchymal transition (EMT)-related marker genes.

FINDINGS

We identified a distal enhancer within the COPD 4q31 GWAS locus that regulates HHIP transcription at baseline and after TGFβ treatment in a SMAD3-dependent, but Hedgehog-independent manner in human bronchial epithelial cells. The distal enhancer also maintains chromatin topological domains near 4q31 locus and HHIP gene. Reduced HHIP expression led to increased EMT induced by TGFβ in human bronchial epithelial cells.

INTERPRETATION

A distal enhancer regulates HHIP expression both under homeostatic condition and upon TGFβ treatment in human bronchial epithelial cells. The interaction between HHIP and TGFβ signalling possibly contributes to COPD pathogenesis.

FUNDING

Supported by NIH grants R01HL127200, R01HL148667 and R01HL162783 (to X. Z).

New whole-genome alignment tools are needed for tapping into plant diversity

Genome alignment is one of the most foundational methods for genome sequence studies. With rapid advances in sequencing and assembly technologies, these newly assembled genomes present challenges for alignment tools to meet the increased complexity and scale. Plant genome alignment is technologically challenging because of frequent whole-genome duplications (WGDs) as well as chromosome rearrangements and fractionation, high nucleotide diversity, widespread structural variation, and high transposable element (TE) activity causing large proportions of repeat elements. We summarize classical pairwise and multiple genome alignment (MGA) methods, and highlight techniques that are widely used or are being developed by the plant research community. We also outline the remaining challenges for precise genome alignment and the interpretation of alignment results in plants.

A quantitative genetic model of background selection in humans

Across the human genome, there are large-scale fluctuations in genetic diversity caused by the indirect effects of selection. This "linked selection signal" reflects the impact of selection according to the physical placement of functional regions and recombination rates along chromosomes. Previous work has shown that purifying selection acting against the steady influx of new deleterious mutations at functional portions of the genome shapes patterns of genomic variation. To date, statistical efforts to estimate purifying selection parameters from linked selection models have relied on classic Background Selection theory, which is only applicable when new mutations are so deleterious that they cannot fix in the population. Here, we develop a statistical method based on a quantitative genetics view of linked selection, that models how polygenic additive fitness variance distributed along the genome increases the rate of stochastic allele frequency change. By jointly predicting the equilibrium fitness variance and substitution rate due to both strong and weakly deleterious mutations, we estimate the distribution of fitness effects (DFE) and mutation rate across three geographically distinct human samples. While our model can accommodate weaker selection, we find evidence of strong selection operating similarly across all human samples. Although our quantitative genetic model of linked selection fits better than previous models, substitution rates of the most constrained sites disagree with observed divergence levels. We find that a model incorporating selective interference better predicts observed divergence in conserved regions, but overall our results suggest uncertainty remains about the processes generating fitness variation in humans.

Discovery of novel microRNA mimic repressors of ribosome biogenesis

While microRNAs and other non-coding RNAs are the next frontier of novel regulators of mammalian ribosome biogenesis (RB), a systematic exploration of microRNA-mediated RB regulation has not yet been undertaken. We carried out a high-content screen in MCF10A cells for changes in nucleolar number using a library of 2603 mature human microRNA mimics. Following a secondary screen for nucleolar rRNA biogenesis inhibition, we identified 72 novel microRNA negative regulators of RB after stringent hit calling. Hits included 27 well-conserved microRNAs present in MirGeneDB, and were enriched for mRNA targets encoding proteins with nucleolar localization or functions in cell cycle regulation. Rigorous selection and validation of a subset of 15 microRNA hits unexpectedly revealed that most of them caused dysregulated pre-rRNA processing, elucidating a novel role for microRNAs in RB regulation. Almost all hits impaired global protein synthesis and upregulated CDKN1A (p21) levels, while causing diverse effects on RNA Polymerase 1 (RNAP1) transcription and TP53 protein levels. We provide evidence that the MIR-28 siblings, hsa-miR-28-5p and hsa-miR-708-5p, potently target the ribosomal protein mRNA RPS28 via tandem primate-specific 3' UTR binding sites, causing a severe pre-18S pre-rRNA processing defect. Our work illuminates novel microRNA attenuators of RB, forging a promising new path for microRNA mimic chemotherapeutics.

Genomic evidence for climate-linked diversity loss and increased vulnerability of wild barley spanning 28 years of climate warming

The information on how plant populations respond genetically to climate warming is scarce. Here, landscape genomic and machine learning approaches were integrated to assess genetic response of 10 wild barley (Hordeum vulgare ssp. spontaneum; WB) populations in the past and future, using whole genomic sequencing (WGS) data. The WB populations were sampled in 1980 and again in 2008. Phylogeny of accessions was roughly in conformity with sampling sites, which accompanied by admixture/introgressions. The 28-y climate warming resulted in decreased genetic diversity, increased selection pressure, and an increase in deleterious single nucleotide polymorphism (dSNP) numbers, heterozygous deleterious and total deleterious burdens for WB. Genome-environment associations identified some candidate genes belonging to peroxidase family (HORVU2Hr1G057450, HORVU4Hr1G052060 and HORVU4Hr1G057210) and heat shock protein 70 family (HORVU2Hr1G112630). The gene HORVU2Hr1G120170 identified by selective sweep analysis was under strong selection during the climate warming of the 28-y, and its derived haplotypes were fixed by WB when faced with the 28-y increasingly severe environment. Temperature variables were found to be more important than precipitation variables in influencing genomic variation, with an eco-physiological index gdd5 (growing degree-days at the baseline threshold temperature of 5 °C) being the most important determinant. Gradient forest modelling revealed higher predicted genomic vulnerability in Sede Boqer under future climate scenarios at 2041-2070 and 2071-2100. Additionally, estimates of effective population size (Ne) tracing back to 250 years indicated a forward decline in all populations over time. Our assessment about past genetic response and future vulnerability of WB under climate warming is crucial for informing conservation efforts for wild cereals and rational use strategies.

Graphylo: A deep learning approach for predicting regulatory DNA and RNA sites from whole-genome multiple alignments

This study focuses on enhancing the prediction of regulatory functional sites in DNA and RNA sequences, a crucial aspect of gene regulation. Current methods, such as motif overrepresentation and machine learning, often lack specificity. To address this issue, the study leverages evolutionary information and introduces Graphylo, a deep-learning approach for predicting transcription factor binding sites in the human genome. Graphylo combines Convolutional Neural Networks for DNA sequences with Graph Convolutional Networks on phylogenetic trees, using information from placental mammals' genomes and evolutionary history. The research demonstrates that Graphylo consistently outperforms both single-species deep learning techniques and methods that incorporate inter-species conservation scores on a wide range of datasets. It achieves this by utilizing a species-based attention model for evolutionary insights and an integrated gradient approach for nucleotide-level model interpretability. This innovative approach offers a promising avenue for improving the accuracy of regulatory site prediction in genomics.

Dissecting the sequence and structural determinants guiding m6A deposition and evolution via inter- and intra-species hybrids

BACKGROUND

N6-methyladenosine (m6A) is the most abundant mRNA modification, and controls mRNA stability. m6A distribution varies considerably between and within species. Yet, it is unclear to what extent this variability is driven by changes in genetic sequences ('cis') or cellular environments ('trans') and via which mechanisms.

RESULTS

Here we dissect the determinants governing RNA methylation via interspecies and intraspecies hybrids in yeast and mammalian systems, coupled with massively parallel reporter assays and m6A-QTL reanalysis. We find that m6A evolution and variability is driven primarily in 'cis', via two mechanisms: (1) variations altering m6A consensus motifs, and (2) variation impacting mRNA secondary structure. We establish that mutations impacting RNA structure - even when distant from an m6A consensus motif - causally dictate methylation propensity. Finally, we demonstrate that allele-specific differences in m6A levels lead to allele-specific changes in gene expression.

CONCLUSIONS

Our findings define the determinants governing m6A evolution and diversity and characterize the consequences thereof on gene expression regulation.

Emergence of saliva protein genes in the secretory calcium-binding phosphoprotein (SCPP) locus and accelerated evolution in primates

Genes within the secretory calcium-binding phosphoprotein (SCPP) family evolved in conjunction with major evolutionary milestones: the formation of a calcified skeleton in vertebrates, the emergence of tooth enamel in fish, and the introduction of lactation in mammals. The SCPP gene family also contains genes expressed primarily and abundantly in human saliva. Here, we explored the evolution of the saliva-related SCPP genes by harnessing currently available genomic and transcriptomic resources. Our findings provide insights into the expansion and diversification of SCPP genes, notably identifying previously undocumented convergent gene duplications. In primate genomes, we found additional duplication and diversification events that affected genes coding for proteins secreted in saliva. These saliva-related SCPP genes exhibit signatures of positive selection in the primate lineage while the other genes in the same locus remain conserved. We found that regulatory shifts and gene turnover events facilitated the accelerated gain of salivary expression. Collectively, our results position the SCPP gene family as a hotbed of evolutionary innovation, suggesting the potential role of dietary and pathogenic pressures in the adaptive diversification of the saliva composition in primates, including humans.

Widespread incomplete lineage sorting and introgression shaped adaptive radiation in the Gossypium genus

Cotton (Gossypium) stands as a crucial economic crop, serving as the primary source of natural fiber for the textile sector. However, the evolutionary mechanisms driving speciation within the Gossypium genus remain unresolved. In this investigation, we leveraged 25 Gossypium genomes and introduced four novel assemblies-G. harknessii, G. gossypioides, G. trilobum, and G. klotzschianum (Gklo)-to delve into the speciation history of this genus. Notably, we encountered intricate phylogenies potentially stemming from introgression. These complexities are further compounded by incomplete lineage sorting (ILS), a factor likely to have been instrumental in shaping the swift diversification of cotton. Our focus subsequently shifted to the rapid radiation episode during a concise period in Gossypium evolution. For a recently diverged lineage comprising G. davidsonii, Gklo, and G. raimondii, we constructed a finely detailed ILS map. Intriguingly, this analysis revealed the non-random distribution of ILS regions across the reference Gklo genome. Moreover, we identified signs of robust natural selection influencing specific ILS regions. Noteworthy variations pertaining to speciation emerged between the closely related sister species Gklo and G. davidsonii. Approximately 15.74% of speciation structural variation genes and 12.04% of speciation-associated genes were estimated to intersect with ILS signatures. These findings enrich our understanding of the role of ILS in adaptive radiation, shedding fresh light on the intricate speciation history of the Gossypium genus.

N-glycosylation as a eukaryotic protective mechanism against protein aggregation

The tendency for proteins to form aggregates is an inherent part of every proteome and arises from the self-assembly of short protein segments called aggregation-prone regions (APRs). While posttranslational modifications (PTMs) have been implicated in modulating protein aggregation, their direct role in APRs remains poorly understood. In this study, we used a combination of proteome-wide computational analyses and biophysical techniques to investigate the potential involvement of PTMs in aggregation regulation. Our findings reveal that while most PTM types are disfavored near APRs, N-glycosylation is enriched and evolutionarily selected, especially in proteins prone to misfolding. Experimentally, we show that N-glycosylation inhibits the aggregation of peptides in vitro through steric hindrance. Moreover, mining existing proteomics data, we find that the loss of N-glycans at the flanks of APRs leads to specific protein aggregation in Neuro2a cells. Our findings indicate that, among its many molecular functions, N-glycosylation directly prevents protein aggregation in higher eukaryotes.

Cross-Species Comparative DNA Methylation Reveals Novel Insights into Complex Trait Genetics among Cattle, Sheep, and Goats

The cross-species characterization of evolutionary changes in the functional genome can facilitate the translation of genetic findings across species and the interpretation of the evolutionary basis underlying complex phenotypes. Yet, this has not been fully explored between cattle, sheep, goats, and other mammals. Here, we systematically characterized the evolutionary dynamics of DNA methylation and gene expression in 3 somatic tissues (i.e. brain, liver, and skeletal muscle) and sperm across 7 mammalian species, including 3 ruminant livestock species (cattle, sheep, and goats), humans, pigs, mice, and dogs, by generating and integrating 160 DNA methylation and transcriptomic data sets. We demonstrate dynamic changes of DNA hypomethylated regions and hypermethylated regions in tissue-type manner across cattle, sheep, and goats. Specifically, based on the phylo-epigenetic model of DNA methylome, we identified a total of 25,074 hypomethylated region extension events specific to cattle, which participated in rewiring tissue-specific regulatory network. Furthermore, by integrating genome-wide association studies of 50 cattle traits, we provided novel insights into the genetic and evolutionary basis of complex phenotypes in cattle. Overall, our study provides a valuable resource for exploring the evolutionary dynamics of the functional genome and highlights the importance of cross-species characterization of multiomics data sets for the evolutionary interpretation of complex phenotypes in cattle livestock.

The Density of Regulatory Information Is a Major Determinant of Evolutionary Constraint on Noncoding DNA in Drosophila

Evolutionary analyses have estimated that ∼60% of nucleotides in intergenic regions of the Drosophila melanogaster genome are functionally relevant, suggesting that regulatory information may be encoded more densely in intergenic regions than has been revealed by most functional dissections of regulatory DNA. Here, we approached this issue through a functional dissection of the regulatory region of the gene shavenbaby (svb). Most of the ∼90 kb of this large regulatory region is highly conserved in the genus Drosophila, though characterized enhancers occupy a small fraction of this region. By analyzing the regulation of svb in different contexts of Drosophila development, we found that the regulatory information that drives svb expression in the abdominal pupal epidermis is organized in a different way than the elements that drive svb expression in the embryonic epidermis. While in the embryonic epidermis svb is activated by compact enhancers separated by large inactive DNA regions, svb expression in the pupal epidermis is driven by regulatory information distributed over broader regions of svb cis-regulatory DNA. In the same vein, we observed that other developmental genes also display a dense distribution of putative regulatory elements in their regulatory regions. Furthermore, we found that a large percentage of conserved noncoding DNA of the Drosophila genome is contained within regions of open chromatin. These results suggest that part of the evolutionary constraint on noncoding DNA of Drosophila is explained by the density of regulatory information, which may be greater than previously appreciated.

Identification of novel microRNAs in the embryonic mouse brain using deep sequencing

Many advances in small RNA-seq technology and bioinformatics pipelines have been made recently, permitting the discovery of novel miRNAs in the embryonic day 15.5 (E15.5) mouse brain. We aimed to improve miRNA discovery in this tissue to expand our knowledge of the regulatory networks that underpin normal neurodevelopment, find new candidates for neurodevelopmental disorder aetiology, and deepen our understanding of non-coding RNA evolution. A high-quality small RNA-seq dataset of 458 M reads was generated. An unbiased miRNA discovery pipeline identified fifty putative novel miRNAs, six of which were selected for further validation. A combination of conservation analysis and target functional prediction was used to determine the authenticity of novel miRNA candidates. These findings demonstrate that miRNAs remain to be discovered, particularly if they have the features of other small RNA species.

Sequence Analysis of Six Candidate Genes in Miniature Schnauzers with Primary Hypertriglyceridemia

Miniature Schnauzers are predisposed to primary hypertriglyceridemia (HTG). In this study, we performed whole genome sequencing (WGS) of eight Miniature Schnauzers with primary HTG and screened for risk variants in six HTG candidate genes: LPL, APOC2, APOA5, GPIHBP1, LMF1, and APOE. Variants were filtered to identify those present in ≥2 Miniature Schnauzers with primary HTG and uncommon (<10% allele frequency) in a WGS variant database including 613 dogs from 61 other breeds. Three variants passed filtering: an APOE TATA box deletion, an LMF1 intronic SNP, and a GPIHBP1 missense variant. The APOE and GPIHBP1 variants were genotyped in a cohort of 108 Miniature Schnauzers, including 68 with primary HTG and 40 controls. A multivariable regression model, including age and sex, did not identify an effect of APOE (estimate = 0.18, std. error = 0.14; p = 0.20) or GPIHBP1 genotypes (estimate = -0.26, std. error = 0.42; p = 0.54) on triglyceride concentration. In conclusion, we did not identify a monogenic cause for primary HTG in Miniature Schnauzers in the six genes evaluated. However, if HTG in Miniature Schnauzers is a complex disease resulting from the cumulative effects of multiple variants and environment, the identified variants cannot be ruled out as contributing factors.

Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

The genome of the black-footed cat: Revealing a rich natural history and urgent conservation priorities for small felids

Habitat degradation and loss of genetic diversity are common threats faced by almost all of today's wild cats. Big cats, such as tigers and lions, are of great concern and have received considerable conservation attention through policies and international actions. However, knowledge of and conservation actions for small wild cats are lagging considerably behind. The black-footed cat, Felis nigripes, one of the smallest felid species, is experiencing increasing threats with a rapid reduction in population size. However, there is a lack of genetic information to assist in developing effective conservation actions. A de novo assembly of a high-quality chromosome-level reference genome of the black-footed cat was made, and comparative genomics and population genomics analyses were carried out. These analyses revealed that the most significant genetic changes in the evolution of the black-footed cat are the rapid evolution of sensory and metabolic-related genes, reflecting genetic adaptations to its characteristic nocturnal hunting and a high metabolic rate. Genomes of the black-footed cat exhibit a high level of inbreeding, especially for signals of recent inbreeding events, which suggest that they may have experienced severe genetic isolation caused by habitat fragmentation. More importantly, inbreeding associated with two deleterious mutated genes may exacerbate the risk of amyloidosis, the dominant disease that causes mortality of about 70% of captive individuals. Our research provides comprehensive documentation of the evolutionary history of the black-footed cat and suggests that there is an urgent need to investigate genomic variations of small felids worldwide to support effective conservation actions.

Phylogenomics analyses of all species of Swordtails (Genus Xiphophorus ) highlights hybridization precedes speciation

Hybridization has been recognized as an important driving force for evolution, however studies of the genetic consequence and its cause are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the central American genus Xiphophorus were proposed to have evolved with multiple ancient and ongoing hybridization events, and served as a valuable research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genome resource and its annotation of all 26 Xiphophorus species. On this dataset we resolved the so far conflicting phylogeny. Through comparative genomic analyses we investigated the molecular evolution of genes related to melanoma, for a main sexually selected trait and for the genetic control of puberty timing, which are predicted to be involved in pre-and postzygotic isolation and thus to influence the probability of interspecific hybridization in Xiphophorus . We demonstrate dramatic size-variation of some gene families across species, despite the reticulate evolution and short divergence time. Finally, we clarify the hybridization history in the genus Xiphophorus genus, settle the long dispute on the hybridization origin of two Southern swordtails, highlight hybridizations precedes speciation, and reveal the distribution of hybridization ancestry remaining in the fused genome.

Chromosome-level genome assembly of hadal snailfish reveals mechanisms of deep-sea adaptation in vertebrates

As the deepest vertebrate in the ocean, the hadal snailfish (Pseudoliparis swirei), which lives at a depth of 6,000-8,000 m, is a representative case for studying adaptation to extreme environments. Despite some preliminary studies on this species in recent years, including their loss of pigmentation, visual and skeletal calcification genes, and the role of trimethylamine N-oxide in adaptation to high-hydrostatic pressure, it is still unknown how they evolved and why they are among the few vertebrate species that have successfully adapted to the deep-sea environment. Using genomic data from different trenches, we found that the hadal snailfish may have entered and fully adapted to such extreme environments only in the last few million years. Meanwhile, phylogenetic relationships show that they spread into different trenches in the Pacific Ocean within a million years. Comparative genomic analysis has also revealed that the genes associated with perception, circadian rhythms, and metabolism have been extensively modified in the hadal snailfish to adapt to its unique environment. More importantly, the tandem duplication of a gene encoding ferritin significantly increased their tolerance to reactive oxygen species, which may be one of the important factors in their adaptation to high-hydrostatic pressure.

Pleiotropic influence of DNA methylation QTLs on physiological and ageing traits

DNA methylation is influenced by genetic and non-genetic factors. Here, we chart quantitative trait loci (QTLs) that modulate levels of methylation at highly conserved CpGs using liver methylome data from mouse strains belonging to the BXD family. A regulatory hotspot on chromosome 5 had the highest density of trans-acting methylation QTLs (trans-meQTLs) associated with multiple distant CpGs. We refer to this locus as meQTL.5a. Trans-modulated CpGs showed age-dependent changes and were enriched in developmental genes, including several members of the MODY pathway (maturity onset diabetes of the young). The joint modulation by genotype and ageing resulted in a more 'aged methylome' for BXD strains that inherited the DBA/2J parental allele at meQTL.5a. Further, several gene expression traits, body weight, and lipid levels mapped to meQTL.5a, and there was a modest linkage with lifespan. DNA binding motif and protein-protein interaction enrichment analyses identified the hepatic nuclear factor, Hnf1a (MODY3 gene in humans), as a strong candidate. The pleiotropic effects of meQTL.5a could contribute to variations in body size and metabolic traits, and influence CpG methylation and epigenetic ageing that could have an impact on lifespan.

An intronic LINE-1 regulates IFNAR1 expression in human immune cells

BACKGROUND

Despite their origins as selfish parasitic sequences, some transposons in the human genome have been co-opted to serve as regulatory elements, contributing to the evolution of transcriptional networks. Most well-characterized examples of transposon-derived regulatory elements derive from endogenous retroviruses (ERVs), due to the intrinsic regulatory activity of proviral long terminal repeat regions. However, one subclass of transposable elements, the Long Interspersed Nuclear Elements (LINEs), have been largely overlooked in the search for functional regulatory transposons, and considered to be broadly epigenetically repressed.

RESULTS

We examined the chromatin state of LINEs by analyzing epigenomic data from human immune cells. Many LINEs are marked by the repressive H3K9me3 modification, but a subset exhibits evidence of enhancer activity in human immune cells despite also showing evidence of epigenetic repression. We hypothesized that these competing forces of repressive and activating epigenetic marks might lead to inducible enhancer activity. We investigated a specific L1M2a element located within the first intron of Interferon Alpha/Beta Receptor 1 (IFNAR1). This element shows epigenetic signatures of B cell-specific enhancer activity, despite being repressed by the Human Silencing Hub (HUSH) complex. CRISPR deletion of the element in B lymphoblastoid cells revealed that the element acts as an enhancer that regulates both steady state and interferon-inducible expression of IFNAR1.

CONCLUSIONS

Our study experimentally demonstrates that an L1M2a element was co-opted to function as an interferon-inducible enhancer of IFNAR1, creating a feedback loop wherein IFNAR1 is transcriptionally upregulated by interferon signaling. This finding suggests that other LINEs may exhibit cryptic cell type-specific or context-dependent enhancer activity. LINEs have received less attention than ERVs in the effort to understand the contribution of transposons to the regulatory landscape of cellular genomes, but these are likely important, lineage-specific players in the rapid evolution of immune system regulatory networks and deserve further study.

Interactive visualization of whole eukaryote genome alignments using NCBI's Comparative Genome Viewer (CGV)

We report a new visualization tool for analysis of whole genome assembly-assembly alignments, the Comparative Genome Viewer (CGV) (https://ncbi.nlm.nih.gov/genome/cgv/). CGV visualizes pairwise same-species and cross-species alignments provided by NCBI using assembly alignment algorithms developed by us and others. Researchers can examine the alignments between the two assemblies using two alternate views: a chromosome ideogram-based view or a 2D genome dotplot. Whole genome alignment views expose large structural differences spanning chromosomes, such as inversions or translocations. Users can also navigate to regions of interest, where they can detect and analyze smaller-scale deletions and rearrangements within specific chromosome or gene regions. RefSeq or user-provided gene annotation is displayed in the ideogram view where available. CGV currently provides approximately 700 alignments from over 300 animal, plant, and fungal species. CGV and related NCBI viewers are undergoing active development to further meet needs of the research community in comparative genome visualization.

Chromosome-scale genome of the human blood fluke Schistosoma mekongi and its implications for public health

BACKGROUND

Schistosoma mekongi is a human blood fluke causing schistosomiasis that threatens approximately 1.5 million humans in the world. Nonetheless, the limited available S. mekongi genomic resources have hindered understanding of its biology and parasite-host interactions for disease management and pathogen control. The aim of our study was to integrate multiple technologies to construct a high-quality chromosome-level assembly of the S. mekongi genome.

METHODS

The reference genome for S. mekongi was generated through integrating Illumina, PacBio sequencing, 10 × Genomics linked-read sequencing, and high-throughput chromosome conformation capture (Hi-C) methods. In this study, we conducted de novo assembly, alignment, and gene prediction to assemble and annotate the genome. Comparative genomics allowed us to compare genomes across different species, shedding light on conserved regions and evolutionary relationships. Additionally, our transcriptomic analysis focused on genes associated with parasite-snail interactions in S. mekongi infection. We employed gene ontology (GO) enrichment analysis for functional annotation of these genes.

RESULTS

In the present study, the S. mekongi genome was both assembled into 8 pseudochromosomes with a length of 404 Mb, with contig N50 and scaffold N50 lengths of 1168 kb and 46,759 kb, respectively. We detected that 43% of the genome consists of repeat sequences and predicted 9103 protein-coding genes. We also focused on proteases, particularly leishmanolysin-like metalloproteases (M8), which are crucial in the invasion of hosts by 12 flatworm species. Through phylogenetic analysis, it was discovered that the M8 gene exhibits lineage-specific amplification among the genus Schistosoma. Lineage-specific expansion of M8 was observed in blood flukes. Additionally, the results of the RNA-seq revealed that a mass of genes related to metabolic and biosynthetic processes were up-regulated, which might be beneficial for cercaria production.

CONCLUSIONS

This study delivers a high-quality, chromosome-scale reference genome of S. mekongi, enhancing our understanding of the divergence and evolution of Schistosoma. The molecular research conducted here also plays a pivotal role in drug discovery and vaccine development. Furthermore, our work greatly advances the understanding of host-parasite interactions, providing crucial insights for schistosomiasis intervention strategies.

Long-range RNA structures in the human transcriptome beyond evolutionarily conserved regions

RNA structure has been increasingly recognized as a critical player in the biogenesis and turnover of many transcripts classes. In eukaryotes, the prediction of RNA structure by thermodynamic modeling meets fundamental limitations due to the large sizes and complex, discontinuous organization of eukaryotic genes. Signatures of functional RNA structures can be found by detecting compensatory substitutions in homologous sequences, but a comparative approach is applicable only within conserved sequence blocks. Here, we developed a computational pipeline called PHRIC, which is not limited to conserved regions and relies on RNA contacts derived from RNA in situ conformation sequencing (RIC-seq) experiments. It extracts pairs of short RNA fragments surrounded by nested clusters of RNA contacts and predicts long, nearly perfect complementary base pairings formed between these fragments. In application to a panel of RIC-seq experiments in seven human cell lines, PHRIC predicted ~12,000 stable long-range RNA structures with equilibrium free energy below -15 kcal/mol, the vast majority of which fall outside of regions annotated as conserved among vertebrates. These structures, nevertheless, show some level of sequence conservation and remarkable compensatory substitution patterns in other clades. Furthermore, we found that introns have a higher propensity to form stable long-range RNA structures between each other, and moreover that RNA structures tend to concentrate within the same intron rather than connect adjacent introns. These results for the first time extend the application of proximity ligation assays to RNA structure prediction beyond conserved regions.

Multi-genome comparisons reveal gain-and-loss evolution of the anti-Mullerian hormone receptor type 2 gene, an old master sex determining gene, in Percidae

The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis , Perca schrenkii and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene ( amhr2bY ), previously suggested to be the master sex determining (MSD) gene in P. flavescens . Phylogenetically related and structurally similar a mhr2 duplications ( amhr2b ) were found in P. schrenkii and Sander lucioperca , potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus , this amhr2b duplicate has been lost while it was subject to amplification in S. lucioperca . Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens . In P. fluviatilis , a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome-18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variants (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three ( c18h1orf198 , hsdl1 , tbc1d32 ) with higher expression in testis than ovary. Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Distinct and shared endothermic strategies in the heat producing tissues of tuna and other teleosts

Although most fishes are ectothermic, some, including tuna and billfish, achieve endothermy through specialized heat producing tissues that are modified muscles. How these heat producing tissues evolved, and whether they share convergent molecular mechanisms, remain unresolved. Here, we generated a high-quality genome from the mackerel tuna (Euthynnus affinis) and investigated the heat producing tissues of this fish by single-nucleus and bulk RNA sequencing. Compared with other teleosts, tuna-specific genetic variation is strongly associated with muscle differentiation. Single-nucleus RNA-seq revealed a high proportion of specific slow skeletal muscle cell subtypes in the heat producing tissues of tuna. Marker genes of this cell subtype are associated with the relative sliding of actin and myosin, suggesting that tuna endothermy is mainly based on shivering thermogenesis. In contrast, cross-species transcriptome analysis indicated that endothermy in billfish relies mainly on non-shivering thermogenesis. Nevertheless, the heat producing tissues of the different species do share some tissue-specific genes, including vascular-related and mitochondrial genes. Overall, although tunas and billfishes differ in their thermogenic strategies, they share similar expression patterns in some respects, highlighting the complexity of convergent evolution.

Inferring Historical Introgression with Deep Learning

Resolving phylogenetic relationships among taxa remains a challenge in the era of big data due to the presence of genetic admixture in a wide range of organisms. Rapidly developing sequencing technologies and statistical tests enable evolutionary relationships to be disentangled at a genome-wide level, yet many of these tests are computationally intensive and rely on phased genotypes, large sample sizes, restricted phylogenetic topologies, or hypothesis testing. To overcome these difficulties, we developed a deep learning-based approach, named ERICA, for inferring genome-wide evolutionary relationships and local introgressed regions from sequence data. ERICA accepts sequence alignments of both population genomic data and multiple genome assemblies, and efficiently identifies discordant genealogy patterns and exchanged regions across genomes when compared with other methods. We further tested ERICA using real population genomic data from Heliconius butterflies that have undergone adaptive radiation and frequent hybridization. Finally, we applied ERICA to characterize hybridization and introgression in wild and cultivated rice, revealing the important role of introgression in rice domestication and adaptation. Taken together, our findings demonstrate that ERICA provides an effective method for teasing apart evolutionary relationships using whole genome data, which can ultimately facilitate evolutionary studies on hybridization and introgression.

How genomic insights into the evolutionary history of clouded leopards inform their conservation

Clouded leopards (Neofelis spp.), a morphologically and ecologically distinct lineage of big cats, are severely threatened by habitat loss and fragmentation, targeted hunting, and other human activities. The long-held poor understanding of their genetics and evolution has undermined the effectiveness of conservation actions. Here, we report a comprehensive investigation of the whole genomes, population genetics, and adaptive evolution of Neofelis. Our results indicate the genus Neofelis arose during the Pleistocene, coinciding with glacial-induced climate changes to the distributions of savannas and rainforests, and signatures of natural selection associated with genes functioning in tooth, pigmentation, and tail development, associated with clouded leopards' unique adaptations. Our study highlights high-altitude adaptation as the main factor driving nontaxonomic population differentiation in Neofelis nebulosa. Population declines and inbreeding have led to reduced genetic diversity and the accumulation of deleterious variation that likely affect reproduction of clouded leopards, highlighting the urgent need for effective conservation efforts.

Genomic Signatures Associated with Transitions to Viviparity in Cyprinodontiformes

The transition from oviparity to viviparity has occurred independently over 150 times across vertebrates, presenting one of the most compelling cases of phenotypic convergence. However, whether the repeated, independent evolution of viviparity is driven by redeployment of similar genetic mechanisms and whether these leave a common signature in genomic divergence remains largely unknown. Although recent investigations into the evolution of viviparity have demonstrated striking similarity among the genes and molecular pathways involved across disparate vertebrate groups, quantitative tests for genome-wide convergent have provided ambivalent answers. Here, we investigate the potential role of molecular convergence during independent transitions to viviparity across an order of ray-finned freshwater fish (Cyprinodontiformes). We assembled de novo genomes and utilized publicly available genomes of viviparous and oviparous species to test for molecular convergence across both coding and noncoding regions. We found no evidence for an excess of molecular convergence in amino acid substitutions and in rates of sequence divergence, implying independent genetic changes are associated with these transitions. However, both statistical power and biological confounds could constrain our ability to detect significant correlated evolution. We therefore identified candidate genes with potential signatures of molecular convergence in viviparous Cyprinodontiformes lineages. Motif enrichment and gene ontology analyses suggest transcriptional changes associated with early morphogenesis, brain development, and immunity occurred alongside the evolution of viviparity. Overall, however, our findings indicate that independent transitions to viviparity in these fish are not strongly associated with an excess of molecular convergence, but a few genes show convincing evidence of convergent evolution.