RefSeq: an update on mammalian reference sequences

Cis-regulatory control of transcriptional timing and noise in response to estrogen

Cis-regulatory elements control transcription levels, temporal dynamics, and cell-cell variation or transcriptional noise. However, the combination of regulatory features that control these different attributes is not fully understood. Here, we used single-cell RNA-seq during an estrogen treatment time course and machine learning to identify predictors of expression timing and noise. We found that genes with multiple active enhancers exhibit faster temporal responses. We verified this finding by showing that manipulation of enhancer activity changes the temporal response of estrogen target genes. Analysis of transcriptional noise uncovered a relationship between promoter and enhancer activity, with active promoters associated with low noise and active enhancers linked to high noise. Finally, we observed that co-expression across single cells is an emergent property associated with chromatin looping, timing, and noise. Overall, our results indicate a fundamental tradeoff between a gene's ability to quickly respond to incoming signals and maintain low variation across cells.

A tradeoff evolution between acoustic fat bodies and skull muscles in toothed whales

Toothed whales have developed specialized echolocation abilities that are crucial for underwater activities. Acoustic fat bodies, including the melon, extramandibular fat body, and intramandibular fat body, are vital for echolocation. This study explores the transcriptome of acoustic fat bodies in toothed whales, revealing some insight into their evolutionary origins and ecological significance. Comparative transcriptome analysis of acoustic fat bodies and related tissues in a harbor porpoise and a Pacific white-sided dolphin reveals that acoustic fat bodies possess characteristics of both muscle and adipose tissue, occupying an intermediate position. The melon and extramandibular fat body exhibit specific muscle-related functions, implying an evolutionary connection between acoustic fat bodies and muscle tissue. Furthermore, we suggested that the melon and extramandibular fat body originate from intramuscular adipose tissue, a component of white adipose tissue. The extramandibular fat body has been identified as an evolutionary homolog of the masseter muscle, supported by the specific expression of MYH16, a pivotal protein in masticatory muscles. The intramandibular fat body, located within the mandibular foramen, shows possibilities of the presence of several immune-related functions, likely due to its proximity to bone marrow. Furthermore, this study sheds light on leucine modification in the catabolic pathway, which leads to the accumulation of isovaleric acid in acoustic fat bodies. Swallowing without chewing, a major toothed whale feeding ecology adaptation, makes the masticatory muscle redundant and leads to the formation of the extramandibular fat body. We propose that the intramuscular fat enlargement in facial muscles, which influences acoustic fat body development, is potentially related to the substantial reorganization of head morphology in toothed whales during aquatic adaptation.

Genetic architecture of trait variance in craniofacial morphology

The genetic architecture of trait variance has long been of interest in genetics and evolution. One of the earliest attempts to understand this architecture was presented in Lerner's Genetic Homeostasis (1954). Lerner proposed that heterozygotes should be better able to tolerate environmental perturbations because of functional differences between the alleles at a given locus, with each allele optimal for slightly different environments. This greater robustness to environmental variance, he argued, would result in smaller trait variance for heterozygotes. The evidence for Lerner's hypothesis has been inconclusive. To address this question using modern genomic methods, we mapped loci associated with differences in trait variance (vQTL) on 1,101 individuals from the F34 of an advanced intercross between LG/J and SM/J mice. We also mapped epistatic interactions for these vQTL in order to understand the influence of epistasis for the architecture of trait variance. We did not find evidence supporting Lerner's hypothesis, that heterozygotes tend to have smaller trait variances than homozygotes. We further show that the effects of most mapped loci on trait variance are produced by epistasis affecting trait means and that those epistatic effects account for about a half of the differences in genotypic-specific trait variances. Finally, we propose a model where the different interactions between the additive and dominance effects of the vQTL and their epistatic partners can explain Lerner's original observations but can also be extended to include other conditions where heterozygotes are not the least variable genotype.

MVAR: A Mouse Variation Registry

The Mouse Variation Registry (MVAR) resource is a scalable registry of mouse single nucleotide variants and small indels and variant annotation. The resource accepts data in standard Variant Call Format (VCF) and assesses the uniqueness of the submitted variants via a canonicalization process. Novel variants are assigned a unique, persistent MVAR identifier; variants that are equivalent to an existing variant in the resource are associated with the existing identifier. Annotations for variant type, molecular consequence, impact, and genomic region in the context of specific transcripts and protein sequences are generated using Ensembl's Variant Effect Predictor (VEP) and Jannovar. Access to the data and annotations in MVAR are supported via an Application Programming Interface (API) and web application. Researchers can search the resource by gene symbol, genomic region, variant (expressed in Human Genome Variation Society syntax), refSNP identifiers, or MVAR identifiers. Tabular search results can be filtered by variant annotations (variant type, molecular consequence, impact, variant region) and viewed according to variant distribution across mouse strains. The registry currently comprises more than 99 million canonical single nucleotide variants for 581 strains of mice. MVAR is accessible from https://mvar.jax.org.

Calcitriol ameliorates motor deficits and prolongs survival of Chrne-deficient mouse, a model for congenital myasthenic syndrome, by inducing Rspo2

Signal transduction at the neuromuscular junction (NMJ) is compromised in a diverse array of diseases including congenital myasthenic syndromes (CMS). Germline mutations in CHRNE encoding the acetylcholine receptor (AChR) ε subunit are the most common cause of CMS. An active form of vitamin D, calcitriol, binds to vitamin D receptor (VDR) and regulates gene expressions. We found that calcitriol enhanced MuSK phosphorylation, AChR clustering, and myotube twitching in co-cultured C2C12 myotubes and NSC34 motor neurons. RNA-seq analysis of co-cultured cells showed that calcitriol increased the expressions of Rspo2, Rapsn, and Dusp6. ChIP-seq of VDR revealed that VDR binds to a region approximately 15 ​kbp upstream to Rspo2. Biallelic deletion of the VDR-binding site of Rspo2 by CRISPR/Cas9 in C2C12 myoblasts/myotubes nullified the calcitriol-mediated induction of Rspo2 expression and MuSK phosphorylation. We generated Chrne knockout (Chrne KO) mouse by CRISPR/Cas9. Intraperitoneal administration of calcitriol markedly increased the number of AChR clusters, as well as the area, the intensity, and the number of synaptophysin-positive synaptic vesicles, in Chrne KO mice. In addition, calcitriol ameliorated motor deficits and prolonged survival of Chrne KO mice. In the skeletal muscle, calcitriol increased the gene expressions of Rspo2, Rapsn, and Dusp6. We propose that calcitriol is a potential therapeutic agent for CMS and other diseases with defective neuromuscular signal transmission.

Cis-regulatory control of transcriptional timing and noise in response to estrogen

Cis-regulatory elements control transcription levels, temporal dynamics, and cell-cell variation or transcriptional noise. However, the combination of regulatory features that control these different attributes is not fully understood. Here, we used single cell RNA-seq during an estrogen treatment time course and machine learning to identify predictors of expression timing and noise. We find that genes with multiple active enhancers exhibit faster temporal responses. We verified this finding by showing that manipulation of enhancer activity changes the temporal response of estrogen target genes. Analysis of transcriptional noise uncovered a relationship between promoter and enhancer activity, with active promoters associated with low noise and active enhancers linked to high noise. Finally, we observed that co-expression across single cells is an emergent property associated with chromatin looping, timing, and noise. Overall, our results indicate a fundamental tradeoff between a gene's ability to quickly respond to incoming signals and maintain low variation across cells.

Quality assessment of enzymatic methyl-seq library constructed using crude cell lysate

Enzymatic methyl-seq (EM-seq), an enzyme-based method, identifies genome-wide DNA methylation, which enables us to obtain reliable methylome data from purified genomic DNA by avoiding bisulfite-induced DNA damage. However, the loss of DNA during purification hinders the methylome analysis of limited samples. The crude DNA extraction method is the quickest and minimal sample loss approach for obtaining useable DNA without requiring additional dissolution and purification. However, it remains unclear whether crude DNA can be used directly for EM-seq library construction. In this study, we aimed to assess the quality of EM-seq libraries prepared directly using crude DNA. The crude DNA-derived libraries provided appropriate fragment sizes and concentrations for sequencing similar to those of the purified DNA-derived libraries. However, the sequencing results of crude samples exhibited lower reference sequence mapping efficiencies than those of the purified samples. Additionally, the lower-input crude DNA-derived sample exhibited a marginally lower cytosine-to-thymine conversion efficiency and hypermethylated pattern around gene regulatory elements than the higher-input crude DNA- or purified DNA-derived samples. In contrast, the methylation profiles of the crude and purified samples exhibited a significant correlation. Our findings indicate that crude DNA can be used as a raw material for EM-seq library construction.

Sperm DNA methylation dynamics after chemotherapy: a longitudinal study of a patient with testicular germ cell tumor treatment

BACKGROUND

An important issue for young men affected by testicular germ cell tumor (TGCT) is how TGCT and its treatment will affect, transiently or permanently, their future reproductive health. Previous studies have reported that xenobiotics can induce changes on human sperm epigenome and have the potential to promote epigenetic alterations in the offspring.

OBJECTIVES

Here, we report the first longitudinal DNA methylation profiling of frozen sperm from a TGCT patient before and up to 2 years after a bleomycin, etoposide, and cisplatin (BEP) chemotherapy.

MATERIALS AND METHODS

A TGCT was diagnosed in a 30-year-old patient. A cryopreservation of spermatozoa was proposed before adjuvant BEP treatment. Semen samples were collected before and after chemotherapy at 6, 9, 12, and 24 months. The DNA methylation status was determined by RRBS to detect DNA differentially methylated regions (DMRs).

RESULTS

The analysis revealed that among the 74 DMRs showing modified methylation status 6 months after therapy, 17 remained altered 24 months after treatment. We next associated DMRs with differentially methylated genes (DMGs), which were subsequently intersected with loci known to be important or expressed during early development.

DISCUSSION AND CONCLUSION

The consequences of the cancer treatment on the sperm epigenome during the recovery periods are topical issues of increasing significance as epigenetic modifications to the paternal genome may have deleterious effects on the offspring. The altered methylated status of these DMGs important for early development might modify their expression pattern and thus affect their function during key stages of embryogenesis, potentially leading to developmental disorders or miscarriages.

Integrase-associated niche differentiation of endogenous large DNA viruses in crustaceans

IMPORTANCE

Crustacean genomes harbor sequences originating from a family of large DNA viruses called nimaviruses, but it is unclear why they are present. We show that endogenous nimaviruses selectively insert into repetitive sequences within the host genome, and this insertion specificity was correlated with different types of integrases, which are DNA recombination enzymes encoded by the nimaviruses themselves. This suggests that endogenous nimaviruses have colonized various genomic niches through the acquisition of integrases with different insertion specificities. Our results point to a novel survival strategy of endogenous large DNA viruses colonizing the host genomes. These findings may clarify the evolution and spread of nimaviruses in crustaceans and lead to measures to control and prevent the spread of pathogenic nimaviruses in aquaculture settings.

MiniDBG: A Novel and Minimal De Bruijn Graph for Read Mapping

The De Bruijn graph (DBG) has been widely used in the algorithms for indexing or organizing read and reference sequences in bioinformatics. However, a DBG model that can locate each node, edge and path on sequence has not been proposed so far. Recently, DBG has been used for representing reference sequences in read mapping tasks. In this process, it is not a one-to-one correspondence between the paths of DBG and the substrings of reference sequence. This results in the false path on DBG, which means no substrings of reference producing the path. Moreover, if a candidate path of a read is true, we need to locate it and verify the candidate on sequence. To solve these problems, we proposed a DBG model, called MiniDBG, which stores the position lists of a minimal set of edges. With the position lists, MiniDBG can locate any node, edge and path efficiently. We also proposed algorithms for generating MiniDBG based on an original DBG and algorithms for locating edges or paths on sequence. We designed and ran experiments on real datasets for comparing them with BWT-based and position list-based methods. The experimental results show that MiniDBG can locate the edges and paths efficiently with lower memory costs.

RIP-PEN-seq identifies a class of kink-turn RNAs as splicing regulators

A kink-turn (K-turn) is a three-dimensional RNA structure that exists in all three primary phylogenetic domains. In this study, we developed the RIP-PEN-seq method to identify the full-length sequences of RNAs bound by the K-turn binding protein 15.5K and discovered a previously uncharacterized class of RNAs with backward K-turn motifs (bktRNAs) in humans and mice. All bktRNAs share two consensus sequence motifs at their fixed terminal position and have complex folding properties, expression and evolution patterns. We found that a highly conserved bktRNA1 guides the methyltransferase fibrillarin to install RNA methylation of U12 small nuclear RNA in humans. Depletion of bktRNA1 causes global splicing dysregulation of U12-type introns by impairing the recruitment of ZCRB1 to the minor spliceosome. Most bktRNAs regulate the splicing of local introns by interacting with the 15.5K protein. Taken together, our findings characterize a class of small RNAs and uncover another layer of gene expression regulation that involves crosstalk among bktRNAs, RNA splicing and RNA methylation.

Chromosome-level genome assembly of chub mackerel (Scomber japonicus) from the Indo-Pacific Ocean

Chub mackerels (Scomber japonicus) are a migratory marine fish widely distributed in the Indo-Pacific Ocean. They are globally consumed for their high Omega-3 content, but their population is declining due to global warming. Here, we generated the first chromosome-level genome assembly of chub mackerel (fScoJap1) using the Vertebrate Genomes Project assembly pipeline with PacBio HiFi genomic sequencing and Arima Hi-C chromosome contact data. The final assembly is 828.68 Mb with 24 chromosomes, nearly all containing telomeric repeats at their ends. We annotated 31,656 genes and discovered that approximately 2.19% of the genome contained DNA transposon elements repressed within duplicated genes. Analyzing 5-methylcytosine (5mC) modifications using HiFi reads, we observed open/close chromatin patterns at gene promoters, including the FADS2 gene involved in Omega-3 production. This chromosome-level reference genome provides unprecedented opportunities for advancing our knowledge of chub mackerels in biology, industry, and conservation.

World Trade Center Exposure, DNA Methylation Changes, and Cancer: A Review of Current Evidence

Introduction: Known carcinogens in the dust and fumes from the destruction of the World Trade Center (WTC) towers on 9 November 2001 included metals, asbestos, and organic pollutants, which have been shown to modify epigenetic status. Epigenome-wide association analyses (EWAS) using uniform (Illumina) methodology have identified novel epigenetic profiles of WTC exposure. Methods: We reviewed all published data, comparing differentially methylated gene profiles identified in the prior EWAS studies of WTC exposure. This included DNA methylation changes in blood-derived DNA from cases of cancer-free "Survivors" and those with breast cancer, as well as tissue-derived DNA from "Responders" with prostate cancer. Emerging molecular pathways related to the observed DNA methylation changes in WTC-exposed groups were explored and summarized. Results: WTC dust exposure appears to be associated with DNA methylation changes across the genome. Notably, WTC dust exposure appears to be associated with increased global DNA methylation; direct dysregulation of cancer genes and pathways, including inflammation and immune system dysregulation; and endocrine system disruption, as well as disruption of cholesterol homeostasis and lipid metabolism. Conclusion: WTC dust exposure appears to be associated with biologically meaningful DNA methylation changes, with implications for carcinogenesis and development of other chronic diseases.

The Complete Sequence and Comparative Analysis of Ape Sex Chromosomes

Apes possess two sex chromosomes-the male-specific Y and the X shared by males and females. The Y chromosome is crucial for male reproduction, with deletions linked to infertility. The X chromosome carries genes vital for reproduction and cognition. Variation in mating patterns and brain function among great apes suggests corresponding differences in their sex chromosome structure and evolution. However, due to their highly repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the state-of-the-art experimental and computational methods developed for the telomere-to-telomere (T2T) human genome, we produced gapless, complete assemblies of the X and Y chromosomes for five great apes (chimpanzee, bonobo, gorilla, Bornean and Sumatran orangutans) and a lesser ape, the siamang gibbon. These assemblies completely resolved ampliconic, palindromic, and satellite sequences, including the entire centromeres, allowing us to untangle the intricacies of ape sex chromosome evolution. We found that, compared to the X, ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements. This divergence on the Y arises from the accumulation of lineage-specific ampliconic regions and palindromes (which are shared more broadly among species on the X) and from the abundance of transposable elements and satellites (which have a lower representation on the X). Our analysis of Y chromosome genes revealed lineage-specific expansions of multi-copy gene families and signatures of purifying selection. In summary, the Y exhibits dynamic evolution, while the X is more stable. Finally, mapping short-read sequencing data from >100 great ape individuals revealed the patterns of diversity and selection on their sex chromosomes, demonstrating the utility of these reference assemblies for studies of great ape evolution. These complete sex chromosome assemblies are expected to further inform conservation genetics of nonhuman apes, all of which are endangered species.

Cost-effective DNA methylation profiling by FML-seq

Current methods for profiling DNA methylation require costly reagents, sequencing, and labor time. We introduce fragmentation at methylated loci and sequencing (FML-seq), a sequencing library protocol that greatly reduces all these costs. Relative to other techniques tested on the same human cell lines, FML-seq produces similar measurements of absolute and differential cytosine methylation at a fraction of the price. FML-seq enables inexpensive, high-throughput experimental designs for large-scale epigenetics research projects.

Astyanax mexicanus surface and cavefish chromosome-scale assemblies for trait variation discovery

The ability of organisms to adapt to sudden extreme environmental changes produces some of the most drastic examples of rapid phenotypic evolution. The Mexican Tetra, Astyanax mexicanus, is abundant in the surface waters of northeastern Mexico, but repeated colonizations of cave environments have resulted in the independent evolution of troglomorphic phenotypes in several populations. Here, we present three chromosome-scale assemblies of this species, for one surface and two cave populations, enabling the first whole-genome comparisons between independently evolved cave populations to evaluate the genetic basis for the evolution of adaptation to the cave environment. Our assemblies represent the highest quality of sequence completeness with predicted protein-coding and non-coding gene metrics far surpassing prior resources and, to our knowledge, all long-read assembled teleost genomes, including zebrafish. Whole genome synteny alignments show highly conserved gene order among cave forms in contrast to a higher number of chromosomal rearrangements when compared to other phylogenetically close or distant teleost species. By phylogenetically assessing gene orthology across distant branches of amniotes, we discover gene orthogroups unique to A. mexicanus. When compared to a representative surface fish genome, we find a rich amount of structural sequence diversity, defined here as the number and size of insertions and deletions as well as expanding and contracting repeats across cave forms. These new more complete genomic resources ensure higher trait resolution for comparative, functional, developmental, and genetic studies of drastic trait differences within a species.

Antiviral responses in a Jamaican fruit bat intestinal organoid model of SARS-CoV-2 infection

Bats are natural reservoirs for several zoonotic viruses, potentially due to an enhanced capacity to control viral infection. However, the mechanisms of antiviral responses in bats are poorly defined. Here we established a Jamaican fruit bat (JFB, Artibeus jamaicensis) intestinal organoid model of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection. Upon infection with SARS-CoV-2, increased viral RNA and subgenomic RNA was detected, but no infectious virus was released, indicating that JFB organoids support only limited viral replication but not viral reproduction. SARS-CoV-2 replication was associated with significantly increased gene expression of type I interferons and inflammatory cytokines. Interestingly, SARS-CoV-2 also caused enhanced formation and growth of JFB organoids. Proteomics revealed an increase in inflammatory signaling, cell turnover, cell repair, and SARS-CoV-2 infection pathways. Collectively, our findings suggest that primary JFB intestinal epithelial cells mount successful antiviral interferon responses and that SARS-CoV-2 infection in JFB cells induces protective regenerative pathways.

Genome-Wide Association Study of Fiber Diameter in Alpacas

The aim of this study was the identification of candidate genomic regions associated with fiber diameter in alpacas. DNA samples were collected from 1011 female Huacaya alpacas from two geographical Andean regions in Peru (Pasco and Puno), and three alpaca farms within each region. The samples were genotyped using an Affymetrix Custom Alpaca genotyping array containing 76,508 SNPs. After the quality controls, 960 samples and 51,742 SNPs were retained. Three association study methodologies were performed. The GWAS based on a linear model allowed us to identify 11 and 35 SNPs (-log10(p-values) > 4) using information on all alpacas and alpacas with extreme values of fiber diameter, respectively. The haplotype and marker analysis method allowed us to identify nine haplotypes with standardized haplotype heritability higher than six standard deviations. The selection signatures based on cross-population extended haplotype homozygosity (XP-EHH) allowed us to identify 180 SNPs with XP-EHH values greater than |3|. Four candidate regions with adjacent SNPs identified via two association methods of analysis are located on VPA6, VPA9, VPA29 and one chromosomally unassigned scaffold. This study represents the first analysis of alpaca whole genome association with fiber diameter, using a recently assembled alpaca SNP microarray.

Gene-based burden scores identify rare variant associations for 28 blood biomarkers

BACKGROUND

A relevant part of the genetic architecture of complex traits is still unknown; despite the discovery of many disease-associated common variants. Polygenic risk score (PRS) models are based on the evaluation of the additive effects attributable to common variants and have been successfully implemented to assess the genetic susceptibility for many phenotypes. In contrast, burden tests are often used to identify an enrichment of rare deleterious variants in specific genes. Both kinds of genetic contributions are typically analyzed independently. Many studies suggest that complex phenotypes are influenced by both low effect common variants and high effect rare deleterious variants. The aim of this paper is to integrate the effect of both common and rare functional variants for a more comprehensive genetic risk modeling.

METHODS

We developed a framework combining gene-based scores based on the enrichment of rare functionally relevant variants with genome-wide PRS based on common variants for association analysis and prediction models. We applied our framework on UK Biobank dataset with genotyping and exome data and considered 28 blood biomarkers levels as target phenotypes. For each biomarker, an association analysis was performed on full cohort using gene-based scores (GBS). The cohort was then split into 3 subsets for PRS construction and feature selection, predictive model training, and independent evaluation, respectively. Prediction models were generated including either PRS, GBS or both (combined).

RESULTS

Association analyses of the cohort were able to detect significant genes that were previously known to be associated with different biomarkers. Interestingly, the analyses also revealed heterogeneous effect sizes and directionality highlighting the complexity of the blood biomarkers regulation. However, the combined models for many biomarkers show little or no improvement in prediction accuracy compared to the PRS models.

CONCLUSION

This study shows that rare variants play an important role in the genetic architecture of complex multifactorial traits such as blood biomarkers. However, while rare deleterious variants play a strong role at an individual level, our results indicate that classical common variant based PRS might be more informative to predict the genetic susceptibility at the population level.

Two Novel Variants in MITF and PAX3 Associated With Splashed White Phenotypes in Horses

Mutations causing depigmentation are relatively common in Equus caballus (horse). Over 40 alleles in multiple genes are associated with increased white spotting (as of February 2023). The splashed white phenotype, a coat spotting pattern described as appearing like the horse has been splashed with white paint, was previously associated with variants in the PAX3 and MITF genes. Both genes encode transcription factors known to control melanocyte migration and pigmentation. We report two novel mutations, a stop-gain mutation in PAX3 (XM_005610643.3:c.927C>T, ECA6:11,196,181, EquCab3.0) and a missense mutation in a binding domain of MITF (NM_001163874.1:c.993A>T, ECA16:21,559,940, EquCab3.0), each with a strong association with increased depigmentation in Pura Raza Española horses (P = 1.144E-11, N = 30, P = 4.441E-16, N = 39 respectively). Using a quantitative method to score depigmentation, the PAX3 and MITF mutations were found to have average white scores of 25.50 and 24.45, respectively, compared to the average white coat spotting score of 1.89 in the control set. The functional impact for each mutation was predicted to be moderate to extreme (I-TASSER, SMART, Variant Effect Predictor, SIFT). We propose to designate the MITF mutant allele as Splashed White 9 and the PAX3 mutant allele as Splashed White 10 per convention.

The complete sequence of a human Y chromosome

The human Y chromosome has been notoriously difficult to sequence and assemble because of its complex repeat structure that includes long palindromes, tandem repeats and segmental duplications1-3. As a result, more than half of the Y chromosome is missing from the GRCh38 reference sequence and it remains the last human chromosome to be finished4,5. Here, the Telomere-to-Telomere (T2T) consortium presents the complete 62,460,029-base-pair sequence of a human Y chromosome from the HG002 genome (T2T-Y) that corrects multiple errors in GRCh38-Y and adds over 30 million base pairs of sequence to the reference, showing the complete ampliconic structures of gene families TSPY, DAZ and RBMY; 41 additional protein-coding genes, mostly from the TSPY family; and an alternating pattern of human satellite 1 and 3 blocks in the heterochromatic Yq12 region. We have combined T2T-Y with a previous assembly of the CHM13 genome4 and mapped available population variation, clinical variants and functional genomics data to produce a complete and comprehensive reference sequence for all 24 human chromosomes.

Proteomics as a Metrological Tool to Evaluate Genome Annotation Accuracy Following De Novo Genome Assembly: A Case Study Using the Atlantic Bottlenose Dolphin (Tursiops truncatus)

The last decade has witnessed dramatic improvements in whole-genome sequencing capabilities coupled to drastically decreased costs, leading to an inundation of high-quality de novo genomes. For this reason, the continued development of genome quality metrics is imperative. Using the 2016 Atlantic bottlenose dolphin NCBI RefSeq annotation and mass spectrometry-based proteomic analysis of six tissues, we confirmed 10,402 proteins from 4711 protein groups, constituting nearly one-third of the possible predicted proteins. Since the identification of larger proteins with more identified peptides implies reduced database fragmentation and improved gene annotation accuracy, we propose the metric NP10, which attempts to capture this quality improvement. The NP10 metric is calculated by first stratifying proteomic results by identifying the top decile (or 10th 10-quantile) of identified proteins based on the number of peptides per protein and then returns the median molecular weight of the resulting proteins. When using the 2016 versus 2012 Tursiops truncatus genome annotation to search this proteomic data set, there was a 21% improvement in NP10. This metric was further demonstrated by using a publicly available proteomic data set to compare human genome annotations from 2004, 2013 and 2016, which showed a 33% improvement in NP10. These results demonstrate that proteomics may be a useful metrological tool to benchmark genome accuracy, though there is a need for reference proteomic datasets across species to facilitate the evaluation of new de novo and existing genome.

A Chromosome-Level Reference Genome for the Black-Legged Kittiwake (Rissa tridactyla), a Declining Circumpolar Seabird

Amidst the current biodiversity crisis, the availability of genomic resources for declining species can provide important insights into the factors driving population decline. In the early 1990s, the black-legged kittiwake (Rissa tridactyla), a pelagic gull widely distributed across the arctic, subarctic, and temperate zones, suffered a steep population decline following an abrupt warming of sea surface temperature across its distribution range and is currently listed as Vulnerable by the International Union for the Conservation of Nature. Kittiwakes have long been the focus for field studies of physiology, ecology, and ecotoxicology and are primary indicators of fluctuating ecological conditions in arctic and subarctic marine ecosystems. We present a high-quality chromosome-level reference genome and annotation for the black-legged kittiwake using a combination of Pacific Biosciences HiFi sequencing, Bionano optical maps, Hi-C reads, and RNA-Seq data. The final assembly spans 1.35 Gb across 32 chromosomes, with a scaffold N50 of 88.21 Mb and a BUSCO completeness of 97.4%. This genome assembly substantially improves the quality of a previous draft genome, showing an approximately 5× increase in contiguity and a more complete annotation. Using this new chromosome-level reference genome and three more chromosome-level assemblies of Charadriiformes, we uncover several lineage-specific chromosome fusions and fissions, but find no shared rearrangements, suggesting that interchromosomal rearrangements have been commonplace throughout the diversification of Charadriiformes. This new high-quality genome assembly will enable population genomic, transcriptomic, and phenotype-genotype association studies in a widely studied sentinel species, which may provide important insights into the impacts of global change on marine systems.

Refining the genomic determinants underlying escape from X-chromosome inactivation

X-chromosome inactivation (XCI) epigenetically silences one X chromosome in every cell in female mammals. Although the majority of X-linked genes are silenced, in humans 20% or more are able to escape inactivation and continue to be expressed. Such escape genes are important contributors to sex differences in gene expression, and may impact the phenotypes of X aneuploidies; yet the mechanisms regulating escape from XCI are not understood. We have performed an enrichment analysis of transcription factor binding on the X chromosome, providing new evidence for enriched factors at the transcription start sites of escape genes. The top escape-enriched transcription factors were detected at the RPS4X promoter, a well-described human escape gene previously demonstrated to escape from XCI in a transgenic mouse model. Using a cell line model system that allows for targeted integration and inactivation of transgenes on the mouse X chromosome, we further assessed combinations of RPS4X promoter and genic elements for their ability to drive escape from XCI. We identified a small transgenic construct of only 6 kb capable of robust escape from XCI, establishing that gene-proximal elements are sufficient to permit escape, and highlighting the additive effect of multiple elements that work together in a context-specific fashion.

5'UTR Variant in KIT Associated with White Spotting in Horses

Mutations in KIT, a gene that influences melanoblast migration and pigmentation, often result in mammalian white spotting. As of February 2023, over 30 KIT variants associated with white spotting were documented in Equus caballus (horse). Here we report an association of increased white spotting on the skin and coat with a variant in the 5'UTR of KIT (rs1149701677: g.79,618,649A>C). Horses possessing at least one alternate allele demonstrate phenotypic characteristics similar to other KIT mutations: clear borders around unpigmented regions on the body, face, and limbs. Using a quantitative measure of depigmentation, we observed an average white score of 10.70 among individuals with rs1149701677, while the average score of the control, homozygous reference sample was 2.23 (p=1.892e-11, n=109, t-test). The rs1149701677 site has a cross-species conservation score of 3.4, one of the highest scores across the KIT 5'UTR, implying regulatory importance for this site. Ensembl also predicted a "moderately impactful" functional effect for the rs1149701677 variant. We propose that this single nucleotide variant likely alters the regulation of KIT, which in turn may disrupt melanoblast migration causing an increase in white spotting on the coat. Alternatively, the rs1149701677 variant may be in linkage with another nearby variant with an as-yet-undiscovered functional impact. We propose to term this new allele "Holiday White" or W35 based on conventional nomenclature.

Evolutionary states and trajectories characterized by distinct pathways stratify patients with ovarian high grade serous carcinoma

Ovarian high-grade serous carcinoma (HGSC) is typically diagnosed at an advanced stage, with multiple genetically heterogeneous clones existing in the tumors long before therapeutic intervention. Herein we integrate clonal composition and topology using whole-genome sequencing data from 510 samples of 148 patients with HGSC in the prospective, longitudinal, multiregion DECIDER study. Our results reveal three evolutionary states, which have distinct features in genomics, pathways, and morphological phenotypes, and significant association with treatment response. Nested pathway analysis suggests two evolutionary trajectories between the states. Experiments with five tumor organoids and three PI3K inhibitors support targeting tumors with enriched PI3K/AKT pathway with alpelisib. Heterogeneity analysis of samples from multiple anatomical sites shows that site-of-origin samples have 70% more unique clones than metastatic tumors or ascites. In conclusion, these analysis and visualization methods enable integrative tumor evolution analysis to identify patient subtypes using data from longitudinal, multiregion cohorts.

Sequencing of 19,219 exomes identifies a low-frequency variant in FKBP5 promoter predisposing to high myopia in a Han Chinese population

High myopia (HM) is one of the leading causes of visual impairment and blindness worldwide. Here, we report a whole-exome sequencing (WES) study in 9,613 HM cases and 9,606 controls of Han Chinese ancestry to pinpoint HM-associated risk variants. Single-variant association analysis identified three newly identified -genetic loci associated with HM, including an East Asian ancestry-specific low-frequency variant (rs533280354) in FKBP5. Multi-ancestry meta-analysis with WES data of 2,696 HM cases and 7,186 controls of European ancestry from the UK Biobank discerned a newly identified European ancestry-specific rare variant in FOLH1. Functional experiments revealed a mechanism whereby a single G-to-A transition at rs533280354 disrupted the binding of transcription activator KLF15 to the promoter of FKBP5, resulting in decreased transcription of FKBP5. Furthermore, burden tests showed a significant excess of rare protein-truncating variants among HM cases involved in retinal blood vessel morphogenesis and neurotransmitter transport.

Twist exome capture allows for lower average sequence coverage in clinical exome sequencing

BACKGROUND

Exome and genome sequencing are the predominant techniques in the diagnosis and research of genetic disorders. Sufficient, uniform and reproducible/consistent sequence coverage is a main determinant for the sensitivity to detect single-nucleotide (SNVs) and copy number variants (CNVs). Here we compared the ability to obtain comprehensive exome coverage for recent exome capture kits and genome sequencing techniques.

RESULTS

We compared three different widely used enrichment kits (Agilent SureSelect Human All Exon V5, Agilent SureSelect Human All Exon V7 and Twist Bioscience) as well as short-read and long-read WGS. We show that the Twist exome capture significantly improves complete coverage and coverage uniformity across coding regions compared to other exome capture kits. Twist performance is comparable to that of both short- and long-read whole genome sequencing. Additionally, we show that even at a reduced average coverage of 70× there is only minimal loss in sensitivity for SNV and CNV detection.

CONCLUSION

We conclude that exome sequencing with Twist represents a significant improvement and could be performed at lower sequence coverage compared to other exome capture techniques.

High resolution genomes of multiple Xiphophorus species provide new insights into microevolution, hybrid incompatibility, and epistasis

Because of diverged adaptative phenotypes, fish species of the genus Xiphophorus have contributed to a wide range of research for a century. Existing Xiphophorus genome assemblies are not at the chromosomal level and are prone to sequence gaps, thus hindering advancement of the intra- and inter-species differences for evolutionary, comparative, and translational biomedical studies. Herein, we assembled high-quality chromosome-level genome assemblies for three distantly related Xiphophorus species, namely, X. maculatus, X. couchianus, and X. hellerii Our overall goal is to precisely assess microevolutionary processes in the clade to ascertain molecular events that led to the divergence of the Xiphophorus species and to progress understanding of genetic incompatibility to disease. In particular, we measured intra- and inter-species divergence and assessed gene expression dysregulation in reciprocal interspecies hybrids among the three species. We found expanded gene families and positively selected genes associated with live bearing, a special mode of reproduction. We also found positively selected gene families are significantly enriched in nonpolymorphic transposable elements, suggesting the dispersal of these nonpolymorphic transposable elements has accompanied the evolution of the genes, possibly by incorporating new regulatory elements in support of the Britten-Davidson hypothesis. We characterized inter-specific polymorphisms, structural variants, and polymorphic transposable element insertions and assessed their association to interspecies hybridization-induced gene expression dysregulation related to specific disease states in humans.

Resequencing of a Pekin duck breeding population provides insights into the genomic response to short-term artificial selection

BACKGROUND

Short-term, intense artificial selection drives fast phenotypic changes in domestic animals and leaves imprints on their genomes. However, the genetic basis of this selection response is poorly understood. To better address this, we employed the Pekin duck Z2 pure line, in which the breast muscle weight was increased nearly 3-fold after 10 generations of breeding. We denovo assembled a high-quality reference genome of a female Pekin duck of this line (GCA_003850225.1) and identified 8.60 million genetic variants in 119 individuals among 10 generations of the breeding population.

RESULTS

We identified 53 selected regions between the first and tenth generations, and 93.8% of the identified variations were enriched in regulatory and noncoding regions. Integrating the selection signatures and genome-wide association approach, we found that 2 regions covering 0.36 Mb containing UTP25 and FBRSL1 were most likely to contribute to breast muscle weight improvement. The major allele frequencies of these 2 loci increased gradually with each generation following the same trend. Additionally, we found that a copy number variation region containing the entire EXOC4 gene could explain 1.9% of the variance in breast muscle weight, indicating that the nervous system may play a role in economic trait improvement.

CONCLUSIONS

Our study not only provides insights into genomic dynamics under intense artificial selection but also provides resources for genomics-enabled improvements in duck breeding.

Divergent sensory and immune gene evolution in sea turtles with contrasting demographic and life histories

Sea turtles represent an ancient lineage of marine vertebrates that evolved from terrestrial ancestors over 100 Mya. The genomic basis of the unique physiological and ecological traits enabling these species to thrive in diverse marine habitats remains largely unknown. Additionally, many populations have drastically declined due to anthropogenic activities over the past two centuries, and their recovery is a high global conservation priority. We generated and analyzed high-quality reference genomes for the leatherback (Dermochelys coriacea) and green (Chelonia mydas) turtles, representing the two extant sea turtle families. These genomes are highly syntenic and homologous, but localized regions of noncollinearity were associated with higher copy numbers of immune, zinc-finger, and olfactory receptor (OR) genes in green turtles, with ORs related to waterborne odorants greatly expanded in green turtles. Our findings suggest that divergent evolution of these key gene families may underlie immunological and sensory adaptations assisting navigation, occupancy of neritic versus pelagic environments, and diet specialization. Reduced collinearity was especially prevalent in microchromosomes, with greater gene content, heterozygosity, and genetic distances between species, supporting their critical role in vertebrate evolutionary adaptation. Finally, diversity and demographic histories starkly contrasted between species, indicating that leatherback turtles have had a low yet stable effective population size, exhibit extremely low diversity compared with other reptiles, and harbor a higher genetic load compared with green turtles, reinforcing concern over their persistence under future climate scenarios. These genomes provide invaluable resources for advancing our understanding of evolution and conservation best practices in an imperiled vertebrate lineage.

A chromosome-level reference genome and pangenome for barn swallow population genomics

Insights into the evolution of non-model organisms are limited by the lack of reference genomes of high accuracy, completeness, and contiguity. Here, we present a chromosome-level, karyotype-validated reference genome and pangenome for the barn swallow (Hirundo rustica). We complement these resources with a reference-free multialignment of the reference genome with other bird genomes and with the most comprehensive catalog of genetic markers for the barn swallow. We identify potentially conserved and accelerated genes using the multialignment and estimate genome-wide linkage disequilibrium using the catalog. We use the pangenome to infer core and accessory genes and to detect variants using it as a reference. Overall, these resources will foster population genomics studies in the barn swallow, enable detection of candidate genes in comparative genomics studies, and help reduce bias toward a single reference genome.

Genome-wide analysis of copy-number variation in humans with cleft lip and/or cleft palate identifies COBLL1, RIC1, and ARHGEF38 as clefting genes

Cleft lip with or without cleft palate (CL/P) is a common birth defect with a complex, heterogeneous etiology. It is well established that common and rare sequence variants contribute to the formation of CL/P, but the contribution of copy-number variants (CNVs) to cleft formation remains relatively understudied. To fill this knowledge gap, we conducted a large-scale comparative analysis of genome-wide CNV profiles of 869 individuals from the Philippines and 233 individuals of European ancestry with CL/P with three primary goals: first, to evaluate whether differences in CNV number, amount of genomic content, or amount of coding genomic content existed within clefting subtypes; second, to assess whether CNVs in our cohort overlapped with known Mendelian clefting loci; and third, to identify unestablished Mendelian clefting genes. Significant differences in CNVs across cleft types or in individuals with non-syndromic versus syndromic clefts were not observed; however, several CNVs in our cohort overlapped with known syndromic and non-syndromic Mendelian clefting loci. Moreover, employing a filtering strategy relying on population genetics data that rare variants are on the whole more deleterious than common variants, we identify several CNV-associated gene losses likely driving non-syndromic clefting phenotypes. By prioritizing genes deleted at a rare frequency across multiple individuals with clefts yet enriched in our cohort of individuals with clefts compared to control subjects, we identify COBLL1, RIC1, and ARHGEF38 as clefting genes. CRISPR-Cas9 mutagenesis of these genes in Xenopus laevis and Danio rerio yielded craniofacial dysmorphologies, including clefts analogous to those seen in human clefting disorders.

A phenotype driven integrative framework uncovers molecular mechanisms of a rare hereditary thrombophilia

Antithrombin resistance is a rare subtype of hereditary thrombophilia caused by prothrombin gene variants, leading to thrombotic disorders. Recently, the Prothrombin Belgrade variant has been reported as a specific variant that leads to antithrombin resistance in two Serbian families with thrombosis. However, due to clinical data scarcity and the inapplicability of traditional genome-wide association studies (GWAS), a broader perspective on molecular and phenotypic mechanisms associated with the Prothrombin Belgrade variant is yet to be uncovered. Here, we propose an integrative framework to address the lack of genomic samples and support the genomic signal from the full genome sequences of five heterozygous subjects by integrating it with subjects' phenotypes and the genes' molecular interactions. Our goal is to identify candidate thrombophilia-related genes for which our subjects possess germline variants by focusing on the resulting gene clusters of our integrative framework. We applied a Non-negative Matrix Tri-Factorization-based method to simultaneously integrate different data sources, taking into account the observed phenotypes. In other words, our data-integration framework reveals gene clusters involved with this rare disease by fusing different datasets. Our results are in concordance with the current literature about antithrombin resistance. We also found candidate disease-related genes that need to be further investigated. CD320, RTEL1, UCP2, APOA5 and PROZ participate in healthy-specific or disease-specific subnetworks involving thrombophilia-annotated genes and are related to general thrombophilia mechanisms according to the literature. Moreover, the ADRA2A and TBXA2R subnetworks analysis suggested that their variants may have a protective effect due to their connection with decreased platelet activation. The results show that our method can give insights into antithrombin resistance even if a small amount of genetic data is available. Our framework is also customizable, meaning that it applies to any other rare disease.

Clustered variants in the 5' coding region of TRA2B cause a distinctive neurodevelopmental syndrome

PURPOSE

Transformer2 proteins (Tra2α and Tra2β) control splicing patterns in human cells, and no human phenotypes have been associated with germline variants in these genes. The aim of this work was to associate germline variants in the TRA2B gene to a novel neurodevelopmental disorder.

METHODS

A total of 12 individuals from 11 unrelated families who harbored predicted loss-of-function monoallelic variants, mostly de novo, were recruited. RNA sequencing and western blot analyses of Tra2β-1 and Tra2β-3 isoforms from patient-derived cells were performed. Tra2β1-GFP, Tra2β3-GFP and CHEK1 exon 3 plasmids were transfected into HEK-293 cells.

RESULTS

All variants clustered in the 5' part of TRA2B, upstream of an alternative translation start site responsible for the expression of the noncanonical Tra2β-3 isoform. All affected individuals presented intellectual disability and/or developmental delay, frequently associated with infantile spasms, microcephaly, brain anomalies, autism spectrum disorder, feeding difficulties, and short stature. Experimental studies showed that these variants decreased the expression of the canonical Tra2β-1 isoform, whereas they increased the expression of the Tra2β-3 isoform, which is shorter and lacks the N-terminal RS1 domain. Increased expression of Tra2β-3-GFP were shown to interfere with the incorporation of CHEK1 exon 3 into its mature transcript, normally incorporated by Tra2β-1.

CONCLUSION

Predicted loss-of-function variants clustered in the 5' portion of TRA2B cause a new neurodevelopmental syndrome through an apparently dominant negative disease mechanism involving the use of an alternative translation start site and the overexpression of a shorter, repressive Tra2β protein.

A chromosome-level genome assembly reveals genomic characteristics of the American mink (Neogale vison)

Availability of a contiguous chromosome-level genome assembly is the foundational step to develop genome-based studies in American mink (Neogale vison). The main objective of this study was to provide a high quality chromosome-level genome assembly for American mink. An initial draft of the genome assembly was generated using 2,884,047 PacBio long reads. Integration of Hi-C data into the initial draft led to an assembly with 183 scaffolds and scaffold N50 of 220 Mb. This gap-free genome assembly of American mink (ASM_NN_V1) had a length of 2.68 Gb in which about 98.6% of the whole genome was covered by 15 chromosomes. In total, 25,377 genes were predicted across the American mink genome using the NCBI Eukaryotic Genome Annotation Pipeline. In addition, gene orthology, demographic history, synteny blocks, and phylogenetic relationships were studied in connection with the genomes of other related Carnivora. Furthermore, population-based statistics of 100 sequenced mink were presented using the newly assembled genome. Remarkable improvements were observed in genome contiguity, the number of scaffolds, and annotation compared to the first draft of mink genome assembly (NNQGG.v01). This high-quality genome assembly will support the development of efficient breeding strategies as well as conservation programs for American mink.

Antiviral response mechanisms in a Jamaican Fruit Bat intestinal organoid model of SARS-CoV-2 infection

Bats are natural reservoirs for several zoonotic viruses, potentially due to an enhanced capacity to control viral infection. However, the mechanisms of antiviral responses in bats are poorly defined. Here we established a Jamaican fruit bat (JFB) intestinal organoid model of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection. JFB organoids were susceptible to SARS-CoV-2 infection, with increased viral RNA and subgenomic RNA detected in cell lysates and supernatants. Gene expression of type I interferons and inflammatory cytokines was induced in response to SARS-CoV-2 but not in response to TLR agonists. Interestingly, SARS-CoV-2 did not lead to cytopathic effects in JFB organoids but caused enhanced organoid growth. Proteomic analyses revealed an increase in inflammatory signaling, cell turnover, cell repair, and SARS-CoV-2 infection pathways. Collectively, our findings suggest that primary JFB intestinal epithelial cells can mount a successful antiviral interferon response and that SARS-CoV-2 infection in JFB cells induces protective regenerative pathways.

Role of lncRNA LIPE-AS1 in adipogenesis

Recent studies have identified long non-coding RNAs (lncRNAs) as potential regulators of adipogenesis. In this study, we have characterized a lncRNA, LIPE-AS1, that spans genes CEACAM1 to LIPE in man with conservation of genomic organization and tissue expression between mouse and man. Tissue-specific expression of isoforms of the murine lncRNA were found in liver and adipose tissue, one of which, designated mLas-V3, overlapped the Lipe gene encoding hormone-sensitive lipase in both mouse and man suggesting that it may have a functional role in adipose tissue. Knock down of expression of mLas-V3 using anti-sense oligos (ASOs) led to a significant decrease in the differentiation of the OP9 pre-adipocyte cell line through the down regulation of the major adipogenic transcription factors Pparg and Cebpa. Knock down of mLas-V3 induced apoptosis during the differentiation of OP9 cells as shown by expression of active caspase-3, a change in the localization of LIP/LAP isoforms of C/EBPβ, and expression of the cellular stress induced factors CHOP, p53, PUMA, and NOXA. We conclude that mLas-V3 may play a role in protecting against stress associated with adipogenesis, and its absence leads to apoptosis.

Transcriptomic analysis reveals distinct mechanisms of adaptation of a polar picophytoplankter under ocean acidification conditions

Human emissions of carbon dioxide are causing irreversible changes in our oceans and impacting marine phytoplankton, including a group of small green algae known as picochlorophytes. Picochlorophytes grown in natural phytoplankton communities under future predicted levels of carbon dioxide have been demonstrated to thrive, along with redistribution of the cellular metabolome that enhances growth rate and photosynthesis. Here, using next-generation sequencing technology, we measured levels of transcripts in a picochlorophyte Chlorella, isolated from the sub-Antarctic and acclimated under high and current ambient CO2 levels, to better understand the molecular mechanisms involved with its ability to acclimate to elevated CO2. Compared to other phytoplankton taxa that induce broad transcriptomic responses involving multiple parts of their cellular metabolism, the changes observed in Chlorella focused on activating gene regulation involved in different sets of pathways such as light harvesting complex binding proteins, amino acid synthesis and RNA modification, while carbon metabolism was largely unaffected. Triggering a specific set of genes could be a unique strategy of small green phytoplankton under high CO2 in polar oceans.

csORF-finder: an effective ensemble learning framework for accurate identification of multi-species coding short open reading frames

Short open reading frames (sORFs) refer to the small nucleic fragments no longer than 303 nt in length that probably encode small peptides. To date, translatable sORFs have been found in both untranslated regions of messenger ribonucleic acids (RNAs; mRNAs) and long non-coding RNAs (lncRNAs), playing vital roles in a myriad of biological processes. As not all sORFs are translated or essentially translatable, it is important to develop a highly accurate computational tool for characterizing the coding potential of sORFs, thereby facilitating discovery of novel functional peptides. In light of this, we designed a series of ensemble models by integrating Efficient-CapsNet and LightGBM, collectively termed csORF-finder, to differentiate the coding sORFs (csORFs) from non-coding sORFs in Homo sapiens, Mus musculus and Drosophila melanogaster, respectively. To improve the performance of csORF-finder, we introduced a novel feature encoding scheme named trinucleotide deviation from expected mean (TDE) and computed all types of in-frame sequence-based features, such as i-framed-3mer, i-framed-CKSNAP and i-framed-TDE. Benchmarking results showed that these features could significantly boost the performance compared to the original 3-mer, CKSNAP and TDE features. Our performance comparisons showed that csORF-finder achieved a superior performance than the state-of-the-art methods for csORF prediction on multi-species and non-ATG initiation independent test datasets. Furthermore, we applied csORF-finder to screen the lncRNA datasets for identifying potential csORFs. The resulting data serve as an important computational repository for further experimental validation. We hope that csORF-finder can be exploited as a powerful platform for high-throughput identification of csORFs and functional characterization of these csORFs encoded peptides.

Cytotoxic T Cell Expression of Leukocyte-Associated Immunoglobulin-Like Receptor-1 (LAIR-1) in Viral Hepatitis C-Mediated Hepatocellular Carcinoma

Virus-related hepatocellular carcinoma (HCC) pathogenesis involves liver inflammation, therefore, despite successful treatment, hepatitis C virus (HCV) may progress to HCC from initiated liver cirrhosis. Cytotoxic T cells (Tcs) are known to be involved in HCV-related cirrhotic complications and HCC pathogenesis. The inhibitory checkpoint leukocyte-associated immunoglobulin-like receptor-1 (LAIR-1) is expressed on Tcs. Therefore, we aimed to determine whether the Tc expression level of LAIR-1 is associated with HCC progression and to evaluate LAIR-1 expression as a noninvasive biomarker for HCC progression in the context of liver cirrhosis related to HCV genotype 4 (G4) in Egyptian patients’ peripheral venous blood liquid biopsy. A total of 64 patients with HCC and 37 patients with liver cirrhosis were enrolled in this case-controlled study, and their LAIR-1 expression on Tc related to the progression of liver cirrhosis was examined and compared to that of the apparently healthy control group (n = 20). LAIR-1 expression was analyzed using flow cytometry. Results: The HCC group had significantly higher LAIR-1 expression on Tc and percentage of Tc positive for LAIR-1 (LAIR-1+Tc%) than the HCV G4-related liver cirrhosis group. LAIR-1+Tc% was correlated with the HCC surrogate tumor marker AFP (r = 0.367, p = 0.001) and insulin resistance and inflammation prognostic ratios/indices. A receiver operating characteristic (ROC) curve revealed that adding LAIR-1+Tc% to AFP can distinguish HCC transformation in the Egyptian patients’ cohort. Upregulated LAIR-1 expression on Tc could be a potential screening noninvasive molecular marker for chronic inflammatory HCV G4 related liver cirrhosis. Moreover, LAIR-1 expression on Tc may be one of the players involved in the progression of liver cirrhosis to HCC.

EfGD: the Erianthus fulvus genome database

Erianthus fulvus (TaxID: 154759) is a valuable germplasm resource in sugarcane breeding and research and has excellent agronomic traits, such as drought resistance, cold resistance, barren tolerance and high brix. With a stable chromosome number (2n = 20) and a small genome (0.9 Gb), it is an ideal candidate for research on sugarcane. Next-generation sequencing technology has enabled a growing number of studies to focus on genomics. Due to the large amount of omics data available, a centralized platform is necessary for ensuring the consistency, independence and maintainability of these large-scale datasets through storage, analysis and integration. Here, we present a comprehensive database for the E. fulvus genome, EfGD. By using the new high-quality reference genome and its annotations, the EfGD provides the largest whole-genome sequencing reference dataset for E. fulvus, which archives 27 165 protein-coding genes and 55 564 488 SNPs from 202 newly resequenced genomes. Furthermore, we created a user-friendly graphical interface for visualizing genomic diversity, population structure and evolution and provided other tools on an open platform.

Database URL: https://efgenome.ynau.edu.cn

Genome-wide association meta-analysis of 88,250 individuals highlights pleiotropic mechanisms of five ocular diseases in UK Biobank

BACKGROUND

Ocular diseases may exhibit common clinical symptoms and epidemiological comorbidity. However, the extent of pleiotropic mechanisms across ocular diseases remains unclear. We aim to examine shared genetic etiology in age-related macular degeneration (AMD), diabetic retinopathy (DR), glaucoma, retinal detachment (RD), and myopia.

METHODS

We analyzed genome-wide association analyses for the five ocular diseases in 43,877 cases and 44,373 controls of European ancestry from UK Biobank, estimated their genetic relationships (LDSC, GNOVA, and Genomic SEM), and identified pleiotropic loci (ASSET and METASOFT).

FINDINGS

The genetic correlation of common SNPs revealed a meaningful genetic structure within these diseases, identifying genetic correlations between AMD, DR, and glaucoma. Cross-trait meta-analysis identified 23 pleiotropic loci associated with at least two ocular diseases and 14 loci unique to individual disorders (non-pleiotropic). We found that the genes associated with these shared genetic loci are involved in neuron differentiation (P = 8.80 × 10^-6) and eye development systems (P = 3.86 × 10^-5), and single cell RNA sequencing data reveals their heightened gene expression from multipotent progenitors to other differentiated retinal cells during retina developmental process.

INTERPRETATION

These results highlighted the potential common genetic architectures among these ocular diseases and can deepen the understanding of the molecular mechanisms underlying the related diseases.

FUNDING

The National Natural Science Foundation of China (61871294), Zhejiang Provincial Natural Science Foundation of China (LR19C060001), and the Scientific Research Foundation for Talents of Wenzhou Medical University (QTJ18023).

Web-Based Protein Interactions Calculator Identifies Likely Proteome Coevolution with Alzheimer’s Disease-Associated Proteins

Protein–protein functional interactions arise from either transitory or permanent biomolecular associations and often lead to the coevolution of the interacting residues. Although mutual information has traditionally been used to identify coevolving residues within the same protein, its application between coevolving proteins remains largely uncharacterized. Therefore, we developed the Protein Interactions Calculator (PIC) to efficiently identify coevolving residues between two protein sequences using mutual information. We verified the algorithm using 2102 known human protein interactions and 233 known bacterial protein interactions, with a respective 1975 and 252 non-interacting protein controls. The average PIC score for known human protein interactions was 4.5 times higher than non-interacting proteins (p = 1.03 × 10⁻¹⁰⁸) and 1.94 times higher in bacteria (p = 1.22 × 10⁻³⁵). We then used the PIC scores to determine the probability that two proteins interact. Using those probabilities, we paired 37 Alzheimer’s disease-associated proteins with 8608 other proteins and determined the likelihood that each pair interacts, which we report through a web interface. The PIC had significantly higher sensitivity and residue-specific resolution not available in other algorithms. Therefore, we propose that the PIC can be used to prioritize potential protein interactions, which can lead to a better understanding of biological processes and additional therapeutic targets belonging to protein interaction groups.

Drug-induced epigenomic plasticity reprograms circadian rhythm regulation to drive prostate cancer towards androgen-independence

In prostate cancer, androgen receptor (AR)-targeting agents are very effective in various disease stages. However, therapy resistance inevitably occurs and little is known about how tumor cells adapt to bypass AR suppression. Here, we performed integrative multi-omics analyses on tissues isolated before and after 3 months of AR-targeting enzalutamide monotherapy from high-risk prostate cancer patients enrolled in a neoadjuvant clinical trial. Transcriptomic analyses demonstrated that AR inhibition drove tumors towards a neuroendocrine-like disease state. Additionally, epigenomic profiling revealed massive enzalutamide-induced reprogramming of pioneer factor FOXA1 - from inactive chromatin sites towards active cis-regulatory elements that dictate pro-survival signals. Notably, treatment-induced FOXA1 sites were enriched for circadian clock component ARNTL. Post-treatment ARNTL levels associated with poor outcome, and ARNTL knockout strongly decreased prostate cancer cell growth. Our data highlight a remarkable cistromic plasticity of FOXA1 following AR-targeted therapy, and revealed an acquired dependency on circadian regulator ARNTL, a novel candidate therapeutic target.

Mutations in RPS19 may affect ribosome function and biogenesis in Diamond Blackfan Anemia

Ribosomes, the cellular organelles translating the genetic code to proteins, are assemblies of RNA chains and many proteins (RPs) arranged in precise fine-tuned interwoven structures. Mutated ribosomal genes cause ribosomopathies, including Diamond Blackfan Anemia (DBA, a rare heterogeneous red-cell aplasia connected to ribosome malfunction) or failed biogenesis. Combined bioinformatical, structural, and predictive analyses of potential consequences of possibly expressed mutations in eS19, the protein product of the highly mutated RPS19, suggests that mutations in its exposed surface could alter its positioning during assembly and consequently prevent biogenesis, implying a natural selective strategy to avoid malfunctions in ribosome assembly. A search for RPS19 pseudogenes indicated >90% sequence identity with the wild type, hinting at its expression in cases of absent or truncated gene products.

APOBEC3A regulates transcription from interferon-stimulated response elements

APOBEC3A (A3A) is a cytidine deaminase that inactivates a variety of viruses through introduction of lethal mutations to the viral genome. Additionally, A3A can suppress HIV-1 transcription in a deaminase-independent manner by binding to the long terminal repeat of proviral HIV-1. However, it is unknown whether A3A targets additional host genomic loci for repression. In this study, we found that A3A suppresses gene expression by binding TTTC doublets that are in close proximity to each other. However, one TTTC motif is sufficient for A3A binding. Because TTTC doublets are present in interferon (IFN)-stimulated response elements (ISRE), we hypothesized that A3A may impact IFN-stimulated gene (ISG) expression. After scanning the human genome for TTTC doublet occurrences, we discovered that these motifs are enriched in the proximal promoters of genes associated with antiviral responses and type I IFN (IFN-I) signaling. As a proof of principle, we examined whether A3A can impact ISG15 expression. We found that A3A binding to the ISRE inhibits phosphorylated STAT-1 binding and suppresses ISG15 induction in response to IFN-I treatment. Consistent with these data, our RNA-sequencing analyses indicate that A3A loss results in increased IFN-I–dependent induction of several ISGs. This study revealed that A3A plays an unexpected role in ISG regulation and suggests that A3A contributes to a negative feedback loop during IFN signaling.

Downregulation of CRTC1 Is Involved in CUMS-Induced Depression-Like Behavior in the Hippocampus and Its RNA Sequencing Analysis

Chronic stress is an important risk factor for mood disorders including depression. The decreased level of CREB (cAMP-responsive element binding)-regulated transcription coactivator 1 (CRTC1) expression in hippocampus may be involved in depression-like behavior in some stress-induced depression models. But the mechanism of CRTC1 in mediating depression-like behavior remains unknown. In this study, chronic unpredictable mild stress (CUMS)-treated mice showed depression-like behavior accompanied by the downregulation of CRTC1 in the hippocampus. Adeno-associated virus (AAV)-CRTC1-mediated overexpression of CRTC1 in the hippocampus by stereotactic brain injection could significantly prevent depression-like behavior in CUMS-treated mice. The above data reveal that the downregulation of hippocampal CRTC1 expression participates in CUMS-induced depression-like behavior. In order to explore the key targets regulated by CRTC1, AAV-mediated CRTC1 short hairpin (shRNA) was constructed to achieve knockdown of CRTC1 in the hippocampus, and then the hippocampi were collected for RNA-sequencing (RNA-seq). The RNA-seq data show that upregulated genes were enriched in stress and immune system-associated GO terms and pathways such as response to stress and external stimulus and regulation of immune response and that downregulated genes were enriched in neural activity such as synaptic transmission and cognitive behavior. We further provided RT-qPCR data that the inflammation-related factors including Gpr84, Tlr2, Lyz2, and Icam1 were significantly upregulated in the hippocampus of both CUMS- and CRTC1 shRNA-induced models, some of them were also validated in protein levels by Western blotting. We propose a hypothesis that CUMS induces downregulation of CRTC1, which might lead to depression-like behavior via neuroinflammation pathway. This study provides new explanation for the inflammatory hypothesis of depression and some clues for exploring the molecular mechanism of CRTC1 regulation.

Human disease prediction from microbiome data by multiple feature fusion and deep learning

Human disease prediction from microbiome data has broad implications in metagenomics. It is rare for the existing methods to consider abundance profiles from both known and unknown microbial organisms, or capture the taxonomic relationships among microbial taxa, leading to significant information loss. On the other hand, deep learning has shown unprecedented advantages in classification tasks for its feature-learning ability. However, it encounters the opposite situation in metagenome-based disease prediction since high-dimensional low-sample-size metagenomic datasets can lead to severe overfitting; and black-box model fails in providing biological explanations. To circumvent the related problems, we developed MetaDR, a comprehensive machine learning-based framework that integrates various information and deep learning to predict human diseases. Experimental results indicate that MetaDR achieves competitive prediction performance with a reduction in running time, and effectively discovers the informative features with biological insights.

Establishment and characterization of the third non-functional human pancreatic neuroendocrine tumor cell line

The mechanisms of neuroendocrine tumor (NET) were still poorly understood, largely due to the lack of preclinical models of neuroendocrine neoplasms. Herein, we established and characterized SPNE1 cell lines from primary pancreatic NET tissue obtained from a 44-year-old female. Neuroendocrine character of SPNE1 was compared with existing non-functional cell lines BON1 and QGP1, and the results indicated expressions of multiple NET-specific markers in SPNE1 were higher relative to BON1 and QGP1. The growth character measured by Ki67 labeling index, cell cycle analysis, and 3D matrigel spheroid essay indicated that the proliferative rate of SPNE1 was lower than that of BON1 and QGP1. SPNE1 also was characterized with cancer stemness because of the higher proportion of CD44 + and CD117 + subpopulations relative to BON1, whereas it was similar to that of QGP1. Interestingly, SPNE1 highly expressed somatostatin receptors (SSTR2 and SSTR5) and angiogenic factors (VEGF1). SPNE1 had sensitive response to the four clinical treatments including tyrosine kinase inhibitor (TKI), mTOR inhibitors, somatostatin analogs (SSA), chemotherapy, which was similar to the BON1 and QGP1. Subcutaneous transplantations of SPNE1 also present the tumorigenicity, and neuroendocrine marker expression of xenograft tumors resembled the original human NET tissue. Then, we found a total of 8 common mutation in BON1, QGP1 and SPNE1 included CROCC, FAM135A, GPATCH4, CTBP2, FBXL14, HERC2, HYDIN, and PABPC3 using whole-exome sequencing (WES), and more neuroendocrine-related functional processes were enriched based on the private mutation genes in SPNE1, such as neuron migration, insulin secretion, and neuron to neuron synapse. In brief, SPNE1 could be used as a relevant model to study pancreatic NET biology and to develop novel treatment options.