Search

Genomic rearrangements, encompassing mutational changes in the genome such as insertions, deletions or inversions, are essential for genetic diversity. These rearrangements are typically orchestrated by enzymes that are involved in fundamental DNA repair processes, such as homologous recombination, or in the transposition of foreign genetic material by viruses and mobile genetic elements 1,2 . Here we report that IS110 insertion sequences, a family of minimal and autonomous mobile genetic elements, express a structured non-coding RNA that binds specifically to their encoded recombinase. This bridge RNA contains two internal loops encoding nucleotide stretches that base-pair with the target DNA and the donor DNA, which is the IS110 element itself. We demonstrate that the target-binding and donor-binding loops can be independently reprogrammed to direct sequence-specific recombination between two DNA molecules. This modularity enables the insertion of DNA into genomic target sites, as well as programmable DNA excision and inversion. The IS110 bridge recombination system expands the diversity of nucleic-acid-guided systems beyond CRISPR and RNA interference, offering a unified mechanism for the three fundamental DNA rearrangements-insertion, excision and inversion-that are required for genome design., (© 2024. The Author(s).)

The human mannose-binding lectin (MBL) gene encodes a polymorphic protein that plays a crucial role in the innate immune response. Human MBL deficiency is associated with immunodeficiencies, and its variants have been linked to autoimmune and infectious diseases. Despite this significance, gene studies concerning MBL sequencing are uncommon in Malaysia. Therefore, we aimed to preliminary described the human MBL sequencing dataset based on the Kelantan population. Blood samples were collected from 30 unrelated individuals and underwent DNA extraction, genotyping, and sequencing. The sequencing data generated 886 bp, which were deposited in GenBank (ON619541-ON619546). Allelic variants were identified and translated into six MBL haplotypes: HYPA, HYPB, LYPB, LXPB, HXPA, and LXPA. An evolutionary tree was constructed using the haplotype sequences. These findings contribute to the expansion of MBL information within the country, providing a valuable baseline for future research exploring the association between the gene and targeted diseases., (© 2024. The Author(s).)

Long-read genome sequencing (lrGS) is a promising method in genetic diagnostics. Here we investigate the potential of lrGS to detect a disease-associated chromosomal translocation between 17p13 and the 19 centromere. We constructed two sets of phased and non-phased de novo assemblies; (i) based on lrGS only and (ii) hybrid assemblies combining lrGS with optical mapping using lrGS reads with a median coverage of 34X. Variant calling detected both structural variants (SVs) and small variants and the accuracy of the small variant calling was compared with those called with short-read genome sequencing (srGS). The de novo and hybrid assemblies had high quality and contiguity with N50 of 62.85 Mb, enabling a near telomere to telomere assembly with less than a 100 contigs per haplotype. Notably, we successfully identified the centromeric breakpoint of the translocation. A concordance of 92% was observed when comparing small variant calling between srGS and lrGS. In summary, our findings underscore the remarkable potential of lrGS as a comprehensive and accurate solution for the analysis of SVs and small variants. Thus, lrGS could replace a large battery of genetic tests that were used for the diagnosis of a single symptomatic translocation carrier, highlighting the potential of lrGS in the realm of digital karyotyping., (© 2024. The Author(s).)

Major depressive disorder (MDD) is a complex and potentially debilitating illness whose etiology and pathology remains unclear. Non-coding RNAs have been implicated in MDD, where they display differential expression in the brain and the periphery. In this study, we quantified small nucleolar RNA (snoRNA) expression by small RNA sequencing in the lateral habenula (LHb) of individuals with MDD (n = 15) and psychiatrically-healthy controls (n = 15). We uncovered five snoRNAs that exhibited differential expression between MDD and controls (FDR < 0.01). Specifically, SNORA69 showed increased expression in MDD and was technically validated via RT-qPCR. We further investigated the expression of Snora69 in the LHb and peripheral blood of an unpredicted chronic mild stress (UCMS) mouse model of depression. Snora69 was specifically up-regulated in mice that underwent the UCMS paradigm. SNORA69 is known to guide pseudouridylation onto 5.8S and 18S rRNAs. We quantified the relative abundance of pseudouridines on 5.8S and 18S rRNA in human post-mortem LHb samples and found increased abundance of pseudouridines in the MDD group. Overall, our findings indicate the importance of brain snoRNAs in the pathology of MDD. Future studies characterizing SNORA69's role in MDD pathology is warranted., (© 2024. The Author(s).)

Athetis lepigone is an emerging highly polyphagous insect pest reported to cause crop damage in several European and Asian countries. However, our understanding of its genetic adaptation mechanisms has been limited due to lack of high-quality genetic resources. In this study, we present a chromosomal-level genome of A. lepigone, representing the first species in the genus of Athetis. We employed PacBio long-read sequencing and Hi-C technologies to generate 612.49 Mb genome assembly which contains 42.43% repeat sequences with a scaffold N50 of 20.9 Mb. The contigs were successfully clustered into 31 chromosomal-size scaffolds with 37% GC content. BUSCO assessment revealed a genome completeness of 97.4% with 96.3 identified as core Arthropoda single copy orthologs. Among the 17,322 genes that were predicted, 15,965 genes were functionally annotated, representing a coverage of 92.17%. Furthermore, we revealed 106 P450, 37 GST, 27 UGT, and 74 COE gene families in the genome of A. lepigone. This genome provides a significant and invaluable genomic resource for further research across the entire genus of Athetis., (© 2024. The Author(s).)

Abies koreana E.H.Wilson is an endangered evergreen coniferous tree that is native to high altitudes in South Korea and susceptible to the effects of climate change. Hybridization and reticulate evolution have been reported in the genus; therefore, multigene datasets from nuclear and cytoplasmic genomes are needed to better understand its evolutionary history. Using the Illumina NovaSeq 6000 and Oxford Nanopore Technologies (ONT) PromethION platforms, we generated complete mitochondrial (1,174,803 bp) and plastid (121,341 bp) genomes from A. koreana. The mitochondrial genome is highly dynamic, transitioning from cis- to trans-splicing and breaking conserved gene clusters. In the plastome, the ONT reads revealed two structural conformations of A. koreana. The short inverted repeats (1186 bp) of the A. koreana plastome are associated with different structural types. Transcriptomic sequencing revealed 1356 sites of C-to-U RNA editing in the 41 mitochondrial genes. Using A. koreana as a reference, we additionally produced nuclear and organelle genomic sequences from eight Abies species and generated multiple datasets for maximum likelihood and network analyses. Three sections (Balsamea, Momi, and Pseudopicea) were well grouped in the nuclear phylogeny, but the phylogenomic relationships showed conflicting signals in the mitochondrial and plastid genomes, indicating a complicated evolutionary history that may have included introgressive hybridization. The obtained data illustrate that phylogenomic analyses based on sequences from differently inherited organelle genomes have resulted in conflicting trees. Organelle capture, organelle genome recombination, and incomplete lineage sorting in an ancestral heteroplasmic individual can contribute to phylogenomic discordance. We provide strong support for the relationships within Abies and new insights into the phylogenomic complexity of this genus., (© 2024. The Author(s).)

Natural killer (NK) cells play essential roles in the tumor development, diagnosis, and prognosis of tumors. In this study, we aimed to establish a reliable signature based on marker genes in NK cells, thus providing a new perspective for assessing immunotherapy and the prognosis of patients with gastric cancer (GC). We analyzed a total of 1560 samples retrieved from the public database. We performed a comprehensive analysis of single-cell RNA-sequencing (scRNA-seq) data of gastric cancer and identified 377 marker genes for NK cells. By performing Cox regression analysis, we established a 12-gene NK cell-associated signature (NKCAS) for the Cancer Genome Atlas (TCGA) cohort, that assigned GC patients into a low-risk group (LRG) or a high-risk group (HRG). In the TCGA cohort, the areas under curve (AUC) value were 0.73, 0.81, and 0.80 at 1, 3, and 5 years. External validation of the predictive ability for the signature was then validated in the Gene Expression Omnibus (GEO) cohorts (GSE84437). The expression levels of signature genes were measured and validated in GC cell lines by real-time PCR. Moreover, NKCAS was identified as an independent prognostic factor by multivariate analysis. We combined this with a variety of clinicopathological characteristics (age, M stage, and tumor grade) to construct a nomogram to predict the survival outcomes of patients. Moreover, the LRG showed higher immune cell infiltration, especially CD8+ T cells and NK cells. The risk score was negatively associated with inflammatory activities. Importantly, analysis of the independent immunotherapy cohort showed that the LRG had a better prognosis and immunotherapy response when compared with the HRG. The identification of NK cell marker genes in this study suggests potential therapeutic targets. Additionally, the developed predictive signatures and nomograms may aid in the clinical management of GC., (© 2024. The Author(s).)

Diplopoda is one of the most diverse and important groups of soil arthropods, but little research has been done on their phylogenetic relationship and evolution. Here, we sequenced and annotated the complete mitochondrial genomes of Spirobolus grahami. The total mitogenome of S. grahami was typical circular, double-stranded molecules, with 14,875 bp in length, including 13 protein-coding genes, 22 tRNAs, two rRNAs, and one control region. Base composition analysis suggested that the mitochondrial sequences were biased toward A and T, with A + T content of 58.68%. The mitogenomes of S. grahami exhibited negative AT and positive GC skews. Most of the 13 PCGs had ATN as the start codon, except COX1 start with CGA, and most PCGs ended with the T stop codon. The dN/dS values for most PCGs were lower than 1, suggesting that purifying selection was likely the main driver of mitochondrial PCG evolution. Phylogenetic analyses based on 13 PCGs using BI and ML methods support the classification of genus Spirobolus and Tropostreptus. Glomeridesmus spelaeus is distantly related to the other Diplopoda species., (© 2024. The Author(s).)

Balantidium ctenopharyngodoni is identified as the sole ciliate species that exclusively resides within the hindgut of grass carp with high prevalence and intensity. In this study, the successful cultivation of B. ctenopharyngodoni enabled us to collect enough cells for genome sequencing. Consequently, we acquired a high-quality genome assembly spanning 68.66 Mb, encompassing a total of 22,334 nanochromosomes. Furthermore, we predicted 29,348 protein-coding genes, and 95.5% of them was supported by the RNA-seq data. The trend of GC content in the subtelomeric regions of single-gene chromosomes was similar to other ciliates containing nanochromosomes. A large number of genes encoding carbohydrate-binding modules with affinities for starch and peptidoglycans was identified. The identification of mitochondrion-related organelles (MROs) within genome indicates its well-suited adaptation to the anaerobic conditions in the hindgut environment. In summary, our results will offer resources for understanding the genetic basis and molecular adaptations of balantidia to hindgut of herbivorous fish., (© 2024. The Author(s).)

The coffee leaf miner (Leucoptera coffeella) is one of the major pests of coffee crops in the neotropical regions, and causes major economic losses. Few molecular data are available to identify this pest and advances in the knowledge of the genome of L. coffeella will contribute to improving pest identification and also clarify taxonomy of this microlepidoptera. L. coffeella DNA was extracted and sequenced using PacBio HiFi technology. Here we report the complete L. coffeella circular mitochondrial genome (16,407 bp) assembled using Aladin software. We found a total of 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and an A + T rich-region and a D-loop. The L. coffeella mitochondrial gene organization is highly conserved with similarities to lepidopteran mitochondrial gene rearrangements (trnM-trnI-trnQ). We concatenated the 13 PCG to construct a phylogenetic tree and inferred the relationship between L. coffeella and other lepidopteran species. L. coffeella is found in the Lyonetiidae clade together with L. malifoliella and Lyonetia clerkella, both leaf miners. Interestingly, this clade is assigned in the Yponomeutoidea superfamily together with Gracillariidae, and both superfamilies displayed species with leaf-mining feeding habits., (© 2024. The Author(s).)

The Amplified Luminescent Proximity Homogenous Assay-linked Immunosorbent Assay (AlphaLISA) is known for detecting various protein targets; however, its ability to detect nucleic acid sequences is not well established. Here, the capabilities of the AlphaLISA technology were expanded to include direct detection of DNA (aka: oligo-Alpha) and was applied to the detection of Listeria monocytogenes. Parameters were defined that allowed the newly developed oligo-Alpha to differentiate L. monocytogenes from other Listeria species through the use of only a single nucleotide polymorphism within the 16S rDNA region. Investigations into the applicability of this assay with different matrices demonstrated its utility in both milk and juice. One remarkable feature of the oligo-Alpha is that greater sensitivity could be achieved through the use of multiple acceptor oligos compared to only a single acceptor oligo, even when only a single donor oligo was employed. Additional acceptor oligos were easily incorporated into the assay and a tenfold change in the detection limit was readily achieved, with detection limits of 250 attomole of target being recorded. In summary, replacement of antibodies with oligonucleotides allows us to take advantage of genotypic difference(s), which both expands its repertoire of biological markers and furthers its use as a diagnostic tool., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Transposon directed insertion-site sequencing (TraDIS), a variant of transposon insertion sequencing commonly known as Tn-Seq, is a high-throughput assay that defines essential bacterial genes across diverse growth conditions. However, the variability between laboratory environments often requires laborious, time-consuming modifications to its protocol. In this technical study, we aimed to refine the protocol by identifying key parameters that can impact the complexity of mutant libraries. Firstly, we discovered that adjusting electroporation parameters including transposome concentration, transposome assembly conditions, and cell densities can significantly improve the recovery of viable mutants for different Escherichia coli strains. Secondly, we found that post-electroporation conditions, such as recovery time and the use of different mediums for selecting mutants may also impact the complexity of viable mutants in the library. Finally, we developed a simplified sequencing library preparation workflow based on a Nextera-TruSeq hybrid design where ~ 80% of sequenced reads correspond to transposon-DNA junctions. The technical improvements presented in our study aim to streamline TraDIS protocols, making this powerful technique more accessible for a wider scientific audience., (© 2024. The Author(s).)

Many lines of evidence indicate that red jungle fowl (RJF) is the primary ancestor of domestic chickens. Although multiple versions of RJF (galgal2-galgal5 and GRCg6a) and commercial chickens (GRCg7b/w and Huxu) genomes have been assembled since 2004, no high-quality indigenous chicken genomes have been assembled, hampering the understanding of chicken domestication and evolution. To fill the gap, we sequenced the genomes of four indigenous chickens with distinct morphological traits in southwest China, using a combination of short, long and Hi-C reads. We assembled each genome (~1.0 Gb) into 42 chromosomes with chromosome N50 90.5-90.9 Mb, amongst the highest quality of chicken genome assemblies. To provide resources for gene annotation and functional analysis, we also sequenced transcriptomes of 10 tissues for each of the four chickens. Moreover, we corrected many mis-assemblies and assembled missing micro-chromosomes 29 and 34-39 for GRCg6a. Our assemblies, sequencing data and the correction of GRCg6a can be valuable resources for studying chicken domestication and evolution., (© 2024. The Author(s).)

Sharks have thrived in the oceans for 400 million years, experienced five extinctions and evolved into today's apex predators. However, enormous genome size, poor karyotyping and limited tissue sampling options are the bottlenecks in shark research. Sharks of the family Orectolobiformes act as model species in transcriptome research with exceptionally high reproductive fecundity, catch prominence and oviparity. The present study illustrates a de novo transcriptome for an adult grey bamboo shark, Chiloscyllium griseum (Chondrichthyes; Hemiscyllidae) using paired-end RNA sequencing. Around 150 million short Illumina reads were obtained from five different tissues and assembled using the Trinity assembler. 70,647 hits on Uniprot by BLASTX was obtained after the transcriptome annotation. The data generated serve as a basis for transcriptome-based population genetic studies and open up new avenues in the field of comparative transcriptomics and conservation biology., (© 2024. The Author(s).)

A major goal of evolutionary genetics is to understand the genetic and molecular mechanisms underlying adaptation. Previous work has established that changes in gene regulation may contribute to adaptive evolution, but most studies have focused on mRNA abundance and only a few studies have investigated the role of post-transcriptional processing. Here, we use a combination of exome sequences and short-read RNA-Seq data from wild house mice (Mus musculus domesticus) collected along a latitudinal transect in eastern North America to identify candidate genes for local adaptation through alternative splicing. First, we identified alternatively spliced transcripts that differ in frequency between mice from the northern-most and southern-most populations in this transect. We then identified the subset of these transcripts that exhibit clinal patterns of variation among all populations in the transect. Finally, we conducted association studies to identify cis-acting splicing quantitative trait loci (cis-sQTL), and we identified cis-sQTL that overlapped with previously ascertained targets of selection from genome scans. Together, these analyses identified a small set of alternatively spliced transcripts that may underlie environmental adaptation in house mice. Many of these genes have known phenotypes associated with body size, a trait that varies clinally in these populations. We observed no overlap between these genes and genes previously identified by changes in mRNA abundance, indicating that alternative splicing and changes in mRNA abundance may provide separate molecular mechanisms of adaptation., (© 2023. The Author(s), under exclusive licence to The Genetics Society.)

Numerous contiguous gene deletion syndromes causing neurodevelopmental disorders have previously been defined using cytogenetics for which only in the current genomic era the disease-causing genes have become elucidated. One such example is deletion at Xq22.2, previously associated with a neurodevelopmental disorder which has more recently been found to be caused by de novo loss-of-function variants in TCEAL1. So far, a single study reported six unrelated individuals with this monogenetic disorder, presenting with syndromic features including developmental delay especially affecting expressive speech, intellectual disability, autistic-like behaviors, hypotonia, gait abnormalities and mild facial dysmorphism, in addition to ocular, gastrointestinal, and immunologic abnormalities. Here we report on four previously undescribed individuals, including two adults, with de novo truncating variants in TCEAL1, identified through trio exome or genome sequencing, further delineating the phenotype of the TCEAL1-related disorder. Whereas overall we identify similar features compared to the original report, we also highlight features in our adult individuals including hyperphagia, obesity, and endocrine abnormalities including hyperinsulinemia, hyperandrogenemia, and polycystic ovarian syndrome. X chromosome inactivation and RNA-seq studies further provide functional insights in the molecular mechanisms. Together this report expands the phenotypic and molecular spectrum of the TCEAL1-related disorder which will be useful for counseling of newly identified individuals and their families., (© 2024. The Author(s).)

Denisova Cave in the Siberian Altai (Russia) is a key site for understanding the complex relationships between hominin groups that inhabited Eurasia in the Middle and Late Pleistocene epoch. DNA sequenced from human remains found at this site has revealed the presence of a hitherto unknown hominin group, the Denisovans.sup.1,2, and high-coverage genomes from both Neanderthal and Denisovan fossils provide evidence for admixture between these two populations.sup.3. Determining the age of these fossils is important if we are to understand the nature of hominin interaction, and aspects of their cultural and subsistence adaptations. Here we present 50 radiocarbon determinations from the late Middle and Upper Palaeolithic layers of the site. We also report three direct dates for hominin fragments and obtain a mitochondrial DNA sequence for one of them. We apply a Bayesian age modelling approach that combines chronometric (radiocarbon, uranium series and optical ages), stratigraphic and genetic data to calculate probabilistically the age of the human fossils at the site. Our modelled estimate for the age of the oldest Denisovan fossil suggests that this group was present at the site as early as 195,000 years ago (at 95.4% probability). All Neanderthal fossils--as well as Denisova 11, the daughter of a Neanderthal and a Denisovan.sup.4--date to between 80,000 and 140,000 years ago. The youngest Denisovan dates to 52,000-76,000 years ago. Direct radiocarbon dating of Upper Palaeolithic tooth pendants and bone points yielded the earliest evidence for the production of these artefacts in northern Eurasia, between 43,000 and 49,000 calibrated years before present (taken as ad 1950). On the basis of current archaeological evidence, it may be assumed that these artefacts are associated with the Denisovan population. It is not currently possible to determine whether anatomically modern humans were involved in their production, as modern-human fossil and genetic evidence of such antiquity has not yet been identified in the Altai region.Bayesian modelling of chronometric, stratigraphic and genetic data from Denisova Cave provides a chronological framework for understanding Neanderthal and Denisovan presence at the site, as well as interactions between these groups., Author(s): Katerina Douka [sup.1] [sup.2] , Viviane Slon [sup.3] , Zenobia Jacobs [sup.4] [sup.5] , Christopher Bronk Ramsey [sup.2] , Michael V. Shunkov [sup.6] [sup.7] , Anatoly P. Derevianko [sup.6] [...]

Chameleons are well-known lizards with unique morphology and physiology, but their sex determination has remained poorly studied. Madagascan chameleons of the genus Furcifer have cytogenetically distinct Z and W sex chromosomes and occasionally Z 1 Z 1 Z 2 Z 2 /Z 1 Z 2 W multiple neo-sex chromosomes. To identify the gene content of their sex chromosomes, we microdissected and sequenced the sex chromosomes of F. oustaleti (ZZ/ZW) and F. pardalis (Z 1 Z 1 Z 2 Z 2 /Z 1 Z 2 W). In addition, we sequenced the genomes of a male and a female of F. lateralis (ZZ/ZW) and F. pardalis and performed a comparative coverage analysis between the sexes. Despite the notable heteromorphy and distinctiveness in heterochromatin content, the Z and W sex chromosomes share approximately 90% of their gene content. This finding demonstrates poor correlation of the degree of differentiation of sex chromosomes at the cytogenetic and gene level. The test of homology based on the comparison of gene copy number variation revealed that female heterogamety with differentiated sex chromosomes remained stable in the genus Furcifer for at least 20 million years. These chameleons co-opted for the role of sex chromosomes the same genomic region as viviparous mammals, lacertids and geckos of the genus Paroedura, which makes these groups excellent model for studies of convergent and divergent evolution of sex chromosomes., (© 2024. The Author(s).)

There are four species in the Crossoptilon genus inhibiting at from very low to very high altitudes across China, and they are in varying levels of danger of extinction. To better understand the genetic basis of adaptation to high altitudes and genetic changes due to bottleneck, we assembled the genome (~1.02 Gb) of a white eared pheasant (WT) (Crossoptilon crossoptilon) inhibiting at high altitudes (3,000~7,000 m) in northwest of Yunnan province, China, using a combination of Illumina short reads, PacBio long reads and Hi-C reads, with a contig N50 of 19.63 Mb and only six gaps. To further provide resources for gene annotation as well as functional and population genetics analyses, we sequenced transcriptomes of 20 major tissues of the WT individual and re-sequenced another 10 WT individuals and a blue eared pheasant (Crossoptilon auritum) individual inhabiting at intermediate altitudes (1,500~3,000 m). Our assembled WT genome, transcriptome data, and DNA sequencing data can be valuable resources for studying the biology, evolution and developing conservation strategies of these endangered species., (© 2024. The Author(s).)

Numerous methods for bulk RNA sequence deconvolution have been developed to identify cellular targets of diseases by understanding the composition of cell types in disease-related tissues. However, issues of heterogeneity in gene expression between subjects and the shortage of reference single-cell RNA sequence data remain to achieve accurate bulk deconvolution. In our study, we investigated whether a new data generative method named sc-CMGAN and benchmarking generative methods (Copula, CTGAN and TVAE) could solve these issues and improve the bulk deconvolutions. We also evaluated the robustness of sc-CMGAN using three deconvolution methods and four public datasets. In almost all conditions, the generative methods contributed to improved deconvolution. Notably, sc-CMGAN outperformed the benchmarking methods and demonstrated higher robustness. This study is the first to examine the impact of data augmentation on bulk deconvolution. The new generative method, sc-CMGAN, is expected to become one of the powerful tools for the preprocessing of bulk deconvolution., (© 2024. The Author(s).)

Chemotherapeutic drugs will affect the process of spermatogenesis. However, most current studies on the effects of chemotherapeutic drugs on spermatogenesis are based on mouse models, with a shortage of human body evidence. In addition, the mechanism of chemotherapeutic drugs causing spermatogenesis disorder is not clear. Therefore, we have collected the testicular tissues of an inguinal-lipoma patient whose testes were resected after chemotherapy and a patient who had normal spermatogenesis disorder and underwent single-nucleus RNA sequencing (snRNA-Seq). After quality control, we obtained a total of 27,957 high-quality cells, including 18,612 normal cells and 9,345 drug-treated cells, which were all used in analyzing the mechanism of chemotherapeutic drugs causing spermatogenesis disorder. This study has provided data resources and references for exploring the mechanism of chemotherapeutic drugs causing spermatogenesis disorder with the insight of protecting the spermatogenic abilities of male tumor patients receiving chemotherapy., (© 2024. The Author(s).)

Pathogen lineage nomenclature systems are a key component of effective communication and collaboration for researchers and public health workers. Since February 2021, the Pango dynamic lineage nomenclature for SARS-CoV-2 has been sustained by crowdsourced lineage proposals as new isolates were sequenced. This approach is vulnerable to time-critical delays as well as regional and personal bias. Here we developed a simple heuristic approach for dividing phylogenetic trees into lineages, including the prioritization of key mutations or genes. Our implementation is efficient on extremely large phylogenetic trees consisting of millions of sequences and produces similar results to existing manually curated lineage designations when applied to SARS-CoV-2 and other viruses including chikungunya virus, Venezuelan equine encephalitis virus complex and Zika virus. This method offers a simple, automated and consistent approach to pathogen nomenclature that can assist researchers in developing and maintaining phylogeny-based classifications in the face of ever-increasing genomic datasets., (© 2024. The Author(s).)

SPARC is an extracellular Ca 2+ -binding, secreted glycoprotein that plays a dynamic role in the growth and development of organisms. This study aimed to describe the isolation, characterization, and expression analysis of HdhSPARC in Pacific abalone (Haliotis discus hannai) to infer its potential functional role. The isolated HdhSPARC was 1633 bp long, encoding a polypeptide of 284 amino acid residues. Structurally, the SPARC protein in abalone is comprised of three biological domains. However, the structure of this protein varied between vertebrates and invertebrates, as suggested by their distinct clustering patterns in phylogenetic analysis. In early development, HdhSPARC was variably expressed, and higher expression was found in veliger larvae. Moreover, HdhSPARC was highly expressed in juvenile abalone with rapid growth compared to their slower-growing counterparts. Among the testicular development stages, the growth stage exhibited higher HdhSPARC expression. HdhSPARC was also upregulated during muscle remodeling and shell biomineralization, as well as in response to different stressors such as heat shock, LPS, and H 2 O 2 exposure. However, this gene was downregulated in Cd-exposed abalone. The present study first comprehensively characterized the HdhSPARC gene, and its spatio-temporal expressions were analyzed along with its responses to various stressors., (© 2024. The Author(s).)

The pink stem borer, Sesamia inferens (Walker), is a significant polyphagous pest historically restricted to regions south of N34° latitude. However, with changes in global climate and farming practices, the distribution of this moth has progressively exceeded its traditional limit of 34° N and encompassed most regions in North China. The genetic adaptations of S. inferens remain incompletely understood due to the lack of high-quality genome resources. Here, we sequenced the genome of S. inferens using PacBio and Hi-C technology, yielding a genome assembly of 865.04 Mb with contig N50 of 1.23 Mb. BUSCO analysis demonstrated this genome assembly has a high-level completeness of 96.1% gene coverage. In total, 459.72 Mb repeat sequences (53.14% of the assembled genome) and 20858 protein-coding genes were identified. We used the Hi-C technique to anchor 1135 contigs to 31 chromosomes, yielding a chromosome-level genome assembly with a scaffold N50 of 29.99 Mb. In conclusion, our high-quality genome assembly provided valuable resource that exploring the genetic characteristics of local adaptation and developing an efficient control strategy., (© 2024. The Author(s).)

Immunotherapy has emerged as a promising modality for addressing advanced or conventionally drug-resistant malignancies. When it comes to lung adenocarcinoma (LUAD), T cells have demonstrated significant influence on both antitumor activity and the tumor microenvironment. However, their specific contributions remain largely unexplored. This investigation aimed to delineate molecular subtypes and prognostic indicators founded on T cell marker genes, thereby shedding light on the significance of T cells in LUAD prognosis and precision treatment. The cellular phenotypes were identified by scrutinizing the single-cell data obtained from the GEO repository. Subsequently, T cell marker genes derived from single-cell sequencing analyses were integrated with differentially expressed genes from the TCGA repository to pinpoint T cell-associated genes. Utilizing Cox analysis, molecular subtypes and prognostic signatures were established and subsequently verified using the GEO dataset. The ensuing molecular and immunological distinctions, along with therapy sensitivity between the two sub-cohorts, were examined via the ESTIMATE, CIBERSORT, and ssGSEA methodologies. Compartmentalization, somatic mutation, nomogram development, chemotherapy sensitivity prediction, and potential drug prediction analyses were also conducted according to the risk signature. Additionally, real-time qPCR and the HPA database corroborated the mRNA and protein expression patterns of signature genes in LUAD tissues. In summary, this research yielded an innovative T cell marker gene-based signature with remarkable potential to prognosis and anticipate immunotherapeutic outcomes in LUAD patients., (© 2024. The Author(s).)

Movement is an important behavior observed in a wide range of taxa. Previous studies have examined genes controlling movement using wing polymorphic insects and genes controlling wing size. However, few studies have investigated genes controlling movement activity rather than morphological traits. In the present study, we conducted RNA sequencing using populations with higher (WL) and lower (WS) mobility established by artificial selection in the red flour beetle Tribolium castaneum and compared gene expression levels between selected populations with two replicate lines. As a result, we found significant differences between the selected populations in 677 genes expressed in one replicate line and 1198 genes expressed in another replicate line, of which 311 genes were common to the two replicate lines. Furthermore, quantitative PCR focusing on 6 of these genes revealed that neuropeptide F receptor gene (NpF) was significantly more highly expressed in the WL population than in the WS population, which was common to the two replicate lines. We discuss differences in genes controlling movement between walking activity and wing polymorphism., (© 2024. The Author(s).)

Circular RNAs (circRNAs) are linked to cancer, but it's still not clear what role they play in prostatic cancer. Through high-throughput sequencing, the goal of this study was to compare how circRNAs are expressed at different stages of prostate cancer. 12 patients attending the Department of Urology at the Second Affiliated Hospital of Nanchang University between June 2020 and October 2021 were used for RNA sequencing, and 14 patients were used for real-time fluorescent quantitative PCR (qRT-PCR). The expression profiles of prostate cancer circRNAs were constructed by sequencing with the help of next-generation high-throughput sequencing technology, and the differentially expressed circRNAs were analyzed by targeting microRNA (miRNA) loci and Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the genes from which circRNAs originated. Finally, the expression of target circRNAs in two prostate tissues was verified by qRT-PCR. Following high-throughput sequencing, 13,047 circRNAs were identified, and 605 circRNAs with significant differential expression were identified, of which 361 circRNAs were up-regulated, and 244 circRNAs were down-regulated. Analysis of circRNA-originated genes using GO and the KEGG enrichment analysis showed that circRNA host genes can regulate and influence multiple signaling pathways in prostate cancer with important biological functions. And the circRNA-miRNA network was constructed. The highest number of differentially expressed circRNA-binding miRNAs were: hsa&#95;circ&#95;000 7582 (52), hsa&#95;circ&#95;000 6198 (37), hsa&#95;circ&#95;000 6759 (28), hsa&#95;circ&#95;000 5675 (25), and hsa&#95;circ&#95;000 2172 (22). Moreover, we further screened out the circRNA (hsa&#95;circ&#95;0005692) that was significantly differentially expressed and common to all groups and verified by qRT-PCR that the expression of the target circRNA (hsa&#95;circ&#95;0005692) was significantly downregulated in prostate cancer compared with benign prostatic hyperplasia (BPH) tissues., (© 2023. The Author(s).)

FOXP3 is a transcription factor that is essential for the development of regulatory T cells, a branch of T cells that suppress excessive inflammation and autoimmunity 1-5 . However, the molecular mechanisms of FOXP3 remain unclear. Here we here show that FOXP3 uses the forkhead domain-a DNA-binding domain that is commonly thought to function as a monomer or dimer-to form a higher-order multimer after binding to T n G repeat microsatellites. The cryo-electron microscopy structure of FOXP3 in a complex with T 3 G repeats reveals a ladder-like architecture, whereby two double-stranded DNA molecules form the two 'side rails' bridged by five pairs of FOXP3 molecules, with each pair forming a 'rung'. Each FOXP3 subunit occupies TGTTTGT within the repeats in a manner that is indistinguishable from that of FOXP3 bound to the forkhead consensus motif (TGTTTAC). Mutations in the intra-rung interface impair T n G repeat recognition, DNA bridging and the cellular functions of FOXP3, all without affecting binding to the forkhead consensus motif. FOXP3 can tolerate variable inter-rung spacings, explaining its broad specificity for T n G-repeat-like sequences in vivo and in vitro. Both FOXP3 orthologues and paralogues show similar T n G repeat recognition and DNA bridging. These findings therefore reveal a mode of DNA recognition that involves transcription factor homomultimerization and DNA bridging, and further implicates microsatellites in transcriptional regulation and diseases., (© 2023. The Author(s).)

Histone post-translational modifications (PTMs) are associated with epigenetic states that form the basis for cell-type-specific gene expression.sup.1,2. Once established, histone PTMs can be maintained by positive feedback involving enzymes that recognize a pre-existing histone modification and catalyse the same modification on newly deposited histones. Recent studies suggest that in wild-type cells, histone PTM-based positive feedback is too weak to mediate epigenetic inheritance in the absence of other inputs.sup.3-7. RNA interference (RNAi)-mediated histone H3 lysine 9 methylation (H3K9me) and heterochromatin formation define a potential epigenetic inheritance mechanism in which positive feedback involving short interfering RNA (siRNA) amplification can be directly coupled to histone PTM positive feedback.sup.8-14. However, it is not known whether the coupling of these two feedback loops can maintain epigenetic silencing independently of DNA sequence and in the absence of enabling mutations that disrupt genome-wide chromatin structure or transcription.sup.15-17. Here, using the fission yeast Schizosaccharomyces pombe, we show that siRNA-induced H3K9me and silencing of a euchromatic gene can be epigenetically inherited in cis during multiple mitotic and meiotic cell divisions in wild-type cells. This inheritance involves the spreading of secondary siRNAs and H3K9me3 to the targeted gene and surrounding areas, and requires both RNAi and H3K9me, suggesting that the siRNA and H3K9me positive-feedback loops act synergistically to maintain silencing. By contrast, when maintained solely by histone PTM positive feedback, silencing is erased by H3K9 demethylation promoted by Epe1, or by interallelic interactions that occur after mating to cells containing an expressed allele even in the absence of Epe1. These findings demonstrate that the RNAi machinery can mediate transgenerational epigenetic inheritance independently of DNA sequence or enabling mutations, and reveal a role for the coupling of the siRNA and H3K9me positive-feedback loops in the protection of epigenetic alleles from erasure. In Schizosaccharomyces pombe, histone H3K9 methylation acts synergistically with short interfering RNA to perpetuate gene silencing during multiple mitotic and meiotic cell divisions., Author(s): Ruby Yu [sup.1] , Xiaoyi Wang [sup.1] , Danesh Moazed [sup.1] Author Affiliations: (1) Department of Cell Biology and Howard Hughes Medical Institute, Harvard Medical School, Boston, USA Main [...]

The rete testis (RT) is a region of the mammalian testis that plays an important role in testicular physiology. The RT epithelium consists of cells sharing some well-known gene markers with supporting Sertoli cells (SCs). However, little is known about the differences in gene expression between these two cell populations. Here, we used fluorescence-activated cell sorting (FACS) to obtain pure cultures of neonatal RT cells and SCs and identified differentially expressed genes (DEGs) between these cell types using RNA sequencing (RNA-seq). We then compared our data with the RNA-seq data of other studies that examined RT cells and SCs of mice of different ages and generated a list of DEGs permanently upregulated in RT cells throughout testis development and in culture, which included 86 genes, and a list of 79 DEGs permanently upregulated in SCs. The analysis of studies on DMRT1 function revealed that nearly half of the permanent DEGs could be regulated by this SC upregulated transcription factor. We suggest that useful cell lineage markers and candidate genes for the specification of both RT cells and SCs may be present among these permanent DEGs., (© 2023. The Author(s).)

Mitochondrial DNAs (mtDNAs) appear in almost all eukaryotic species and are useful molecular markers for phylogenetic studies and species identification. Kinetoplast DNAs (kDNAs) are structurally complex circular mtDNA networks in kinetoplastids, divided into maxicircles and minicircles. Despite several kDNAs of many Leishmania species being examined, the kDNAs of the new species, Leishmania orientalis (formerly named Leishmania siamensis) strain PCM2, have not been explored. This study aimed to investigate the maxicircle and minicircle DNAs of L. orientalis strain PCM2 using hybrid genome sequencing technologies and bioinformatic analyses. The kDNA sequences were isolated and assembled using the SPAdes hybrid assembler from the Illumina short-read and PacBio long-read data. Circular contigs of the maxicircle and minicircle DNAs were reconstructed and confirmed by BLASTn and rKOMICs programs. The kDNA genome was annotated by BLASTn before the genome comparison and phylogenetic analysis by progressiveMauve, MAFFT, and MEGA programs. The maxicircle of L. orientalis strain PCM2 (18,215 bp) showed 99.92% similarity and gene arrangement to Leishmania enriettii strain LEM3045 maxicircle with variation in the 12s rRNA gene and divergent region. Phylogenetics of the whole sequence, coding regions, divergent regions, and 12s rRNA gene also confirmed this relationship and subgenera separation. The identified 105 classes of minicircles (402-1177 bp) were clustered monophyletically and related to the Leishmania donovani minicircles. The kinetoplast maxicircle and minicircle DNAs of L. orientalis strain PCM2 contained a unique conserved region potentially useful for specific diagnosis of L. orientalis and further exploration of this parasite population genetics in Thailand and related regions., (© 2023. The Author(s).)

Understanding the influence of genetic variations in olfactory receptor (OR) genes on the olfaction-influenced phenotypes such as behaviors, reproduction, and feeding is important in animal biology. However, our understanding of the complexity of the OR subgenome is limited. In this study, we analyzed 1120 typing results of 20 representative OR genes belonging to 13 OR families on 14 pig chromosomes from 56 individuals belonging to seven different breeds using a sequence-based OR typing method. We showed that the presence of copy number variations, conservation of locus-specific diversity, abundance of breed-specific alleles, presence of a loss-of-function allele, and low-level purifying selection in pig OR genes could be common characteristics of OR genes in mammals. The observed nucleotide sequence diversity of pig ORs was higher than that of dogs. To the best of our knowledge, this is the first report on the individual- or population-level characterization of a large number of OR family genes in livestock species., (© 2023. The Author(s).)

The public health emergency caused by the Covid-19 outbreak in March 2020 encouraged worldwide initiatives to monitor the genetic diversity and features of the SARS-CoV-2 circulating variants, mainly based on the genomic surveillance. However, due to the impossibility to carry out extensive sequencing in resource-limited hospitals, other PCR-based strategies could be applied to efficiently monitor the circulating variants without the need to greatly expand the sequencing capacity. In our case, overpassing the technical limitations inherent to a second level hospital, we were able to characterize the weekly distribution of SARS-CoV-2 by the RT-qPCR amplification patterns visualization, single nucleotide polymorphism genotyping, and sequencing of randomly selected samples. All these molecular approaches allowed us to trace the epidemiology of SARS-CoV-2 viruses circulating in Ibiza and Formentera (Balearic Islands, Spain) during the third to the sixth pandemic waves (January 2021-July 2022), in which three major lineages that were considered as VOCs (Alpha, Delta, and Omicron), and many other non-VOC variants were detected and tracked., (© 2023. Springer Nature Limited.)

The Altai falcon from Central Asia always attracted the attention of humans. Long considered a totemic bird in its native area, modern falconers still much appreciated this large-bodied and mighty bird of prey due to its rarity and unique look. The peculiar body characteristics halfway between the saker falcon (Falco cherrug) and the gyrfalcon (F. rusticolus) triggered debates about its contentious taxonomy. The weak phylogenetic signal associated with traditional genetic methods could not resolve this uncertainty. Here, we address the controversial evolutionary origin of Altai falcons by means of a genome-wide approach, Restriction-site Associated DNA sequencing, using sympatric eastern sakers falcons, allopatric western saker falcons and gyrfalcons as outgroup. This approach provided an unprecedented insight into the phylogenetic relationships of the studied populations by delivering 17,095 unlinked SNPs shedding light on the polyphyletic nature of Altai falcons within eastern sakers. Thus we concluded that the former must correspond to a low taxonomic rank, probably an ecotype or form of the latter. Also, we found that eastern sakers are paraphyletic without gyrfalcons, thus, these latter birds are best regarded as the direct sister lineage of the eastern sakers. This evolutionary relationship, corroborated also by re-analyzing the dataset with the inclusion of outgroup samples (F. biarmicus and F. peregrinus), put eastern sakers into a new light as the potential ancestral genetic source of high latitude and altitude adaptation in descendent populations. Finally, conservation genomic values hint at the stable genetic background of the studied saker populations., (© 2023. Springer Nature Limited.)

Congenital disorders of glycosylation (CDG) are rare genetic disorders with a spectrum of clinical manifestations caused by abnormal N-glycosylation of secreted and cell surface proteins. Over 130 genes are implicated and next generation sequencing further identifies potential disease drivers in affected individuals. However, functional testing of these variants is challenging, making it difficult to distinguish pathogenic from non-pathogenic events. Using proximity labelling, we identified OST48 as a protein that transiently interacts with lysyl oxidase (LOX), a secreted enzyme that cross-links the fibrous extracellular matrix. OST48 is a non-catalytic component of the oligosaccharyltransferase (OST) complex, which transfers glycans to substrate proteins. OST48 is encoded by DDOST, and 43 variants of DDOST are described in CDG patients, of which 34 are classified as variants of uncertain clinical significance (VUS). We developed an assay based on LOX N-glycosylation that confirmed two previously characterised DDOST variants as pathogenic. Notably, 39 of the 41 remaining variants did not have impaired activity, but we demonstrated that p.S243F and p.E286del were functionally impaired, consistent with a role in driving CDG in those patients. Thus, we describe a rapid assay for functional testing of clinically relevant CDG variants to complement genome sequencing and support clinical diagnosis of affected individuals., (© 2023. Springer Nature Limited.)

Pancreatic cancer (PC) is known for its high degree of heterogeneity and exceptionally adverse outcome. While disulfidptosis is the most recently identified form of cell death, the predictive and therapeutic value of disulfidptosis-related genes (DRGs) for PC remains unknown. RNA sequencing data with the follow-up information, were retrieved from the TCGA and ICGC databases. Consensus clustering analysis was conducted on patient data using R software. Subsequently, the LASSO regression analysis was conducted to create a prognostic signature for foreseeing the outcome of PC. Differences in relevant pathways, mutational landscape, and tumor immune microenvironment were compared between PC samples with different risk levels. Finally, we experimentally confirmed the impact of DSG3 on the invasion and migration abilities of PC cells. All twenty DRGs were found to be hyperexpressed in PC tissues, and fourteen of them significantly associated with PC survival. Using consensus clustering analysis based on these DRGs, four DRclusters were identified. Additionally, altogether 223 differential genes were evaluated between clusters, indicating potential biological differences between them. Four gene clusters (geneClusters) were recognized according to these genes, and a 10-gene prognostic signature was created. High-risk patients were found to be primarily enriched in signaling pathways related to the cell cycle and p53. Furthermore, the rate of mutations was markedly higher in high-risk patients, besides important variations were present in terms of immune microenvironment and chemotherapy sensitivity among patients with different risk levels. DSG3 could appreciably enhance the invasion and migration of PC cells. This work, based on disulfidoptosis-related genes (DRGs), holds the promise of classifying PC patients and predicting their prognosis, mutational landscape, immune microenvironment, and drug therapy. These insights could boost an improvement in a better comprehension of the role of DRGs in PC as well as provide new opportunities for prognostic prediction and more effective treatment strategies., (© 2023. Springer Nature Limited.)

Prostate cancer (PCa) stands as a prominent contributor to morbidity and mortality among males on a global scale. Cancer-associated fibroblasts (CAFs) are considered to be closely connected to tumour growth, invasion, and metastasis. We explored the role and characteristics of CAFs in PCa through bioinformatics analysis and built a CAFs-based risk model to predict prognostic treatment and treatment response in PCa patients. First, we downloaded the scRNA-seq data for PCa from the GEO. We extracted bulk RNA-seq data for PCa from the TCGA and GEO and adopted "ComBat" to remove batch effects. Then, we created a Seurat object for the scRNA-seq data using the package "Seurat" in R and identified CAF clusters based on the CAF-related genes (CAFRGs). Based on CAFRGs, a prognostic model was constructed by univariate Cox, LASSO, and multivariate Cox analyses. And the model was validated internally and externally by Kaplan-Meier analysis, respectively. We further performed GO and KEGG analyses of DEGs between risk groups. Besides, we investigated differences in somatic mutations between different risk groups. We explored differences in the immune microenvironment landscape and ICG expression levels in the different groups. Finally, we predicted the response to immunotherapy and the sensitivity of antitumour drugs between the different groups. We screened 4 CAF clusters and identified 463 CAFRGs in PCa scRNA-seq. We constructed a model containing 10 prognostic CAFRGs by univariate Cox, LASSO, and multivariate Cox analysis. Somatic mutation analysis revealed that TTN and TP53 were significantly more mutated in the high-risk group. Finally, we screened 31 chemotherapeutic drugs and targeted therapeutic drugs for PCa. In conclusion, we identified four clusters based on CAFs and constructed a new CAFs-based prognostic signature that could predict PCa patient prognosis and response to immunotherapy and might suggest meaningful clinical options for the treatment of PCa., (© 2023. Springer Nature Limited.)

Carcinoma of unknown primary (CUP) is a type of metastatic cancer with tissue-of-origin (TOO) unidentifiable by traditional methods. CUP patients typically have poor prognosis but therapy targeting the original cancer tissue can significantly improve patients' prognosis. Thus, it's critical to develop accurate computational methods to infer cancer TOO. While qPCR or microarray-based methods are effective in inferring TOO for most cancer types, the overall prediction accuracy is yet to be improved. In this study, we propose a cross-cohort computational framework to trace TOO of 32 cancer types based on RNA sequencing (RNA-seq). Specifically, we employed logistic regression models to select 80 genes for each cancer type to create a combined 1356-gene set, based on transcriptomic data from 9911 tissue samples covering the 32 cancer types with known TOO from the Cancer Genome Atlas (TCGA). The selected genes are enriched in both tissue-specific and tissue-general functions. The cross-validation accuracy of our framework reaches 97.50% across all cancer types. Furthermore, we tested the performance of our model on the TCGA metastatic dataset and International Cancer Genome Consortium (ICGC) dataset, achieving an accuracy of 91.09% and 82.67%, respectively, despite the differences in experiment procedures and pipelines. In conclusion, we developed an accurate yet robust computational framework for identifying TOO, which holds promise for clinical applications. Our code is available at http://github.com/wangbo00129/classifybysklearn ., (© 2023. Springer Nature Limited.)

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 duplications 1-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 finished 4,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 genome 4 and mapped available population variation, clinical variants and functional genomics data to produce a complete and comprehensive reference sequence for all 24 human chromosomes., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)