Search

Microscopy and genomic approaches provide detailed descriptions of the three-dimensional folding of chromosomes and nuclear organization. The fundamental question is how activity of molecules at the nanometer scale can lead to complex and orchestrated spatial organization at the scale of chromosomes and the whole nucleus. At least three key mechanisms can bridge across scales: (1) tethering of specific loci to nuclear landmarks leads to massive reorganization of the nucleus; (2) spatial compartmentalization of chromatin, which is driven by molecular affinities, results in spatial isolation of active and inactive chromatin; and (3) loop extrusion activity of SMC (structural maintenance of chromosome) complexes can explain many features of interphase chromatin folding and underlies key phenomena during mitosis. Interestingly, many features of chromosome organization ultimately result from collective action and the interplay between these mechanisms, and are further modulated by transcription and topological constraints. Finally, we highlight some outstanding questions that are critical for our understanding of nuclear organization and function. We believe many of these questions can be answered in the coming years.

Our understanding of how genomic DNA is tightly packed inside the nucleus, yet is still accessible for vital cellular processes, has grown dramatically over recent years with advances in microscopy and genomics technologies. Computational methods have played a pivotal role in the structural interpretation of experimental data, which helped unravel some organizational principles of genome folding. Here, we give an overview of current computational efforts in mechanistic and data-driven 3D chromatin structure modeling. We discuss strengths and limitations of different methods and evaluate the added value and benefits of computational approaches to infer the 3D structural and dynamic properties of the genome and its underlying mechanisms at different scales and resolution, ranging from the dynamic formation of chromatin loops and topological associated domains to nuclear compartmentalization of chromatin and nuclear bodies.

Amoebozoa include lineages of diverse ecology, behavior, and morphology. They are assumed to encompass members with the largest genome sizes of all living things, yet genomic studies in the group are limited. Trichosphaerium, a polymorphic, multinucleate, marine amoeba with a complicated life cycle, has puzzled experts for over a century. In an effort to explore the genomic diversity and investigate extraordinary behavior observed among the Amoebozoa, we used integrated omics approaches to study this enigmatic marine amoeba. Omics data, including single-cell transcriptomics and cytological data, demonstrate that Trichosphaerium sp. possesses the complete meiosis toolkit genes. These genes are expressed in life stages of the amoeba including medium and large cells. The life cycle of Trichosphaerium sp. involves asexual processes via binary fission and multiple fragmentation of giant cells, as well as sexual-like processes involving genes implicated in sexual reproduction and polyploidization. These findings are in stark contrast to a life cycle previously reported for this amoeba. Despite the extreme morphological plasticity observed in Trichosphaerium, our genomic data showed that populations maintain a species-level intragenomic variation. A draft genome of Trichosphaerium indicates elevated lateral gene transfer (LGT) from bacteria and giant viruses. Gene trafficking in Trichosphaerium is the highest within Amoebozoa and among the highest in microbial eukaryotes.

Mammalian genomes have distinct levels of spatial organization and structure that have been hypothesized to play important roles in transcription regulation. Although much has been learned about these architectural features with ensemble techniques, single-cell studies are showing a new universal trend: Genomes are stochastic and dynamic at every level of organization. Stochastic gene expression, on the other hand, has been studied for years. In this review, we probe whether there is a causative link between the two phenomena. We specifically discuss the functionality of chromatin state, topologically associating domains (TADs), and enhancer biology in light of their stochastic nature and their specific roles in stochastic gene expression. We highlight persistent fundamental questions in this area of research.

Polyploidy is defined as a cell with three or more whole genome sets and enables cell growth across the kingdoms of life. Studies in model organisms have revealed that polyploid cell growth can be required for optimal tissue repair and regeneration. In mammals, polyploid cell growth contributes to repair of many tissues, including the liver, heart, kidney, bladder, and eye, and similar strategies have been identified in Drosophila and zebrafish tissues. This review discusses the heterogeneity and versatility of polyploidy in tissue repair and regeneration. Polyploidy has been shown to restore tissue mass and maintain organ size as well as protect against oncogenic insults and genotoxic stress. Polyploid cells can also serve as a reservoir for new diploid cells in regeneration. The numerous mechanisms to generate polyploid cells provide an unlimited resource for tissues to exploit to undergo repair or regeneration.

The calcareous shell and sessile lifestyle are the representative phenotypes of many molluscs, which happen to be present in barnacles, a group of unique crustaceans. The origin of these phenotypes is unclear, but it may be embodied in the convergent genetics of such distant groups (interphylum). Herein, we perform comprehensive comparative genomics analysis in barnacles and molluscs, and reveal a genome-wide strong convergent molecular evolution between them, including coexpansion of biomineralization and organic matrix genes for shell formation, and origination of lineage-specific orphan genes for settlement. Notably, the expanded biomineralization gene encoding alkaline phosphatase evolves a novel, highly conserved motif that may trigger the origin of barnacle shell formation. Unlike molluscs, barnacles adopt novel organic matrices and cement proteins for shell formation and settlement, respectively, and their calcareous shells have potentially originated from the cuticle system of crustaceans. Therefore, our study corroborates the idea that selection pressures driving convergent evolution may strongly act in organisms inhabiting similar environments regardless of phylogenetic distance. The convergence signatures shed light on the origin of the shell and sessile lifestyle of barnacles and molluscs. In addition, notable non-convergence signatures are also present and may contribute to morphological and functional specificities.

Influenza polymerase (FluPol) plays a key role in the viral infection cycle by transcribing and replicating the viral genome. FluPol is a multifunctional, heterotrimeric enzyme with cap-binding, endonuclease, RNA-dependent RNA polymerase and polyadenylation activities. It performs its functions in the context of the viral ribonucleoprotein particle (RNP), wherein the template viral RNA is coated by multiple copies of viral nucleoprotein. Moreover, it interacts with a number of host proteins that are essential cofactors and, consequently, adaptive mutations in the polymerase are required for crossing the avian-human species barrier. In this review, we show how mechanistic understanding of how FluPol performs its multiple functions has greatly advanced over the last decade through determination of high-resolution structures by X-ray crystallography and cryo-electron microscopy. These have revealed not only the detailed architecture of FluPol but highlighted the remarkably conformational flexibility that is inherent to its functioning as a dynamic RNA synthesis machine. Structural studies are also underpinning current attempts to develop next-generation anti-influenza drugs that directly target FluPol.

Testing the effect of rare variants on phenotypic variation is difficult due to the need for extremely large cohorts to identify associated variants given expected effect sizes. An alternative approach is to investigate the effect of rare genetic variants on low-level genomic traits, such as gene expression or DNA methylation (DNAm), as effect sizes are expected to be larger for low-level compared to higher-order complex traits. Here, we investigate DNAm in healthy ageing populations - the Lothian Birth cohorts of 1921 and 1936 and identify both transient and stable outlying DNAm levels across the genome. We find an enrichment of rare genetic variants within 1kb of DNAm sites in individuals with stable outlying DNAm, implying genetic control of this extreme variation. Using a family-based cohort, the Brisbane Systems Genetics Study, we observed increased sharing of DNAm outliers among more closely related individuals, consistent with these outliers being driven by rare genetic variation. We demonstrated that outlying DNAm levels have a functional consequence on gene expression levels, with extreme levels of DNAm being associated with gene expression levels towards the tails of the population distribution. Overall, this study demonstrates the role of rare variants in the phenotypic variation of low-level genomic traits, and the effect of extreme levels of DNAm on gene expression.

DNA introgressed from one species into another is typically deleterious at many genomic loci in the recipient species. It is therefore purged by selection over time. Here, we use mathematical modeling and whole-genome simulations to study the influence of recombination on the purging of introgressed DNA. We find that aggregate recombination controls the genome-wide rate of purging in the first few generations after admixture, when purging is most rapid. Aggregate recombination is quantified by, the average recombination rate across all locus pairs, and analogous metrics. It is influenced by the number of crossovers (i.e., the map length) and their locations along chromosomes, and by the number of chromosomes and heterogeneity in their size. A comparative prediction of our analysis is that species with fewer chromosomes should purge introgressed DNA more profoundly, and therefore should exhibit a weaker genomic signal of historical introgression. With regard to patterns across the genome, we show that, in heterogametic species with autosomal recombination in both sexes, more purging is expected on sex chromosomes than on autosomes, all else equal. The opposite prediction holds for species without autosomal recombination in the heterogametic sex. Finally, we show that positive genomic correlations between local recombination rate and introgressed ancestry, as recently observed in several taxa, are likely driven not by recombination’s effect in unlinking neutral from deleterious introgressed alleles, but rather by its effect on the rate of purging of the deleterious alleles themselves.Note on this versionAn earlier version of this manuscript had two parts: (1) Calculations of the variance of genetic relatedness between individuals with particular pedigree relationships, taking into account the randomness of recombination and segregation in their pedigree. (2) An investigation of the rate of purging of introgressed DNA following admixture, based in part on results from part (1). Part (1) has since been published as Veller et al. (2020). The present manuscript has been reconfigured to focus on part (2).

The maternal-to-zygotic transition (MZT) is a conserved embryonic process in animals where developmental control shifts from the maternal to zygotic genome. A key step in this transition is zygotic transcription, and deciphering the MZT requires classifying newly transcribed genes. However, due to current technological limitations, this starting point remains a challenge for studying many species. Here we present an alternative approach that characterizes transcriptome changes based solely on RNA-seq data. By combining intron-mapping reads and transcript-level quantification, we characterized transcriptome dynamics during the Drosophila melanogaster MZT. Our approach provides an accessible platform to investigate transcriptome dynamics that can be applied to the MZT in non-model organisms. In addition to classifying zygotically transcribed genes, our analysis revealed that over 300 genes express different maternal and zygotic transcript isoforms due to alternative splicing, polyadenylation, and promoter usage. The vast majority of these zygotic isoforms have the potential to be subject to different regulatory control, and over two-thirds encode different proteins. Thus, our analysis reveals an additional layer of regulation during the MZT, where new zygotic transcripts can generate additional proteome diversity.

Flowering time is an important agronomic trait of Chinese cabbage with late flowering being a primary breeding objective. In our previous work, we obtained Chinese cabbage-cabbage translocation lines that contained several beneficial cabbage genes. Cabbage-specific molecular markers show that these genes were coming from chromosome C01 of cabbage. In this study, we investigated the inheritance of flowering time in a couple of translocation lines and analyzed the transmission rate of molecular markers in the offspring. Consequently, we obtained the late flowering Chinese cabbage-cabbage translocation line ‘AT7-4’ in which the flowering time was later than that of ‘85-1’ by about 7 days under 4-week vernalization. Based on previous studies of the genomes of Chinese cabbage and cabbage, we located the cabbage-specific molecular markers that were closely linked at the top of the chromosome A01 in the F2 mapping population generated by self-crossing F1s derived from a cross between the translocation line ‘AT7-4’ and Chinese cabbage ‘14-36’. Five flowering-related genes in the alien fragment were found by functional annotation and their molecular markers were developed. This study lays the foundation for the future improvement of Chinese cabbage varieties using A-C translocation lines.

Two Streptomyces strains, named N11-26 and DC10-5, were isolated from deep-sea and non-photosynthetic stony coral, respectively. Strain N11-26 produces lobophorin C and divergolides, which are antimicrobial substances. This study aimed to classify these strains and reveal their cryptic potential to synthesize other secondary metabolites, such as polyketides and nonribosomal peptides. Strains N11-26 and DC10-5 showed 16S rRNA gene sequence similarities of 100% and 99.9% to Streptomyces olivaceus NRRL B-3009T, respectively. By digital DNA–DNA hybridization using whole-genome sequences, these strains were classified as Streptomyces olivaceus. Strain N11-26 was closer to the type strain of S. olivaceus than strain DC10-5 and possessed 17 clusters of polyketide synthase (PKS) and/or nonribosomal peptide synthetases (NRPS) genes, whereas strain DC10-5 harbored 19 clusters. Putative products by these gene clusters were predicted by bioinformatic analyses. Although 15 clusters were conserved between the two strains, two and four clusters were specific in strains N11-26 and DC10-5, respectively. This represents a diversity of potential polyketide and nonribosomal peptide compounds between strains of S. olivaceus. To the best of our knowledge, this is the first report annotating all the PKS and NRPS gene clusters in S. olivaceus strains with their putative products to provide useful information for genome mining.

SUMMARY OBJECTIVE: Genome sequencing has been proved to be an excellent tool to monitor the molecular epidemiology of the disease caused by severe acute respiratory syndrome coronavirus 2, i.e., coronavirus disease 2019. Some reports of infected, vaccinated individuals have aroused great interest because they are primarily being infected with circulating variants of concern. To investigate the cases of infected, vaccinated individuals in Salvador, Bahia, Brazil, we performed genomic monitoring to estimate the magnitude of the different variants of concern in these cases. METHODS: Nasopharyngeal swabs from infected (symptomatic and asymptomatic), vaccinated or unvaccinated individuals (n=29), and quantitative reverse transcription polymerase chain reaction cycle threshold value (Ct values) of ≤30 were subjected to viral sequencing using nanopore technology. RESULTS: Our analysis revealed that the Omicron variant was found in 99% of cases and the Delta variant was found in only one case. Infected, fully vaccinated patients have a favorable clinical prognosis; however, within the community, they become viral carriers with the aggravating factor of viral dissemination of variants of concern not neutralized by the currently available vaccines. CONCLUSION: It is important to acknowledge the limitations of these vaccines and to develop new vaccines to emergent variants of concern, as is the case of influenza vaccine; going through new doses of the same coronavirus vaccines is “more of the same.”

A novel bacterial strain, designated GeG2T, was isolated from soils of native Cerrado, a highly biodiverse savanna-like Brazilian biome. 16S rRNA gene sequence analysis of strain GeG2T revealed high sequence identity (100%) to the alphaproteobacterium Novosphingobium rosa, however, comparisons with N. rosa DSM7285T showed several distinctive features, prompting a full characterization of the new strain in terms of growth, morphology, biochemistry and, ultimately, its genome. GeG2T cells were Gram-stain negative bacilli, facultatively anaerobic, motile, positive for catalase and oxidase activities and for starch hydrolysis. Strain GeG2T presented planktonic-sessile dimorphism and cell aggregates surrounded by extracellular matrix and nanometric spherical structures were observed in liquid cultures, suggesting the production of exopolysaccharides (EPS) and outer membrane vesicles (OMVs). Whole genome assembly revealed four circular replicons: a 4.1 Mb chromosome, a 2.7 Mb extrachromosomal megareplicon and two plasmids (212.7 and 68.6 kb). The megareplicon contains few core genes and plasmid-type replication/maintenance systems, consistent with its classification as a chromid. Genome annotation shows a vast repertoire of carbohydrate active enzymes and genes involved in the degradation of aromatic compounds, highlighting the biotechnological potential of the new isolate obtained from Cerrado soils, especially regarding EPS production and biodegradation of recalcitrant compounds. Chemotaxonomic features, including polar lipid and fatty acid profiles, as well as physiological, molecular and whole genome comparisons showed significant differences between strain GeG2T and a N. rosa, clearly indicating that it represents a novel species, for which the name Novosphingobium terrae is proposed. The type strain is GeG2T (=CBMAI 2313T =CBAS 753T).IMPORTANCENovosphingobium is an alphaproteobacterial genus presenting diverse physiological profiles and broad biotechnological applications. However, many aspects regarding the biology of this important bacterial group remain elusive. A novel Novosphingobium strain was isolated from soils of Cerrado, an important Brazilian biome. Despite 100% 16S rRNA gene identity with Novosphingobium rosa, polyphasic characterizations, including physiological, chemotaxonomic, and whole genome- based analyses revealed significant differences between GeG2T and N. rosa DSM7285T, reinforcing resolution limitations in phylogenetic analysis based solely on 16S RNA and highlighting the importance of employing different approaches for the description of bacterial species. Using short and long read sequencing approaches, a high-quality fully resolved genome assembly was generated and one of the largest chromids reported to date was identified. A comprehensive characterization of environmental isolates allows us to better elucidate the diversity and biology of members of this bacterial group with potential biotechnological importance, guiding future bioprospecting efforts and genomic studies.

Pangenomics is emerging as a powerful computational paradigm in bioinformatics. This field uses population-level genome reference structures, typically consisting of a sequence graph, to mitigate reference bias and facilitate analyses that were challenging with previous reference-based methods. In this work, we extend these methods into transcriptomics to analyze sequencing data using the pantranscriptome: a population-level transcriptomic reference. Our novel toolchain can construct spliced pangenome graphs, map RNA-seq data to these graphs, and perform haplotype-aware expression quantification of transcripts in a pantranscriptome. This workflow improves accuracy over state-of-the-art RNA-seq mapping methods, and it can efficiently quantify haplotype-specific transcript expression without needing to characterize a sample’s haplotypes beforehand.

The nuclear envelope (NE) protects but also organizes the eukaryotic genome. In this review we will discuss recent literature on how the NE disassembles and reassembles, how it varies in surface area and protein composition and how this translates into chromatin organization and gene expression in the context of animal development.

Extremophiles have always garnered great interest because of their exotic lifestyles and ability to thrive at the physical limits of life. In hot springs environments, the Cyanidiophyceae red algae are the only photosynthetic eukaryotes able to live under extremely low pH (0-5) and relatively high temperature (35ºC to 63ºC). These extremophiles live as biofilms in the springs, inhabit acid soils near the hot springs, and form endolithic populations in the surrounding rocks. Cyanidiophyceae represent a remarkable source of knowledge about the evolution of extremophilic lifestyles and their genomes encode specialized enzymes that have applied uses. Here we review the evolutionary origin, taxonomy, genome biology, industrial applications, and use of Cyanidiophyceae as genetic models. Currently, Cyanidiophyceae comprise a single order (Cyanidiales), three families, four genera, and nine species, including the well-known Cyanidioschyzon merolae and Galdieria sulphuraria. These algae have small, gene-rich genomes that are analogous to those of prokaryotes they live and compete with. There are few spliceosomal introns and evidence exists for horizontal gene transfer as a driver of local adaptation to gain access to external fixed carbon and to extrude toxic metals. Cyanidiophyceae offer a variety of commercial opportunities such as phytoremediation to detoxify contaminated soils or waters and exploitation of their mixotrophic lifestyles to support the efficient production of bioproducts such as phycocyanin and floridosides. In terms of exobiology, Cyanidiophyceae are an ideal model system for understanding the evolutionary effects of foreign gene acquisition and the interactions between different organisms inhabiting the same harsh environment on the early Earth. Finally, we describe ongoing research with C. merolae genetics and summarize the unique insights they offer to the understanding of algal biology and evolution.

The large yellow croaker (Larimichthys crocea) is an economically important marine species with the highest annual production among the farmed marine fishes in China. However, the aquaculture industry of this species is suffering from severe problems that include weakened disease resistance, decreased growth rate, and reduced meat quality due to frequent inbreeding. Genome editing, which has a huge potential for solving those problems by introducing favorable genetic changes, is not yet available for the large yellow croaker. Here, we pioneered the techniques of embryo microinjection and genome editing using the CRISPR/Cas9 system in this species. Recombinant plasmids encoding green fluorescent protein (GFP) were introduced into the fertilized eggs of L. crocea by microinjection before the chorion had hardened. A high survival rate (40%) and GFP-positive larvae rate (81.8%) were achieved, indicating that the microinjection technique in L. crocea was successfully established. On this basis, Cas9 mRNA and sgRNA targeting the tyrosinase a gene in L. crocea (Lc-tyra) were co-injected into fertilized eggs of L. crocea. Mutant individuals with insertion and deletion mutations of Lc-tyra were detected. These results indicated that the CRISPR/Cas9-based genome editing technology established herein could efficiently introduce mutations at a specific site in the L. crocea genome. This method provides the potential for genetic improvement and functional genomic study in this species. This is the first report on effective CRISPR/Cas9-based genome editing in L. crocea.

International audience; In brown algae, the extracellular matrix (ECM) and its constitutive polymers play crucial roles in specialized functions, including algal growth and development. In this review we offer an integrative view of ECM construction in brown algae. We briefly report the chemical composition of its main constituents, and how these are interlinked in a structural model. We examine the ECM assembly at the tissue and cell level, with consideration on its structure in vivo and on the putative subcellular sites for the synthesis of its main constituents. We further discuss the biosynthetic pathways of two major polysaccharides, alginates and sulfated fucans, and the progress made beyond the candidate genes with the biochemical validation of encoded proteins. Key enzymes involved in the elongation of the glycan chains are still unknown and predictions have been made at the gene level. Here, we offer a re-examination of some glycosyltransferases and sulfotransferases from published genomes. Overall, our analysis suggests novel investigations to be performed at both the cellular and biochemical levels. First, to depict the location of polysaccharide structures in tissues. Secondly, to identify putative actors in the ECM synthesis to be functionally studied in the future.

T cells are essential for mounting defense against various pathogens and malignantly transformed cells. Thymic development and peripheral T cell differentiation are highly orchestrated biological processes that require precise gene regulation. Higher-order genome organization on multiple scales, in the form of chromatin loops, topologically associating domains and compartments, provides pivotal control of T cell gene expression. CTCF and the cohesin machinery are ubiquitously expressed architectural proteins responsible for establishing chromatin structures. Recent studies indicate that transcription factors, such as T lineage–defining Tcf1 and TCR-induced Batf, may have intrinsic ability and/or engage CTCF to shape chromatin architecture. In this article, we summarize current knowledge on the dynamic changes in genome topology that underlie normal or leukemic T cell development, CD4+ helper T cell differentiation, and CD8+ cytotoxic T cell functions. The knowledge lays a solid foundation for elucidating the causative link of spatial chromatin configuration to transcriptional and functional output in T cells.

Germline mutation rates are essential for genetic and evolutionary analyses. Yet, estimating accurate fine-scale mutation rates across the genome is a great challenge, due to relatively few observed mutations and intricate relationships between predictors and mutation rates. Here we present MuRaL (Mutation Rate Learner), a deep learning framework to predict mutation rates at the nucleotide level using only genomic sequences as input. Harnessing human germline variants for comprehensive assessment, we show that MuRaL achieves better predictive performance than current state-of-the-art methods. Moreover, MuRaL can build models with relatively few training mutations and a moderate number of sequenced individuals, and can leverage transfer learning to further reduce data and time demands. We apply MuRaL to produce genome-wide mutation rate maps for four representative species - Homo sapiens, Macaca mulatta, Arabidopsis thaliana and Drosophila melanogaster, demonstrating its high applicability. As an example, we use improved mutation rate estimates to stratify human genes into distinct groups which are enriched for different functions, and highlight that many developmental genes are subject to high mutational burden. The open-source software and generated mutation rate maps can greatly facilitate related research.

Introduction Legume crops are an important source of protein and oil for human health and in fixing atmospheric N2 for soil enrichment. With an objective to accelerate much-needed genetic analyses and breeding applications, draft genome assemblies were generated in several legume crops; many of them are not high quality because they are mainly based on short reads. However, the superior quality of genome assembly is crucial for a detailed understanding of genomic architecture, genome evolution, and crop improvement. Objectives Present study was undertaken with an objective of developing improved chromosome-length genome assemblies in six different legumes followed by their systematic investigation to unravel different aspects of genome organization and legume evolution. Methods We employed in situ Hi-C data to improve the existing draft genomes and performed different evolutionary and comparative analyses using improved genome assemblies. Results We have developed chromosome-length genome assemblies in chickpea, pigeonpea, soybean, subterranean clover, and two wild progenitor species of cultivated groundnut (A. duranensis and A. ipaensis). A comprehensive comparative analysis of these genome assemblies offered improved insights into various evolutionary events that shaped the present-day legume species. We highlighted the expansion of gene families contributing to unique traits such as nodulation in legumes, gravitropism in groundnut, and oil biosynthesis in oilseed legume crops such as groundnut and soybean. As examples, we have demonstrated the utility of improved genome assemblies for enhancing the resolution of “QTL-hotspot” identification for drought tolerance in chickpea and marker-trait associations for agronomic traits in pigeonpea through genome-wide association study. Genomic resources developed in this study are publicly available through an online repository, ‘Legumepedia’. Conclusion This study reports chromosome-length genome assemblies of six legume species and demonstrates the utility of these assemblies in crop improvement. The genomic resources developed here will have significant role in accelerating genetic improvement applications of legume crops.

The increasing availability of high-throughput datasets allows amalgamating research information across a large body of genome regulation studies. Given the recent success of meta-analyses on transcriptional regulators, epigenetic marks, and enhancer:gene associations, we expect that such surveys will continue to provide novel and reproducible insights. However, meta-analyses are severely hampered by the diversity of available data, concurring protocols, an eclectic amount of bioinformatics tools, and myriads of conceivable parameter combinations. Such factors can easily bar life scientists from synthesizing omics data and substantially curb their interpretability. Despite statistical challenges of the method, we would like to emphasize the advantages of joining data from different sources through vote-counting and showcase examples that achieve a simple but highly intuitive data integration.

Minimizers are widely used to sample representative

Peroxisome proliferator-activated receptor (PPAR) plays an important role in the regulation of lipid metabolism and has been widely identified in diverse species. Gymnocypris przewalskii is a native fish of the Qinghai Tibetan Plateau that survives in a chronically cold environment. In the current study, we conducted genome-wide identification of PPAR genes, revealing the existence of seven PPARs in the G. przewalskii genome. Collinearity was observed between two copies of PPARαb and PPARγ in G. przewalskii, suggesting that the additional copy might be gained through whole genome duplication. Both phylogenetic and multiple sequence alignment analyses indicated that PPARs in G. przewalskii were conserved with teleosts. The cold treatment (10 °C and 4 °C) led to the developmental delay of G. przewalskii embryos. Continuous expression of PPARs was observed during the embryonic development of G. przewalskii under normal and cold conditions, with significantly different transcriptional patterns. These results indicated that PPARs participated in the embryonic development of G. przewalskii, and were involved in the cold response during development. The current study proposed a potential role of PPARs in the cold response in the embryonic development of G. przewalskii, which shed light on understanding cold adaptation in Tibetan highland fish.

The recently established prime editor (PE) system is regarded as next-generation gene-editing technology. This methodology can install any base-to-base change as well as insertions and deletions without the requirement for double-stranded break formation or donor DNA templates; thus, it offers more targeting flexibility and greater editing precision than conventional CRISPR-Cas systems or base editors. In this study, we introduce the basic principles of PE and then review its most recent progress in terms of editing versatility, specificity, and efficiency in mammals. Next, we summarize key considerations regarding the selection of PE variants, prime editing guide RNA (pegRNA) design rules, and the efficiency and accuracy evaluation of PE. Finally, we highlight and discuss how PE can assist in a wide range of biological studies and how it can be applied to make precise genomic corrections in animal models, which paves the way for curing human diseases.

The rubber boa, Charina bottae is a semi-fossorial, cold-temperature adapted snake that ranges across the wetter and cooler ecoregions of the California Floristic Province. The rubber boa is 1 of 2 species in the family Boidae native to California and currently has 2 recognized subspecies, the Northern rubber boa C. bottae bottae and the Southern rubber boa C. bottae umbratica. Recent genomic work on C. bottae indicates that these 2 subspecies are collectively composed of 4 divergent lineages that separated during the late Miocene. Analysis of habitat suitability indicates that C. bottae umbratica montane sky-island populations from southern California will lose the majority of their habit over the next 70 yr, and is listed as Threatened under the California Endangered Species Act. Here, we report a new, chromosome-level assembly of C. bottae bottae as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genome strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises 289 scaffolds covering 1,804,944,895 bp, has a contig N50 of 37.3 Mb, a scaffold N50 of 97 Mb, and BUSCO completeness score of 96.3%, and represents the first reference genome for the Boidae snake family. This genome will enable studies of genetic differentiation and connectivity among C. bottae bottae and C. bottae umbratica populations across California and help manage locally endemic lineages as they confront challenges from human-induced climate warming, droughts, and wildfires across California.

Joint genomic prediction (GP) is an attractive method to improve the accuracy of GP by combining information from multiple populations. However, many factors can negatively influence the accuracy of joint GP, such as differences in linkage disequilibrium phasing between single nucleotide polymorphisms (SNPs) and causal variants, minor allele frequencies and causal variants' effect sizes across different populations. The objective of this study was to investigate whether the imputed high-density genotype data can improve the accuracy of joint GP using genomic best linear unbiased prediction (GBLUP), single-step GBLUP (ssGBLUP), multi-trait GBLUP (MT-GBLUP) and GBLUP based on genomic relationship matrix considering heterogenous minor allele frequencies across different populations (wGBLUP). Three traits, including days taken to reach slaughter weight, backfat thickness and loin muscle area, were measured on 67 276 Large White pigs from two different populations, for which 3334 were genotyped by SNP array. The results showed that a combined population could substantially improve the accuracy of GP compared with a single-population GP, especially for the population with a smaller size. The imputed SNP data had no effect for single population GP but helped to yield higher accuracy than the medium-density array data for joint GP. Of the four methods, ssGLBUP performed the best, but the advantage of ssGBLUP decreased as more individuals were genotyped. In some cases, MT-GBLUP and wGBLUP performed better than GBLUP. In conclusion, our results confirmed that joint GP could be beneficial from imputed high-density genotype data, and the wGBLUP and MT-GBLUP methods are promising for joint GP in pig breeding.

The use of genotype information to improve the predictability of Expected Progeny Difference was first implemented in American Angus cattle in 2009 and has now grown to where over 50% of all registered calves are genotyped. Animals with only a genotype now have genetic prediction accuracy equivalent to eight or more progeny records across all traits. Reproductive technologies have also been widely adopted with approximately 50% of all calves born being the result of artificial insemination. Non-surgical embryo transfer started increasing in the mid 1990s with just over 10% of calves born being the result of embryo transfer since 2005. The number of embryos created with in vitro technologies has risen sharply since 2015 and now accounts for close to 30% of all ET calves. Genomics has enabled embryo technologies to be more impactful, as females can be selected with greater accuracy and sires can be used at earlier ages with moderate accuracy. Large numbers of females genotyped each year also increases the number of selection candidates, increasing the selection intensity. Genomics, combined with increased recording, also provides more information on females. This increases the spread in the estimated index values of current dams, identifying more elite dams for selection as embryo donors. The greater scope of female selection also contributes to better inbreeding management. Commercial animals genotyped could be targeted for oocyte harvesting at slaughter, creating opportunities for low cost high value beef embryos to be used in the beef on dairy segment of the industry.

Motivation Years of time-series gene expression studies have built a strong understanding of clock-controlled pathways across species. However, comparatively little is known about how ‘non-clock’ pathways influence clock function. We need a strong understanding of clock-coupled pathways in human tissues to better appreciate the links between disease and clock function. Results We developed a new computational approach to explore candidate pathways coupled to the clock in human tissues. This method, termed LTM, is an in silico screen to infer genetic influences on circadian clock function. LTM uses natural variation in gene expression in human data and directly links gene expression variation to clock strength independent of longitudinal data. We applied LTM to three human skin and one melanoma datasets and found that the cell cycle is the top candidate clock-coupled pathway in healthy skin. In addition, we applied LTM to thousands of tumor samples from 11 cancer types in the TCGA database and found that extracellular matrix organization-related pathways are tightly associated with the clock strength in humans. Further analysis shows that clock strength in tumor samples is correlated with the proportion of cancer-associated fibroblasts and endothelial cells. Therefore, we show both the power of LTM in predicting clock-coupled pathways and classify factors associated with clock strength in human tissues. Availability and implementation LTM is available on GitHub (https://github.com/gangwug/LTMR) and figshare (https://figshare.com/articles/software/LTMR/21217604) to facilitate its use. Supplementary information Supplementary data are available at Bioinformatics online.

The blackspotted croaker (Protonibea diacanthus) is an endangered coastal marine fish. It is also a valuable species that is cultured on the southeast coast of China. While some genetic studies have been conducted to protect this species, genomic resources are lacking. Here, we report a chromosome-scale assembly of P. diacanthus genome by high-depth genome sequencing, assembly, and annotation. The genome scale was 635.69 Mb with contig and scaffold N50 length of 3.33 Mb and 25.60 Mb, respectively. Hi-C scaffolding of the genome resulted in 24 chromosomes of 94.15% total genome. We predicted 23,971 protein-coding genes. In addition, we constructed a phylogenetic tree using 2755 single-copy gene families and identified 462 unique gene families in P. diacanthus genome compared to three other sciaenids. What's more, from the analysis of gene families, we found that several gene families related to innate immunity were significantly expanded in the blackspotted croaker genome compared to other teleost genomes. The high-quality genome can improve our understanding of the molecular mechanisms behind economically valuable traits and provide insights into characteristics of the immune system.

One gene does not equal one RNA or protein. The genomic revolution has revealed numerous different RNA and protein molecules that can be produced from one gene, such as circular RNAs generated by back-splicing, proteins with residues mismatching the genomic encoding because of RNA editing, and proteins extended in the C terminus via stop codon readthrough in translation. Are these diverse products results of exquisite gene regulations or imprecise biological processes? While there are cases where the gene product diversity appears beneficial, genome-scale patterns suggest that much of this diversity arises from nonadaptive, molecular errors. This finding has important implications for studying the functions of diverse gene products and for understanding the fundamental properties and evolution of cellular life.

In this article, a bioinformatics pipeline is proposed that focuses on two types of elements, namely the mobile genetic elements (MGE) and Ribonucleic acids (RNAs). The MGEs are called insertion sequences (ISs) in the prokaryotic domain. The objective of this research work is to study the behaviour of RNAs and MGEs genes, and the effects of their presence around inversions in genome sequences. The proposed pipeline finds the relation between the transposase gene types (e.g., DDE and DEDD) located within insertion sequences according to their IS family and sub-family, and RNAs (tRNA and rRNA) on the one hand, and genomic inversion on the other hand. More precisely, we wonder whether the insertion sequences and RNAs have any effects on the occurrence of inversions. To investigate this, we start by extracting all inversions between every two close isolates of P. aeruginosa from a high supported phylo- genetic tree, and we investigate inversion boundaries, to determine the impact of IS elements and RNAs around these inversions. For RNAs genes, we noticed that a high occurrence of these elements are found in what we called InBoundary, while few to zero occurrences have been recorded in OutBoundary. We remark that the NCGM2.S1 genome has the largest amount of tRNAs and inversions comparing to other genomes. For RNA genes, the tRNAs genes (Arg, Ala, Leu) recorded the highest occurrences around the inversions, while there was few occurrences of rRNA genes, within all genomes. For IS elements, we noticed that the family of IS3 has the highest impact on inversions, but sometimes other type of IS elements can also be found with IS3 around the inversions. In this case, these types of ISs shared the same type of transposase gene (DDE). We finally conclude that the distribution of RNA and IS genes has a big impact on the inversions in bacterial genomes.

Point mutations and structural variants that directly disrupt the coding sequence of MEF2C have been associated with a spectrum of neurodevelopmental disorders (NDDs). However, the impact of MEF2C haploinsufficiency on neurodevelopmental pathways and synaptic processes is not well understood, nor are the complex mechanisms that govern its regulation. To explore the functional changes associated with structural variants that alter MEF2C expression and/or regulation, we generated an allelic series of 204 isogenic human induced pluripotent stem cell (hiPSC)-derived neural stem cells and glutamatergic induced neurons. These neuronal models harbored CRISPR-engineered mutations that involved direct deletion of MEF2C or deletion of the boundary points for topologically associating domains (TADs) and chromatin loops encompassing MEF2C. Systematic profiling of mutation-specific alterations, contrasted to unedited controls that were exposed to the same guide RNAs for each edit, revealed that deletion of MEF2C caused differential expression of genes associated with neurodevelopmental pathways and synaptic function. We also discovered significant reduction in synaptic activity measured by multielectrode arrays (MEAs) in neuronal cells. By contrast, we observed robust buffering against MEF2C regulatory disruption following deletion of a distal 5q14.3 TAD and loop boundary, whereas homozygous loss of a proximal loop boundary resulted in down-regulation of MEF2C expression and reduced electrophysiological activity on MEA that was comparable to direct gene disruption. Collectively, these studies highlight the considerable functional impact of MEF2C deletion in neuronal cells and systematically characterize the complex interactions that challenge a priori predictions of regulatory consequences from structural variants that disrupt three-dimensional genome organization.

Summary Viruses circulating in wild and domestic animals pose a constant threat to human health1. Identifying human genetic factors that protect against zoonotic infections is a health priority. The RD-114 and Type-D retrovirus (RDR) interference group includes infectious viruses that circulate in domestic cats and various Old World monkeys (OWM), and utilize ASCT2 as a common target cell receptor2. While human ASCT2 can mediate RDR infection in cell culture, it is unknown whether humans and other hominoids encode factors that restrict RDR infection in nature2,3. Here we test the hypothesis that Suppressyn, a truncated envelope protein that binds ASCT2 and is derived from a human endogenous retrovirus4,5, restricts RDR infection. Transcriptomics and regulatory genomics reveal that Suppressyn expression initiates in the preimplantation embryo. Loss and gain of function experiments in cell culture show Suppressyn expression is necessary and sufficient to restrict RDR infection. Evolutionary analyses show Suppressyn was acquired in the genome of a common ancestor of hominoids and OWMs, but preserved by natural selection only in hominoids. Restriction assays using modern primate orthologs and reconstructed ancestral genes indicate that Suppressyn antiviral activity has been conserved in hominoids, but lost in most OWM. Thus in humans and other hominoids, Suppressyn acts as a restriction factor against retroviruses with zoonotic capacity. Transcriptomics data predict that other virus-derived proteins with potential antiviral activity lay hidden in the human genome.

Large-scale whole-genome sequencing studies have enabled analysis of noncoding rare-variant (RV) associations with complex human diseases and traits. Variant-set analysis is a powerful approach to study RV association. However, existing methods have limited ability in analyzing the noncoding genome. We propose a computationally efficient and robust noncoding RV association detection framework, STAARpipeline, to automatically annotate a whole-genome sequencing study and perform flexible noncoding RV association analysis, including gene-centric analysis and fixed window-based and dynamic window-based non-gene-centric analysis by incorporating variant functional annotations. In gene-centric analysis, STAARpipeline uses STAAR to group noncoding variants based on functional categories of genes and incorporate multiple functional annotations. In non-gene-centric analysis, STAARpipeline uses SCANG-STAAR to incorporate dynamic window sizes and multiple functional annotations. We apply STAARpipeline to identify noncoding RV sets associated with four lipid traits in 21,015 discovery samples from the Trans-Omics for Precision Medicine (TOPMed) program and replicate several of them in an additional 9,123 TOPMed samples. We also analyze five non-lipid TOPMed traits.

Evolution can repeat itself, resulting in parallel adaptations in independent lineages occupying similar environments. Moreover, parallel evolution sometimes, but not always, uses the same genes. Two main hypotheses have been put forth to explain the probability and extent of parallel evolution. First, parallel evolution is more likely when shared ecologies result in similar patterns of natural selection in different taxa. Second, parallelism is more likely when genomes are similar because of shared standing variation and similar mutational effects in closely related genomes. Here we combine ecological, genomic, experimental and phenotypic data with Bayesian modelling and randomization tests to quantify the degree of parallelism and its relationship with ecology and genetics. Our results show that the extent to which genomic regions associated with climate are parallel among species of Timema stick insects is shaped collectively by shared ecology and genomic background. Specifically, the extent of genomic parallelism decays with divergence in climatic conditions (that is, habitat or ecological similarity) and genomic similarity. Moreover, we find that climate-associated loci are likely subject to selection in a field experiment, overlap with genetic regions associated with cuticular hydrocarbon traits and are not strongly shaped by introgression between species. Our findings shed light on when evolution is most expected to repeat itself.

Psychiatric disorders and other complex traits have a polygenic architecture, often associated with dozens or even hundreds of independent genomic loci. As each of these have a relatively small influence on the trait, the dissection of their biological components is a non-trivial task. For psychiatric disorders in particular, the majority of associated loci lie within non-coding regions of the genome, suggesting that most of the genetic risk for disease originates from the disruption of regulatory sequences. While previously exploration of the heritability of these sequences has focused on variants that modify DNA elements, those that alter cis-acting RNA sequences, such as miRNA binding sites, are also likely to have a significant impact in these disorders. MiRNA have already been shown to be dysregulated in these disorders through both genetic and environmental influence, so it is reasonable to suspect their target genes may also be affected by common variation. In this study, we investigated the distribution of miRNA binding site variants (MBSVs) predicted to alter miRNA binding affinity in psychiatric disorders and observed significant enrichment in schizophrenia, depression, bipolar disorder, and anorexia nervosa. We also observed significant enrichment of MBSVs in genes targeted by several miRNA families, including miR-335-5p, miR-21-5p/590-5p, miR-361-5p, and miR-557 in both schizophrenia and depression, and nominally significant enrichment of MBSV for miR-323b-3p in schizophrenia. We also identified a significant association between MBSVs in gene sets involved in regulation of the synapse and synaptic depression in schizophrenia. While these observations support the role of miRNA in the pathophysiology of psychiatric disorders, we also observed significant association of MBSVs in other complex traits suggesting that MBSVs are an important class of regulatory variants that have functional implications for many disorders.

Genomic analyses of Neanderthals have previously provided insights into their population history and relationship to modern humans1–8, but the social organization of Neanderthal communities remains poorly understood. Here we present genetic data for 13 Neanderthals from two Middle Palaeolithic sites in the Altai Mountains of southern Siberia: 11 from Chagyrskaya Cave9,10 and 2 from Okladnikov Cave11—making this one of the largest genetic studies of a Neanderthal population to date. We used hybridization capture to obtain genome-wide nuclear data, as well as mitochondrial and Y-chromosome sequences. Some Chagyrskaya individuals were closely related, including a father–daughter pair and a pair of second-degree relatives, indicating that at least some of the individuals lived at the same time. Up to one-third of these individuals’ genomes had long segments of homozygosity, suggesting that the Chagyrskaya Neanderthals were part of a small community. In addition, the Y-chromosome diversity is an order of magnitude lower than the mitochondrial diversity, a pattern that we found is best explained by female migration between communities. Thus, the genetic data presented here provide a detailed documentation of the social organization of an isolated Neanderthal community at the easternmost extent of their known range.

The long non-coding RNAs (lncRNAs) other than rRNA and tRNA were earlier assumed to be ‘junk genomic material’. However, recent advancements in genomics methods have highlighted their roles not only in housekeeping but also in the progression of diseases like cancer as well as viral infections. lncRNAs owing to their length, have both short-range and long-range interactions resulting in complex folded structures that recruit various biomolecules enabling lncRNAs to undertake their various biological functions. Using cell lysate pull-down assays increasing number of lnRNAs-interacting proteins are being identified. These interactions can be further exploited to develop targeted novel therapeutic strategies to inhibit lncRNA–protein interactions. This review attempts to succinctly techniques that can identify and characterize the lnRNAs–protein interactions (i.e. affinity, stoichiometry, and thermodynamics). Furthermore, using other sophisticated biophysical techniques, one can also perform size estimations, and determine low-resolution structures. Since these methods study the biomolecules in solution, large-scale structural observations can be performed in real-time. This review attempts to briefly introduce the readers to biochemical and biophysical techniques, such that they can utilize these methods to obtain a holistic characterization of the biomolecules of interest. Additionally, it should be noted that the use of these methods is not limited to the characterization of the interacting molecules but can also be used to determine the efficacy of the therapeutic molecules to disrupt these interactions.

Background: Desmoid tumor, also known as aggressive fibromatosis, is well-characterized by abnormal Wnt/β-catenin signaling. Various therapeutic options, including imatinib, are available to efficiently treat desmoid tumor. However, molecular mechanism of why imatinib works remains poorly understood. Here, we describe the potential roles of NOTCH2 and HES1 in association with clinical response to imatinib as in genome and transcriptome levels. Methods: We identified all somatic mutations in coding and non-coding regions via whole genome sequencing using desmoid tumor samples. To validate the genetic interaction with expression level in desmoid-tumor condition, we utilized large-scale Whole-genome sequencing (WGS) and transcriptome datasets from the Pan-Cancer Analysis of Whole Genomes (PCAWG) project. RNA-sequencing was performed using prospective and retrospective cohort samples to evaluate the expressional relevance with clinical response. Results: Among 20 patients, 4 (20%) had a partial response and 14 (66.7%) had stable disease, 11 of which continued for ≥1 year. With gene-wise functional analyses, we detected significant correlation between recurrent NOTCH2 noncoding mutations and clinical response to imatinib. Based on PCAWG data analyses, NOTCH2 mutations affect its expression levels particularly in the presence of CTNNB1 missense mutations. By analyzing RNA-sequencing with additional desmoid tumor samples, we found that NOTCH2 expression was significantly correlated with HES1 expression. Interestingly, NOTCH2 had no statistical power to discriminate responders and non-responders. Instead, HES1 was differentially expressed with statistical significance between responders and non-responders Conclusions: Imatinib was effective and well tolerated for advanced desmoid tumor treatment. Our results show that HES1, regulated by NOTCH2, as an indicator of sensitivity to imatinib, which provides an important therapeutic consideration for desmoid tumor. Trial Registration: ClinicalTrials.gov, NCT02495519, Registered July 13, 2015- Retrospectively registered, https://www.clinicaltrials.gov/ct2/show/NCT02495519

Emerging evidence shows that some Lactobacillus fermentum strains can contribute to the prevention and treatment of ulcerative colitis (UC). In this study, 105 isolates of L. fermentum strains were separated from fecal samples of populations in different regions in China and their draft genomes were sequenced. Pan-genomic and phylogenetic characterizations of these strains and four model strains (L. fermentum 3872, CECT5716, IFO3956 and VRI003) were performed. Phylogenetic analysis indicated that there was no significant adaptive evolution between the genomes of L. fermentum strains and the geographical location, sex, ethnicity and age of the hosts. Three L. fermentum strains (FWXBH115, FGDLZR121, and FXJCJ61) from different branches of the phylogenetic tree and strain type L. fermentum CECT5716 were selected and their anti-inflammatory and immune modulatory activities in a dextran sulphate sodium (DSS)-induced colitis mouse model were further investigated. Both L. fermentum FXJCJ61 and CECT5716 significantly alleviated UC by reducing all colitis-associated histological indices, maintaining mucosal integrity, and stimulating replenishment of short-chain fatty acids (SCFAs), while the other two strains failed to offer similar protection. The anti-inflammatory mechanisms of L. fermentum FXJCJ61 and CECT5716 were related to the inhibition of nuclear factor kappa-B (NF-κB) signaling pathway activation and enhancement of interleukin 10 (IL-10) production. Comparative genomic analysis of these strains identified candidate genes that may contribute to the anti-inflammatory effects of specific L. fermentum strains.

Sequencing cancer predisposing genes (CPGs) in evocative patients (i.e., patients with personal and family history of multiple/early-onset/unusual cancers) allows follow-up in their relatives to be adapted when a causative pathogenic variant is identified. Unfortunately, many evocative families remain unexplained. Part of this “missing heritability” could be due to CPG dysregulations caused by remote noncoding genomic alterations. Transcription levels are regulated through the ability of promoters to physically interact with their distant cis-regulatory elements. Three-dimensional chromatin contacts, mediated by a dynamic loop extrusion process, are uncovered by chromosome conformation capture (3C) and 3C-derived techniques, which have enabled the discovery of new pathological mechanisms in developmental diseases and cancers. High-penetrance cancer predisposition is caused by germline hereditary alterations otherwise found at the somatic level in sporadic cancers. Thus, data from both developmental diseases and cancers provide information about possible unknown cancer predisposition mechanisms. This mini-review aims to deduce from these data whether abnormal chromatin folding can cause high-penetrance cancer predisposition.

Vertebrate blood coagulation is controlled by a cascade containing more than 20 proteins. The cascade proteins are found in the blood in their zymogen forms and when the cascade is triggered by tissue damage, zymogens are activated and in turn activate their downstream proteins by serine protease activity. In this study, we examined proteomes of 21 chordates, of which 18 are vertebrates, to reveal the modular evolution of the blood coagulation cascade. Additionally, two Arthropoda species were used to compare domain arrangements of the proteins belonging to the hemolymph clotting and the blood coagulation cascades. Within the vertebrate coagulation protein set, almost half of the studied proteins are shared with jawless vertebrates. Domain similarity analyses revealed that there are multiple possible evolutionary trajectories for each coagulation protein. During the evolution of higher vertebrate clades, gene and genome duplications led to the formation of other coagulation cascade proteins.

Cytochrome P450 is a superfamily of enzymes that contain heme as a cofactor that functions as monooxygenases. These proteins are responsible for the oxidation of fatty acids, steroids, and xenobiotics, and the synthesis and breakdown of various hormones. Candida tropicalis has been identified as the leading cause of Candidemia, a bloodstream infection that contributes to the high mortality rate worldwide. So far, no research has reported a detailed study of the cytochrome P450 genes in the genome of Candida tropicalis. Hence, a detailed analysis involving the identification and characterisation of cytochrome P450 genes from the Candida tropicalis genome would provide valuable information about the evolutionary history and functions of the Cytochrome P450 gene family. In this study, 32 Cytochrome P450 genes were identified in Candida tropicalis, which were further grouped into four clades. The results show that most cytochrome P450 families in Candida tropicalis have a close phylogenetic relationship. And diversity might have mainly been a result of gene duplication. Clans and family groupings, alongside their putative functions and subcellular localisation, also showed the possible relationship between the genes’ cellular location and their functions. The results obtained from this study pave the way for a broader and more detailed analysis of the specific functions of each of the Cytochrome P450 genes in Candida tropicalis, as well as possible ways of altering their functions to find a way of curbing/eradicating the negative effect of Candida tropicalis and its contribution to the high rate of mortality worldwide.

The Splendidus clade is the largest clade in Vibrionaceae, and its members are often related to mortality of marine animals with huge economic losses. The molecular bases of their pathogenicity and virulence, however, remain largely unknown. In particular, the complete genome sequences of the Splendidus clade species are rarely registered, which is one of the obstacles to predict core and/or unique genes responsible for their adaptation and pathogenicity, and to perform a fine scale meta-transcriptome during bacterial infection to their hosts. In this study, we obtained the complete genomes of all type strains in the Splendidus clade and revealed that (1) different genome sizes (4.4-5.9 Mb) with V. lentus the biggest and most of them had several big plasmids, likely because of the different features on mobilome elements; (2) the Splendidus clade consists of 19 species except V. cortegadensis, and 3 sub-clades (SC) were defined with the 15 most closely related members as SC1; (3) different carbohydrate degradation preferences may be the result of environmental adaptation; and (4) a broad prediction of virulence factors (VFs) revealed core and species unique VF genes.

The sequence assembly algorithms have rapidly evolved with the vigorous growth of genome sequencing technology over the past two decades. Assembly mainly uses the iterative expansion of overlap relationships between sequences to construct the target genome. The assembly algorithms can be typically classified into several categories, such as the Greedy strategy, Overlap-Layout-Consensus (OLC) strategy, and de Bruijn graph (DBG) strategy. In particular, due to the rapid development of third-generation sequencing (TGS) technology, some prevalent assembly algorithms have been proposed to generate high-quality chromosome-level assemblies. However, due to the genome complexity, the length of short reads, and the high error rate of long reads, contigs produced by assembly may contain misassemblies adversely affecting downstream data analysis. Therefore, several read-based and reference-based methods for misassembly identification have been developed to improve assembly quality. This work primarily reviewed the development of DNA sequencing technologies and summarized sequencing data simulation methods, sequencing error correction methods, various mainstream sequence assembly algorithms, and misassembly identification methods. A large amount of computation makes the sequence assembly problem more challenging, and therefore, it is necessary to develop more efficient and accurate assembly algorithms and alternative algorithms.

Complex structural variants (CSVs) encompass multiple breakpoints and are often missed or misinterpreted by state-of-the-art variant detection algorithms. We developed SVision, a deep-learning based multi-object recognition framework, to automatically detect and characterize CSVs from long-read data. SVision outperforms current variant callers at identifying internal structure of complex events and revealed 80 high-quality CSVs with 25 distinct structures from an individual genome. SVision directly detects CSVs without pattern matching against a database of known structures, allowing sensitive detection of both common and previously uncharacterized complex rearrangements.

The dynamic three-dimensional structures of chromatin and extrachromosomal DNA molecules regulate fundamental cellular processes and beyond. However, the visualization of specific DNA sequences in live cells, especially nonrepetitive sequences accounting for most of the genome, is still vastly challenging. Here, we introduce a robust CRISPR-mediated fluorescence in situ hybridization amplifier (CRISPR FISHer) system, which exploits engineered sgRNA and protein trimerization domain-mediated, phase separation-based exponential assembly of fluorescent proteins in the CRISPR-targeting locus, conferring enhancements in both local brightness and signal-to-background ratio and thus achieving single sgRNA-directed visualization of native nonrepetitive DNA loci in live cells. In one application, by labeling and tracking the broken ends of chromosomal fragments, CRISPR FISHer enables real-time visualization of the entire process of chromosome breakage, separation, and subsequent intra- or inter-chromosomal ends rejoining in a single live cell. Furthermore, CRISPR FISHer allows the movement of small extrachromosomal circular DNAs (eccDNAs) and invading DNAs to be recorded, revealing substantial differences in dynamic behaviors between chromosomal and extrachromosomal loci. With the potential to track any specified self or non-self DNA sequences, CRISPR FISHer dramatically broadens the scope of live-cell imaging in biological events and for biomedical diagnoses.

Inbreeding and loss of genetic diversity are serious problems in local Chinese pig breeds. Runs of homozygosity (ROHs) are contiguous lengths of homozygous genotypes that can provide information on inbreeding levels, mating schemes, selection pressure, and genetic events. The Anqing six-end-white (AQ) pig, an important autochthonous pig breed bred in the Anhui province of China, plays a key role in the local pig industry. In recent decades, AQ pigs have become a closed breed with a small population size in conservation farms, raising the issue of inbreeding decline. In this study, we used 10× resequencing to detect the whole genome of 24 AQ pigs and six Asian wild boars (AWBs). We used the Plink software to analyze the homozygosity of the genome and the distribution of ROHs in the genome based on the resequencing data. We obtained 935.04 Gb of raw data from the resequencing results and more than 45 822 239 SNPs. Additionally, we identified 28 702 ROHs. Compared with AWB, AQ pigs had higher ROH numbers, longer ROH rates, and higher F