Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “genome assembly validation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9Linked to original sources

Loss of maternal PADI6 disrupts DNA methylation and genomic imprinting maintenance in late preimplantation mouse embryos.

BACKGROUND: The maternal-effect protein PADI6, which is part of the subcortical maternal complex, is involved in proper spindle assembly, organelle distribution, ribosome storage, and cytoplasmic lattice organization in mouse oocytes. In humans, variants of PADI6 are associated with female infertility and multilocus imprinting disturbance in offspring. Recently, it was demonstrated that PADI6 plays a role in the storage and cytoplasmic localization of epigenetic factors, including UHRF1 and DNMT1. Moreover, maternal PADI6 depletion leads to defective epigenetic reprogramming and zygotic genome activation but not to an imprinting defect in two-cell mouse embryos. These findings raise the possibility that imprinting disturbances arise later in development. RESULTS: By employing combined single-blastocyst RNA-seq/BS-seq and immunostaining validation in the embryos derived from Padi6P620A-mutant oocytes, we investigated the role of Padi6 in late preimplantation development. We demonstrated that embryos that overcame the two-cell stage block had a dramatic reduction in UHRF1 and DNMT1 protein levels, a decrease in H3K9me3, and whole-genome hypomethylation, including most imprinted loci and repetitive elements, at the blastocyst stage. Furthermore, these maternal mutant embryos showed deregulation of inner cell mass markers and defective blastocyst implantation, but no effect on trophoblast differentiation. CONCLUSION: Our results demonstrate that maternal PADI6 is a key regulator of the stability of epigenetic factors required to maintain repressive marks in late preimplantation mouse embryos. Its deficiency results in genomic imprinting defects that closely resemble those found in human patients and provide a mechanistic explanation for MLID caused by maternal PADI6 variants. Furthermore, the impairment of blastocyst implantation capacity, likely due to dysregulation of inner cell mass differentiation, provides new mechanistic insights into the control of female fertility and embryo development exerted by PADI6.

DNA Methylation↗

Nucleocapsid protein of SARS coronavirus tightly binds to human cyclophilin A.

Severe acute respiratory syndrome coronavirus (SARS-CoV) is responsible for SARS infection. Nucleocapsid protein (NP) of SARS-CoV (SARS_NP) functions in enveloping the entire genomic RNA and interacts with viron structural proteins, thus playing important roles in the process of virus particle assembly and release. Protein-protein interaction analysis using bioinformatics tools indicated that SARS_NP may bind to human cyclophilin A (hCypA), and surface plasmon resonance (SPR) technology revealed this binding with the equilibrium dissociation constant ranging from 6 to 160nM. The probable binding sites of these two proteins were detected by modeling the three-dimensional structure of the SARS_NP-hCypA complex, from which the important interaction residue pairs between the proteins were deduced. Mutagenesis experiments were carried out for validating the binding model, whose correctness was assessed by the observed effects on the binding affinities between the proteins. The reliability of the binding sites derived by the molecular modeling was confirmed by the fact that the computationally predicted values of the relative free energies of the binding for SARS_NP (or hCypA) mutants to the wild-type hCypA (or SARS_NP) are in good agreement with the data determined by SPR. Such presently observed SARS_NP-hCypA interaction model might provide a new hint for facilitating the understanding of another possible SARS-CoV infection pathway against human cell.

Amino Acid Sequence↗

A comparative analysis of HGSC and Celera human genome assemblies and gene sets.

MOTIVATION: Since the simultaneous publication of the human genome assembly by the International Human Genome Sequencing Consortium (HGSC) and Celera Genomics, several comparisons have been made of various aspects of these two assemblies. In this work, we set out to provide a more comprehensive comparative analysis of the two assemblies and their associated gene sets. RESULTS: The local sequence content for both draft genome assemblies has been similar since the early releases, however it took a year for the quality of the Celera assembly to approach that of HGSC, suggesting an advantage of HGSC's hierarchical shotgun (HS) sequencing strategy over Celera's whole genome shotgun (WGS) approach. While similar numbers of ab initio predicted genes can be derived from both assemblies, Celera's Otto approach consistently generated larger, more varied gene sets than the Ensembl gene build system. The presence of a non-overlapping gene set has persisted with successive data releases from both groups. Since most of the unique genes from either genome assembly could be mapped back to the other assembly, we conclude that the gene set discrepancies do not reflect differences in local sequence content but rather in the assemblies and especially the different gene-prediction methodologies.

Databases, Protein↗

Cerebellar gene expression profiling and eQTL analysis in inbred mouse strains selected for ethanol sensitivity.

BACKGROUND: Inbred Long-Sleep (ILS) and Inbred Short-Sleep (ISS) mice exhibit striking differences in a number of alcohol and drug related behaviors. This study examined the expression levels of more than 39,000 transcripts in these strains in the cerebellum, a major target of ethanol's actions in the CNS, to find differentially expressed (DE) candidate genes for these phenotypes. METHODS: Genes that were differentially expressed between the strains were identified using oligonucleotide arrays as well as complimentary DNA arrays. Sequence alignment was used to locate DE genes in the mouse genome assembly. In silico expression QTL (eQTL) mapping was used to identify chromosomal regions likely to control the transcription level of DE genes, and the EASE program identified overrepresented functional themes. The genomic region immediately upstream of the cyclase associated protein homolog 1 (Cap1) gene was directly sequenced from PCR products. RESULTS: Nearly 300 genes were identified as differentially expressed between the cerebella of ILS and ISS. These genes and their corresponding eQTLs map to genomic regions linked to several phenotypes that differ between the ILS and ISS strains, including ethanol preference and cocaine-induced locomotor activation on Chromosomes 4 and 7 respectively. Eight genes were cross-platform validated, four of which are more highly expressed in ILS cerebellum. Three SNPs, one of which disrupts a predicted Sp1 binding site, were found in the upstream region of Cap1, a strong candidate for influencing ethanol phenotypes. CONCLUSIONS: Many of these DE genes are candidates to influence ethanol and drug regulated phenotypes because they either map to ethanol related QTLs in the genome or are linked to them through eQTL mapping. Genes involved in calcium ion binding and transcriptional regulation are overrepresented and therefore these gene classes may influence ethanol behaviors in mice and humans.

Adenylyl Cyclases↗

High-throughput assays using a luciferase-expressing replicon, virus-like particles, and full-length virus for West Nile virus drug discovery.

Many flaviviruses cause significant human disease worldwide. The development of flavivirus chemotherapy requires reliable high-throughput screening (HTS) assays. Although genetic systems have been developed for many flaviviruses, their usage in antiviral HTS assays has not been well explored. Here we compare three cell-based HTS assays for West Nile virus (WNV) drug discovery: (i) an assay that uses a cell line harboring a persistently replicating subgenomic replicon (containing a deletion of viral structural genes), (ii) an assay that uses packaged virus-like particles containing replicon RNA, and (iii) an assay that uses a full-length reporting virus. A Renilla luciferase gene was engineered into the replicon or into the full-length viral genome to monitor viral replication. Potential inhibitors could be identified through suppression of luciferase signals upon compound incubation. The antiviral assays were optimized in a 96-well format, validated with known WNV inhibitors, and proved useful in identifying a new inhibitor(s) through HTS of a compound library. In addition, because each assay encompasses multiple but discrete steps of the viral life cycle, the three systems could potentially be used to discriminate the mode of action of any inhibitor among viral entry (detected by assays ii and iii but not by assay i), replication (including viral translation and RNA synthesis; detected by assays i to iii), and virion assembly (detected by assay iii but not by assays i and ii). The approaches described in this study should be applicable to the development of cell-based assays for other flaviviruses.

Antiviral Agents↗

Machine learning prognostic model and drug survival analysis for lung adenocarcinoma in the context of radiotherapy.

BACKGROUND: Patients with lung adenocarcinoma (LUAD) receiving radiotherapy represent an important but underexplored clinical subgroup. These patients often undergo concomitant pharmacologic treatments, yet the prognostic impact and underlying determinants of such combined regimens remain poorly understood. OBJECTIVE: This retrospective observational study aimed to develop and validate a radiotherapy-specific machine learning prognostic model for LUAD and to compare survival across concomitant pharmacologic regimens. METHODS: In this retrospective observational study, using genomic and clinical data from TCGA, a radiotherapy-specific prognostic model for LUAD was developed and validated through ten machine learning algorithms. Survival analyses were conducted across distinct concomitant pharmacologic strategies, followed by functional enrichment to elucidate molecular mechanisms underlying differential outcomes. RESULTS: Demonstrating robust prognostic abilities, the model efficiently sorted patients into high- and low-risk categories. Both treatment type and risk score independently predicted overall survival, with significant interaction effects. Low-risk patients receiving targeted or combination therapy-mainly erlotinib, gefitinib, or bevacizumab-exhibited substantially improved survival compared with those receiving conventional chemotherapy. Enrichment of "Exogenous peptide presentation," "MHC class II assembly," "Peptide-MHC II assembly," and "Symbiotic interaction" pathways indicated immune modulation and host-tumor crosstalk as key mediators of treatment efficacy. CONCLUSION: This study establishes a radiotherapy-specific prognostic model for lung adenocarcinoma, demonstrating distinct molecular and therapeutic heterogeneity and highlighting the superior survival benefit of targeted combination therapy in low-risk patients.

Humans↗

Inhibition of the atypical kinase WNK1 as a therapeutic strategy in TAL-related T-cell acute lymphoblastic leukemia.

Driver mutations in T-cell acute lymphoblastic leukemia (T-ALL) rarely affect druggable kinases. However, these kinases can be aberrantly activated or repressed as secondary oncogenic events. Thus, integrating unbiased phosphoproteomics with genomic approaches may offer novel opportunities for target discovery and therapeutic interventions. In our study, we identified WNK1 (with no lysine [K]) as a potential target in T-ALL by pairing a list of vulnerable kinases with data from a phosphoproteomic screen of T-ALL cell lines. We subsequently validated WNK1 by loss-of-function-based studies and tested WNK inhibitors in several in vitro and in vivo T-ALL models and clinical T-ALL samples. We showed that therapeutic WNK1 repression promotes polyploidy, resulting in cell proliferation arrest, and morphometric changes, such as incomplete cell division or chromosome segregation through altered mitotic spindle assembly and abscission defects. Furthermore, we found that WNK1 is overexpressed in the TAL1/2-related subgroup, but not in normal thymus or lymph nodes, suggesting a potential translational area for clinical exploitation in poor-prognosis T-ALL carrying PTEN mutations and del(6q). Our work also reports a functional contribution of WNK1 in the leukemia establishment and progression. Structurally WNK1 is an atypical serine/threonine kinase that diverges from canonical kinases by lacking the conserved lysine in subdomain II, instead featuring a cysteine in subdomain I, which is critical for adenosine triphosphate (ATP) binding. This unusual structural configuration creates a distinct ATP-binding pocket with limited sequence similarity to conventional kinases, offering a unique opportunity to develop highly selective small molecules. Targeting this atypical ATP domain could thus provide a therapeutic advantage and broaden the treatment landscape for T-ALL.

WNK Lysine-Deficient Protein Kinase 1↗

Shotgun optical maps of the whole Escherichia coli O157:H7 genome.

We have constructed NheI and XhoI optical maps of Escherichia coli O157:H7 solely from genomic DNA molecules to provide a uniquely valuable scaffold for contig closure and sequence validation. E. coli O157:H7 is a common pathogen found in contaminated food and water. Our approach obviated the need for the analysis of clones, PCR products, and hybridizations, because maps were constructed from ensembles of single DNA molecules. Shotgun sequencing of bacterial genomes remains labor-intensive, despite advances in sequencing technology. This is partly due to manual intervention required during the last stages of finishing. The applicability of optical mapping to this problem was enhanced by advances in machine vision techniques that improved mapping throughput and created a path to full automation of mapping. Comparisons were made between maps and sequence data that characterized sequence gaps and guided nascent assemblies.

Contig Mapping↗

A genomic code for nucleosome positioning.

Eukaryotic genomes are packaged into nucleosome particles that occlude the DNA from interacting with most DNA binding proteins. Nucleosomes have higher affinity for particular DNA sequences, reflecting the ability of the sequence to bend sharply, as required by the nucleosome structure. However, it is not known whether these sequence preferences have a significant influence on nucleosome position in vivo, and thus regulate the access of other proteins to DNA. Here we isolated nucleosome-bound sequences at high resolution from yeast and used these sequences in a new computational approach to construct and validate experimentally a nucleosome-DNA interaction model, and to predict the genome-wide organization of nucleosomes. Our results demonstrate that genomes encode an intrinsic nucleosome organization and that this intrinsic organization can explain approximately 50% of the in vivo nucleosome positions. This nucleosome positioning code may facilitate specific chromosome functions including transcription factor binding, transcription initiation, and even remodelling of the nucleosomes themselves.

Base Sequence↗

Phenotypic and phylogenomic characterization of Lactococcus garvieae isolates from rainbow trout (Oncorhynchus mykiss) in Türkiye.

Lactococcosis is an important bacterial disease of farmed fish and causes substantial economic losses in rainbow trout (Oncorhynchus mykiss) aquaculture. In this study, Lactococcus garvieae isolates recovered from rainbow trout farms in Türkiye were characterized using phenotypic, molecular, and phylogenomic methods. Among 32 presumptive Lactococcus isolates recovered from 127 dead rainbow trout, four were confirmed as L. garvieae and exhibited identical biochemical characteristics, Pulsed Field Gel Electrophoresis (PFGE) profiles, and broad growth tolerance across different pH, salinity, and temperature conditions. All isolates were presumptively classified as resistant to ciprofloxacin and florfenicol, while remaining susceptible to tetracycline and penicillin. Based on the AMR profiles, strain LG2, which exhibited the most susceptible antimicrobial profile among the isolates, was selected for whole-genome sequencing (WGS). WGS of the representative isolate LG2 generated a single 2,214,687-bp chromosomal contig with 38.5% GC content and 99.0% BUSCO completeness. In silico PCR assigned LG2 to serotype I, and the genome contained an intact capsule-associated cps/kps locus. The chromosomal lsa(D) determinant and an mdt(A)-like efflux-associated gene were detected, whereas no plasmid replicons or acquired quinolone or florfenicol resistance genes were identified, indicating discordance between the phenotypic and genomic AMR results. Taxonomic verification of 236 publicly available Lactococcus assemblies yielded 41 verified public L. garvieae genomes, which, together with LG2, formed a 42-genome within-species dataset. LG2 was most closely related to the Turkish isolate OS-37, sharing 99.96% ANI and differing by three core SNPs; both belonged to ST109, whereas the other Turkish isolates belonged to ST139. cgMLST identified a conserved genomic backbone, while pan-genome analysis identified 5,655 gene clusters and an open pan-genome characterized by a large cloud-gene fraction. These findings demonstrate the importance of species verification in Lactococcus population genomics and reveal substantial accessory-genome diversity within L. garvieae. The genomic features of LG2 provide a basis for future pathogenicity and immunogenicity studies, although experimental validation is required. Overall, these findings highlight the importance of local genomic surveillance for understanding L. garvieae population structure and provide a genomic framework for future region-specific vaccine research.

Animals↗

Kynurenine metabolism-related gene signature for prognostic stratification in hepatocellular carcinoma.

BACKGROUND: Hepatocellular carcinoma (HCC) remains a major global health burden with high mortality rates and limited therapeutic options. The identification of reliable biomarkers for early diagnosis and prognosis prediction is urgently needed. Kynurenine metabolism, a critical pathway in immune regulation and tumor progression, has been implicated in various cancers. However, its prognostic value in HCC has not been fully elucidated. This study aimed to develop a prognostic risk model based on kynurenine metabolism-related genes (KMRGs) for HCC patients. METHODS: Transcriptomic and clinical data of HCC patients were retrieved from The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) databases. A prognostic risk model was established using least absolute shrinkage and selection operator (LASSO) and Cox regression analyses. Survival analysis and functional enrichment analysis were conducted to validate the predictive performance of the model and to investigate the underlying mechanisms. ALDH8A1 was ultimately identified as a target gene based on survival analysis, and its impact on tumor cell migration was assessed using the HCC cell line. RESULTS: A prognostic model based on seven KMRGs was established. The high-risk group exhibited significantly worse overall survival compared to the low-risk group. Functional enrichment analysis in high-risk patients highlighted significant enrichment in core biological processes, including spliceosome assembly and ribonucleoprotein complex biogenesis. Furthermore, a nomogram integrating the risk score and clinical pathological features was developed, demonstrating moderate predictive performance for HCC prognosis. CONCLUSIONS: This study successfully constructed a prognostic risk model based on seven KMRGs, providing a valuable tool for predicting clinical outcomes in HCC patients. These findings highlight the potential role of kynurenine metabolism in HCC progression and offer new insights for future therapeutic strategies.

ALDH8A1↗

Identification of transposable elements using multiple alignments of related genomes.

Accurate genome-wide cataloging of transposable elements (TEs) will facilitate our understanding of mobile DNA evolution, expose the genomic effects of TEs on the host genome, and improve the quality of assembled genomes. Using the availability of several nearly complete Drosophila genomes and developments in whole genome alignment methods, we introduce a large-scale comparative method for identifying repetitive mobile DNA regions. These regions are highly enriched for transposable elements. Our method has two main features distinguishing it from other repeat-finding methods. First, rather than relying on sequence similarity to determine the location of repeats, the genomic artifacts of the transposition mechanism itself are systematically tracked in the context of multiple alignments. Second, we can derive bounds on the age of each repeat instance based on the phylogenetic species tree. We report results obtained using both complete and draft sequences of four closely related Drosophila genomes and validate our results with manually curated TE annotations in the Drosophila melanogaster euchromatin. We show the utility of our findings in exploring both transposable elements and their host genomes: In the study of TEs, we offer predictions for novel families, annotate new insertions of known families, and show data that support the hypothesis that all known TE families in D. melanogaster were recently active; in the study of the host, we show how our findings can be used to determine shifts in the eu-heterochromatin junction in the pericentric chromosome regions.

Animals↗

Structure modeling of all identified G protein-coupled receptors in the human genome.

G protein-coupled receptors (GPCRs), encoded by about 5% of human genes, comprise the largest family of integral membrane proteins and act as cell surface receptors responsible for the transduction of endogenous signal into a cellular response. Although tertiary structural information is crucial for function annotation and drug design, there are few experimentally determined GPCR structures. To address this issue, we employ the recently developed threading assembly refinement (TASSER) method to generate structure predictions for all 907 putative GPCRs in the human genome. Unlike traditional homology modeling approaches, TASSER modeling does not require solved homologous template structures; moreover, it often refines the structures closer to native. These features are essential for the comprehensive modeling of all human GPCRs when close homologous templates are absent. Based on a benchmarked confidence score, approximately 820 predicted models should have the correct folds. The majority of GPCR models share the characteristic seven-transmembrane helix topology, but 45 ORFs are predicted to have different structures. This is due to GPCR fragments that are predominantly from extracellular or intracellular domains as well as database annotation errors. Our preliminary validation includes the automated modeling of bovine rhodopsin, the only solved GPCR in the Protein Data Bank. With homologous templates excluded, the final model built by TASSER has a global C(alpha) root-mean-squared deviation from native of 4.6 angstroms, with a root-mean-squared deviation in the transmembrane helix region of 2.1 angstroms. Models of several representative GPCRs are compared with mutagenesis and affinity labeling data, and consistent agreement is demonstrated. Structure clustering of the predicted models shows that GPCRs with similar structures tend to belong to a similar functional class even when their sequences are diverse. These results demonstrate the usefulness and robustness of the in silico models for GPCR functional analysis. All predicted GPCR models are freely available for noncommercial users on our Web site (http://www.bioinformatics.buffalo.edu/GPCR).

Algorithms↗

Pharmacogenomics of phenolic antioxidant butylated hydroxyanisole (BHA) in the small intestine and liver of Nrf2 knockout and C57BL/6J mice.

PURPOSE: The objective of this study was to investigate the pharmacogenomics and the spatial regulation of global gene expression profiles elicited by cancer chemopreventive agent butylated hydroxyanisole (BHA) in mouse small intestine and liver as well as to identify BHA-modulated nuclear factor-E2-related factor 2 (Nrf2)-dependent genes. METHODS: C57BL/6J (+/+; wildtype) and C57BL/6J/Nrf2(-/-; knockout) mice were administered a single 200 mg/kg oral dose of BHA or only vehicle. Both small intestine and liver were collected at 3 h after treatment and total RNA was extracted. Gene expression profiles were analyzed using 45,000 Affymetrix mouse genome 430 2.0 array and GeneSpring 7.2 software. Microarray results were validated by quantitative real-time reverse transcription-PCR analyses. RESULTS: Clusters of genes that were either induced or suppressed more than two fold by BHA treatment compared with vehicle in C57BL/6J/Nrf2(-/-; knockout) and C57BL/6J Nrf2 (+/+; wildtype) mice genotypes were identified. Amongst these, in small intestine and liver, 1,490 and 493 genes respectively were identified as Nrf2-dependent and upregulated, and 1,090 and 824 genes respectively as Nrf2-dependent and downregulated. Based on their biological functions, these genes can be categorized into ubiquitination/proteolysis, apoptosis/cell cycle, electron transport, detoxification, cell growth/differentiation, transcription factors/interacting partners, kinases and phosphatases, transport, biosynthesis/metabolism, RNA/protein processing and nuclear assembly, and DNA replication genes. Phase II detoxification/antioxidant genes as well as novel molecular target genes, including putative interacting partners of Nrf2 such as nuclear corepressors and coactivators, were also identified as Nrf2-dependent genes. CONCLUSIONS: The identification of BHA-regulated and Nrf2-dependent genes not only provides potential novel insights into the gestalt biological effects of BHA on the pharmacogenomics and spatial regulation of global gene expression profiles in cancer chemoprevention, but also points to the pivotal role of Nrf2 in these biological processes.

Animals↗

Toxicogenomics of endoplasmic reticulum stress inducer tunicamycin in the small intestine and liver of Nrf2 knockout and C57BL/6J mice.

This objective of this study was to investigate the toxicogenomics and the spatial regulation of global gene expression profiles elicited by endoplasmic reticulum (ER) stress inducer tunicamycin (TM) in mouse small intestine and liver as well as to identify TM-modulated nuclear factor-E2-related factor 2 (Nrf2)-dependent genes. Gene expression profiles were analyzed using 45,000 Affymetrix mouse genome 430 2.0 array and GeneSpring 7.2 software. Microarray results were validated by quantitative real-time reverse transcription-PCR analyses. Clusters of genes that were either induced or suppressed more than two-fold by TM treatment compared with vehicle in C57BL/6J/Nrf2 (-/-; knockout) and C57BL/6J Nrf2 (+/+; wildtype) mice genotypes were identified. Amongst these, in small intestine and liver, 1291 and 750 genes, respectively, were identified as Nrf2-dependent and upregulated, and 1370 and 943 genes, respectively, as Nrf2-dependent and downregulated. Based on their biological functions, these genes can be categorized into molecular chaperones and heat shock proteins, ubiquitination/proteolysis, apoptosis/cell cycle, electron transport, detoxification, cell growth/differentiation, signaling molecules/interacting partners, kinases and phosphatases, transport, biosynthesis/metabolism, nuclear assembly and processing, and genes related to calcium and glucose homeostasis. Phase II detoxification/antioxidant genes as well as putative interacting partners of Nrf2 such as nuclear corepressors and coactivators, were also identified as Nrf2-dependent genes. The identification of TM-regulated and Nrf2-dependent genes in the unfolded protein response to ER stress not only provides potential novel insights into the gestalt biological effects of TM on the toxicogenomics and spatial regulation of global gene expression profiles in cancer pharmacology and toxicology, but also points to the pivotal role of Nrf2 in these biological processes.

Animals↗

ProtEST: protein multiple sequence alignments from expressed sequence tags.

MOTIVATION: An automatic sequence searching method (ProtEST) is described which constructs multiple protein sequence alignments from protein sequences and translated expressed sequence tags (ESTs). ProtEST is more effective than a simple TBLASTN search of the query against the EST database, as the sequences are automatically clustered, assembled, made non-redundant, checked for sequence errors, translated into protein and then aligned and displayed. RESULTS: A ProtEST search found a non-redundant, translated, error- and length-corrected EST sequence for > 58% of sequences when single sequences from 1407 Pfam-A seed alignments were used as the probe. The average family size of the resulting alignments of translated EST sequences contained > 10 sequences. In a cross-validated test of protein secondary structure prediction, alignments from the new procedure led to an improvement of 3.4% average Q3 prediction accuracy over single sequences. AVAILABILITY: The ProtEST method is available as an Internet World Wide Web service http://barton.ebi.ac.uk/servers/protest.html+ ++ The Wise2 package for protein and genomic comparisons and the ProtESTWise script can be found at http://www.sanger.ac.uk/Software/Wise2 CONTACT: geoff@ebi.ac.uk

Amino Acid Sequence↗

Chromosome-Level Genome Assembly of Eden's Whale Clarifies the Taxonomy and Speciation of Bryde's Whale Complex.

Eden's whale (Balaenoptera edeni), a poorly understood baleen cetacean, has long been shrouded in taxonomic ambiguity due to limited genomic resources, obscuring its distinction from closely related species and its position within the cetacean Tree of Life. In this paper, we present a high-quality chromosomal-level genome of B. edeni and conduct comparative genomic analyses to address long-standing taxonomic confusion and elucidate speciation of balaenopterids. Our phylogenomic analysis and demographic reconstruction reveal that B. edeni is a distinct sister to Bryde's whale (Balaenoptera brydei), sharing a common ancestor that diverged approximately 7.84 million years ago during the late Miocene. Their genetic divergence exceeds typical intraspecific variation in whales, supporting the reinstatement of B. brydei as a valid species. Chromosomal syntenic analyses suggest that macro-fragment inversions contributed to speciation in balaenopterid whales and uncover unexpected large-scale complex genome rearrangements in Bryde's whale, offering novel insights into cetacean genome evolution. Functional enrichment analysis of inverted regions between B. edeni and Balaenoptera musculus indicates their predominant association with metabolism and biosynthesis, as well as responses to various substances, stress, and stimuli. These genomic resources for B. edeni not only lay a critical foundation for comparative genetic and evolutionary research of cetaceans but also advance our understanding of the taxonomy and evolutionary dynamics of the Bryde's whale complex, with broader implications for baleen whale conservation and biodiversity.

Animals↗

Identification of human chromosome 22 transcribed sequences with ORF expressed sequence tags.

Transcribed sequences in the human genome can be identified with confidence only by alignment with sequences derived from cDNAs synthesized from naturally occurring mRNAs. We constructed a set of 250,000 cDNAs that represent partial expressed gene sequences and that are biased toward the central coding regions of the resulting transcripts. They are termed ORF expressed sequence tags (ORESTES). The 250,000 ORESTES were assembled into 81,429 contigs. Of these, 1, 181 (1.45%) were found to match sequences in chromosome 22 with at least one ORESTES contig for 162 (65.6%) of the 247 known genes, for 67 (44.6%) of the 150 related genes, and for 45 of the 148 (30.4%) EST-predicted genes on this chromosome. Using a set of stringent criteria to validate our sequences, we identified a further 219 previously unannotated transcribed sequences on chromosome 22. Of these, 171 were in fact also defined by EST or full length cDNA sequences available in GenBank but not utilized in the initial annotation of the first human chromosome sequence. Thus despite representing less than 15% of all expressed human sequences in the public databases at the time of the present analysis, ORESTES sequences defined 48 transcribed sequences on chromosome 22 not defined by other sequences. All of the transcribed sequences defined by ORESTES coincided with DNA regions predicted as encoding exons by genscan. (http://genes.mit.edu/GENSCAN.html).

Chromosomes, Human, Pair 22↗