Meiotic recombination mirrors patterns of germline replication in mice and humans.

Genetic recombination generates novel trait combinations, and understanding how recombination is distributed across the genome is key to modern genetics. The PRDM9 protein defines recombination hotspots; however, megabase-scale recombination patterning is independent of PRDM9. The single round of DNA replication, which precedes recombination in meiosis, may establish these patterns; therefore, we devised an approach to study meiotic replication that includes robust and sensitive mapping of replication origins. We find that meiotic DNA replication is distinct; reduced origin firing slows replication in meiosis, and a distinctive replication pattern in human males underlies the subtelomeric increase in recombination. We detected a robust correlation between replication and both contemporary and historical recombination and found that replication origin density coupled with chromosome size determines the recombination potential of individual chromosomes. Our findings and methods have implications for understanding the mechanisms underlying DNA replication, genetic recombination, and the landscape of mammalian germline variation.

Isolation, characterization, and application of a novel Pseudomonas fluorescens phage vB_PF_Y1-MI in contaminated milk.

The food industry has incurred substantial losses from contamination by Pseudomonas fluorescens, emphasizing the critical importance of implementing effective control strategies. Phages are potential sterilizers due to their specific killing abilities and the difficulty bacteria face in developing resistance. However, a significant barrier to their development is the lack of diversity among phage types. In this study, we characterized a novel lytic P. fluorescens phage, named vB_PF_Y1-MI. Phage vB_PF_Y1-MI displayed a latent period of nearly 10 min and a high burst size of 1493 PFU/cell. This phage showed good activity over a wide range of temperature (up to 70 °C) and pH (3-12). The genome of phage vB_PF_Y1-MI spans 93,233 bp with a GC content of 45%. It encompasses 174 open-reading frames and 19 tRNA genes, while no lysogeny or virulence-associated genes were detected. Phylogenetic analysis positions it as a novel unassigned evolutionary lineage within the Caudoviricetes class among related dsDNA phages. Our study provides foundational insights into vB_PF_Y1-MI and emphasizes its potential as an effective biological control agent against P. fluorescens. This research offers crucial theoretical groundwork and technical support for subsequent efforts in preventing and controlling P. fluorescens contamination.

Probiotic functional gene explorations in the genome of Limosilactobacillus fermentum GD5MG.

Limosilactobacillus fermentum is an isolate obtained from oral gingival samples of healthy human individuals. The whole genome of Lb. fermentum GD5MG is composed of a circular DNA molecule containing 1,834,134 bp and exhibits a GC content of 52.80 %. The sequencing effort produced 38.6 million reads, each 150 bp in length, resulting in a sequencing depth of 2912.48x. Our examination unveiled a total of 1961 protein-coding genes, 27 rRNA genes, 24 tRNA genes, 3 non-coding RNA genes, and 63 pseudogenes with the use of gene annotations in NCBI Prokaryotic Genome Annotation tool. RAST revealed 1863 coding genes distributed across 209 subsystems, with a predominant involvement in amino acid, carbohydrate, and protein metabolism. Phylogenetic analysis infers that the Lb. fermentum GD5MG shares 281 gene clusters. Furthermore, the genome features showed a single CRISPR locus of 45 bp in length. Three genes associated with adhesion ability (strA, dltD, and dltA) and 26 genes related to acid tolerance, digestive enzyme secretion, and bile salt resistance were identified. Numerous genes associated with oral probiotic properties, comprising adhesion, acid and bile salt tolerance, oxidative stress tolerance, and sugar metabolism, were identified in the genome. Our findings shed light on the genomic characteristics of Lb. fermentum GD5MG, which are probable probiotics with functional benefits in humans.

Comparative mitochondrial genome analysis provides new insights into the classification of Modiolinae.

BACKGROUND: Mitochondrial genomes have become a powerful tool for studying molecular genetics and phylogeny of mollusks. Currently, the position of Modiolinae within Mytilidae and the taxonomic and phylogenetic relationships within Modiolinae were still controversial. This study focuses on the complete mitochondrial genomes of two species: Modiolus modulaides (Röding, 1798) and Modiolus auriculatus Krauss, 1848, which have not been sequenced before. METHODS AND RESULTS: We assembled and characterized the mitochondrial genomes of M. modulaides and M. auriculatus and then analyzed the phylogenetic relationships. The mitochondrial genomes of M. modulaides and M. auriculatus were 15,422 bp and 16,027 bp, respectively. Both of them were composed of 36 functional genes, including 12 protein-coding genes, 22 transfer RNAs, and 2 ribosomal RNAs. All protein-coding genes showed A + T bias, positive GC skews, and negative AT skews in nucleotide composition. Phylogenetic analysis based on the mitochondrial genomes showed that Modiolinae and Bathymodiolinae clustered together to form a sister relationship. Seven Modiolinae species were divided into two clades: L1 (M. modulaides, M. auriculatus and Modiolus philippinarum Hanley, 1843) and L2 [Modiolus modiolus (Linnaeus, 1758), Modiolus kurilensis Bernard, 1983, Modiolus nipponicus (Oyama, 1950), and Modiolus comptus (Sowerby III, 1915)]. The divergence time of the two clades was approximately 105.75 Ma. Furthermore, the transfer RNA gene rearrangement, longer genetic distance, and greater genetic differentiation were confirmed between the L1 and L2 clades, as well as differences in the external characteristics of the shells of the two clades. CONCLUSIONS: Based on the molecular data, it was speculated that species from the L1 clade might belong to other genera or new genera. This study provides molecular information for further taxonomic and phylogenetic studies of Mytilidae.

Complete mitochondrial genome and phylogenetic analysis of Mancinella alouina.

BACKGROUND: The Muricidae family in the Class Gastropoda comprises numerous species with a vast range of morphological features and a worldwide presence. The phylogeny of the Muricidae has been analyzed in previous studies; however, the evolutionary relationships among the main branches of the Muricidae remain unknown. METHODS AND RESULTS: In the present study, the mitochondrial genome of Mancinella alouina was sequenced. The mitochondrial genome was found to be 16,671 bp in length and made up of 37 genes (13 protein-coding genes, 22 transfer RNA and 2 ribosomal RNA genes). The genome has an A-T-rich region (66.5% A + T content) and all of the PCGs use the ATN start codon and the TAG or TAA stop codons. The mitochondrial gene arrangement of Mancinella alouina is similar to that of other Muricidae, except for Ocinebrellus inornatus and Ceratostoma burnetti. On the basis of a flexible molecular clock model, time-calibrated phylogenetic results indicate that the genus Mancinella diverged roughly 18.09 Mya, and that the family Muricidae emerged in the Late Cretaceous. CONCLUSIONS: This study reveals the structural and sequence information features of the mitochondrial genome of Mancinella alouina. This study provides evidence for the relationships within the family Muricidae at the molecular level, and infer the divergence time. The results of phylogenetic analyses strongly support the current classification.

Genomic analysis and phylogenetic characterization of Himalayan snow trout, Schizothorax esocinus based on mitochondrial protein-coding genes.

BACKGROUND: Mitochondrial DNA (mtDNA) has become a significant tool for exploring genetic diversity and delineating evolutionary links across diverse taxa. Within the group of cold-water fish species that are native to the Indian Himalayan region, Schizothorax esocinus holds particular importance due to its ecological significance and is potentially vulnerable to environmental changes. This research aims to clarify the phylogenetic relationships within the Schizothorax genus by utilizing mitochondrial protein-coding genes. METHODS: Standard protocols were followed for the isolation of DNA from S. esocinus. For the amplification of mtDNA, overlapping primers were used, and then subsequent sequencing was performed. The genetic features were investigated by the application of bioinformatic approaches. These approaches covered the evaluation of nucleotide composition, codon usage, selective pressure using nonsynonymous substitution /synonymous substitution (Ka/Ks) ratios, and phylogenetic analysis. RESULTS: The study specifically examined the 13 protein-coding genes of Schizothorax species which belongs to the Schizothoracinae subfamily. Nucleotide composition analysis showed a bias towards A + T content, consistent with other cyprinid fish species, suggesting evolutionary conservation. Relative Synonymous Codon Usage highlighted leucine as the most frequent (5.18%) and cysteine as the least frequent (0.78%) codon. The positive AT-skew and the predominantly negative GC-skew indicated the abundance of A and C. Comparative analysis revealed significant conservation of amino acids in multiple genes. The majority of amino acids were hydrophobic rather than polar. The purifying selection was revealed by the genetic distance and Ka/Ks ratios. Phylogenetic study revealed a significant genetic divergence between S. esocinus and other Schizothorax species with interspecific K2P distances ranging from 0.00 to 8.87%, with an average of 5.76%. CONCLUSION: The present study provides significant contributions to the understanding of mitochondrial genome diversity and genetic evolution mechanisms in Schizothoracinae, hence offering vital insights for the development of conservation initiatives aimed at protecting freshwater fish species.

Comprehensive analysis of the complete mitochondrial genome of Lilium tsingtauense reveals a novel multichromosome structure.

KEY MESSAGE: Lilium tsingtauense mitogenome comprises 27 independent chromosome molecules, it undergoes frequent genomic recombination, and the rate of recombination and mutation between different repetitive sequences affects the formation of multichromosomal structures. Given the extremely large genome of Lily, which likely harbors additional genetic resources, it serves as an ideal material for studying the phylogenetic evolution of organisms. Although the Lilium chloroplast genome has been documented, the sequence of its mitochondrial genome (mitogenome) remains uncharted. Using BGI short reads and Nanopore long reads, we sequenced, assembled, and annotated the mitogenome of Lilium tsingtauense. This effort culminated in the characterization of Lilium's first complete mitogenome. Comparative analysis with other angiosperms revealed the unique multichromosomal structure of the L. tsingtauense mitogenome, spanning 1,125,108 bp and comprising 27 independent circular chromosomes. It contains 36 protein-coding genes, 12 tRNA genes, and 3 rRNA genes, with a GC content of 44.90%. Notably, three chromosomes in the L. tsingtauense mitogenome lack identifiable genes, hinting at the potential existence of novel genes and noncoding elements. The high degree of observed genome fragmentation implies frequent reorganization, with recombination and mutation rates among diverse repetitive sequences likely driving the formation of multichromosomal structures. Our comprehensive analysis, covering genome size, coding genes, structure, RNA editing, repetitive sequences, and sequence migration, sheds light on the evolutionary and molecular biology of multichromosomal mitochondria in Lilium. This high-quality mitogenome of L. tsingtauense not only enriches our understanding of multichromosomal mitogenomes but also establishes a solid foundation for future genome breeding and germplasm innovation in Lilium.

Comparative Analysis of Codon Usage Bias in Six Eimeria Genomes.

The codon usage bias (CUB) of genes encoded by different species' genomes varies greatly. The analysis of codon usage patterns enriches our comprehension of genetic and evolutionary characteristics across diverse species. In this study, we performed a genome-wide analysis of CUB and its influencing factors in six sequenced Eimeria species that cause coccidiosis in poultry: Eimeria acervulina, Eimeria necatrix, Eimeria brunetti, Eimeria tenella, Eimeria praecox, and Eimeria maxima. The GC content of protein-coding genes varies between 52.67% and 58.24% among the six Eimeria species. The distribution trend of GC content at different codon positions follows GC1 > GC3 > GC2. Most high-frequency codons tend to end with C/G, except in E. maxima. Additionally, there is a positive correlation between GC3 content and GC3s/C3s, but a significantly negative correlation with A3s. Analysis of the ENC-Plot, neutrality plot, and PR2-bias plot suggests that selection pressure has a stronger influence than mutational pressure on CUB in the six Eimeria genomes. Finally, we identified from 11 to 15 optimal codons, with GCA, CAG, and AGC being the most commonly used optimal codons across these species. This study offers a thorough exploration of the relationships between CUB and selection pressures within the protein-coding genes of Eimeria species. Genetic evolution in these species appears to be influenced by mutations and selection pressures. Additionally, the findings shed light on unique characteristics and evolutionary traits specific to the six Eimeria species.

Decoding the Chloroplast Genome of Tetrastigma (Vitaceae): Variations and Phylogenetic Selection Insights.

Tetrastigma (Vitaceae) is known for its ornamental, medicinal, and ecological significance. However, the structural and variational characteristics of the Tetrastigma chloroplast genome and their impact on phylogenetic relationships remain underexplored. This study utilized bioinformatics methods to assemble and annotate the chloroplast genomes of 10 Tetrastigma species and compare them with five previously sequenced species. This study analyzed gene composition, simple sequence repeats, and codon usage patterns, revealing a high A/T content, uniquely identified pentanucleotide repeats in five species and several preferred codons. In addition, comparative analyses were conducted of the chloroplast genomes of 15 Tetrastigma species, examining their structural differences and identifying polymorphic hotspots (rps16, rps16-trnQ, trnS, trnD, psbC-trnS-psbZ, accD-psaI, psbE-petL-petG, etc.) suitable for DNA marker development. Furthermore, phylogenetic and selective pressure analyses were performed based on the chloroplast genomes of these 15 Tetrastigma species, validating and elucidating intra-genus relationships within Tetrastigma. Futhermore, several genes under positive selection, such as atpF and accD, were identified, shedding light on the adaptive evolution of Tetrastigma. Utilizing 40 Vitaceae species, the divergence time of Tetrastigma was estimated, clarifying the evolutionary relationships within Tetrastigma relative to other genera. The analysis revealed diverse divergences of Tetrastigma in the Miocene and Pliocene, with possible ancient divergence events before the Eocene. Furthermore, family-level selective pressure analysis identified key features distinguishing Tetrastigma from other genera, showing a higher degree of purifying selection. This research enriches the chloroplast genome data for Tetrastigma and offers new insights into species identification, phylogenetic analysis, and adaptive evolution, enhancing our understanding of the genetic diversity and evolutionary history of these species.

Budding yeast Rap1, but not telomeric DNA, is inhibitory for multiple stages of DNA replication in vitro.

Telomeres are copied and reassembled each cell division cycle through a multistep process called telomere replication. Most telomeric DNA is duplicated semiconservatively during this process, but replication forks frequently pause or stall at telomeres in yeast, mouse and human cells, potentially causing chronic telomere shortening or loss in a single cell cycle. We have investigated the cause of this effect by examining the replication of telomeric templates in vitro. Using a reconstituted assay for eukaryotic DNA replication in which a complete eukaryotic replisome is assembled and activated with purified proteins, we show that budding yeast telomeric DNA is efficiently duplicated in vitro unless the telomere binding protein Rap1 is present. Rap1 acts as a roadblock that prevents replisome progression and leading strand synthesis, but also potently inhibits lagging strand telomere replication behind the fork. Both defects can be mitigated by the Pif1 helicase. Our results suggest that GC-rich sequences do not inhibit DNA replication per se, and that in the absence of accessory factors, telomere binding proteins can inhibit multiple, distinct steps in the replication process.

Differences in firing efficiency, chromatin, and transcription underlie the developmental plasticity of the Arabidopsis DNA replication origins.

Eukaryotic genome replication depends on thousands of DNA replication origins (ORIs). A major challenge is to learn ORI biology in multicellular organisms in the context of growing organs to understand their developmental plasticity. We have identified a set of ORIs of Arabidopsis thaliana and their chromatin landscape at two stages of post-embryonic development. ORIs associate with multiple chromatin signatures including transcription start sites (TSS) but also proximal and distal regulatory regions and heterochromatin, where ORIs colocalize with retrotransposons. In addition, quantitative analysis of ORI activity led us to conclude that strong ORIs have high GC content and clusters of GGN trinucleotides. Development primarily influences ORI firing strength rather than ORI location. ORIs that preferentially fire at early developmental stages colocalize with GC-rich heterochromatin, but at later stages with transcribed genes, perhaps as a consequence of changes in chromatin features associated with developmental processes. Our study provides the set of ORIs active in an organism at the post-embryo stage that should allow us to study ORI biology in response to development, environment, and mutations with a quantitative approach. In a wider scope, the computational strategies developed here can be transferred to other eukaryotic systems.

Interplay between coding and exonic splicing regulatory sequences.

The inclusion of exons during the splicing process depends on the binding of splicing factors to short low-complexity regulatory sequences. The relationship between exonic splicing regulatory sequences and coding sequences is still poorly understood. We demonstrate that exons that are coregulated by any given splicing factor share a similar nucleotide composition bias and preferentially code for amino acids with similar physicochemical properties because of the nonrandomness of the genetic code. Indeed, amino acids sharing similar physicochemical properties correspond to codons that have the same nucleotide composition bias. In particular, we uncover that the TRA2A and TRA2B splicing factors that bind to adenine-rich motifs promote the inclusion of adenine-rich exons coding preferentially for hydrophilic amino acids that correspond to adenine-rich codons. SRSF2 that binds guanine/cytosine-rich motifs promotes the inclusion of GC-rich exons coding preferentially for small amino acids, whereas SRSF3 that binds cytosine-rich motifs promotes the inclusion of exons coding preferentially for uncharged amino acids, like serine and threonine that can be phosphorylated. Finally, coregulated exons encoding amino acids with similar physicochemical properties correspond to specific protein features. In conclusion, the regulation of an exon by a splicing factor that relies on the affinity of this factor for specific nucleotide(s) is tightly interconnected with the exon-encoded physicochemical properties. We therefore uncover an unanticipated bidirectional interplay between the splicing regulatory process and its biological functional outcome.

CtIP maintains stability at common fragile sites and inverted repeats by end resection-independent endonuclease activity.

Chromosomal rearrangements often occur at genomic loci with DNA secondary structures, such as common fragile sites (CFSs) and palindromic repeats. We developed assays in mammalian cells that revealed CFS-derived AT-rich sequences and inverted Alu repeats (Alu-IRs) are mitotic recombination hotspots, requiring the repair functions of carboxy-terminal binding protein (CtBP)-interacting protein (CtIP) and the Mre11/Rad50/Nbs1 complex (MRN). We also identified an endonuclease activity of CtIP that is dispensable for end resection and homologous recombination (HR) at I-SceI-generated "clean" double-strand breaks (DSBs) but is required for repair of DSBs occurring at CFS-derived AT-rich sequences. In addition, CtIP nuclease-defective mutants are impaired in Alu-IRs-induced mitotic recombination. These studies suggest that an end resection-independent CtIP function is important for processing DSB ends with secondary structures to promote HR. Furthermore, our studies uncover an important role of MRN, CtIP, and their associated nuclease activities in protecting CFSs in mammalian cells.

GC ends control topology of DNA G-quadruplexes and their cation-dependent assembly.

GCn and GCnCG, where n = (G2AG4AG2), fold into well-defined, dimeric G-quadruplexes with unprecedented folding topologies in the presence of Na+ ions as revealed by nuclear magnetic resonance spectroscopy. Both G-quadruplexes exhibit unique combination of structural elements among which are two G-quartets, A(GGGG)A hexad and GCGC-quartet. Detailed structural characterization uncovered the crucial role of 5'-GC ends in formation of GCn and GCnCG G-quadruplexes. Folding in the presence of 15NH4+ and K+ ions leads to 3'-3' stacking of terminal G-quartets of GCn G-quadruplexes, while 3'-GC overhangs in GCnCG prevent dimerization. Results of the present study expand repertoire of possible G-quadruplex structures. This knowledge will be useful in DNA sequence design for nanotechnological applications that may require specific folding topology and multimerization properties.

Single-molecule DNA sequencing of widely varying GC-content using nucleotide release, capture and detection in microdroplets.

Despite remarkable progress in DNA sequencing technologies there remains a trade-off between short-read platforms, having limited ability to sequence homopolymers, repeated motifs or long-range structural variation, and long-read platforms, which tend to have lower accuracy and/or throughput. Moreover, current methods do not allow direct readout of epigenetic modifications from a single read. With the aim of addressing these limitations, we have developed an optical electrowetting sequencing platform that uses step-wise nucleotide triphosphate (dNTP) release, capture and detection in microdroplets from single DNA molecules. Each microdroplet serves as a reaction vessel that identifies an individual dNTP based on a robust fluorescence signal, with the detection chemistry extended to enable detection of 5-methylcytosine. Our platform uses small reagent volumes and inexpensive equipment, paving the way to cost-effective single-molecule DNA sequencing, capable of handling widely varying GC-bias, and demonstrating direct detection of epigenetic modifications.

iFLinkC: an iterative functional linker cloning strategy for the combinatorial assembly and recombination of linker peptides with functional domains.

Recent years have witnessed increasing efforts to engineer artificial biological functions through recombination of modular-organized toolboxes of protein scaffolds and parts. A critical, yet frequently neglected aspect concerns the identity of peptide linkers or spacers connecting individual domains which remain poorly understood and challenging to assemble. Addressing these limitations, iFlinkC comprises a highly scalable DNA assembly process that facilitates the combinatorial recombination of functional domains with linkers of varying length and flexibility, thereby overcoming challenges with high GC-content and the repeat nature of linker elements. The capacity of iFLinkC is demonstrated in the construction of synthetic protease switches featuring PDZ-FN3-based affinity clamps and single-chain FKBP12-FRB receptors as allosteric inputs. Library screening experiments demonstrate that linker space is highly plastic as the induction of allosterically regulated protease switches can vary from >150-fold switch-ON to >13-fold switch-OFF solely depending on the identity of the connecting linkers and relative orientation of functional domains. In addition, Pro-rich linkers yield the most potent switches contradicting the conventional use of flexible Gly-Ser linkers. Given the ease and efficiency how functional domains can be readily recombined with any type of linker, iFLinkC is anticipated to be widely applicable to the assembly of any type of fusion protein.

Replisome bypass of transcription complexes and R-loops.

The vast majority of the genome is transcribed by RNA polymerases. G+C-rich regions of the chromosomes and negative superhelicity can promote the invasion of the DNA by RNA to form R-loops, which have been shown to block DNA replication and promote genome instability. However, it is unclear whether the R-loops themselves are sufficient to cause this instability or if additional factors are required. We have investigated replisome collisions with transcription complexes and R-loops using a reconstituted bacterial DNA replication system. RNA polymerase transcription complexes co-directionally oriented with the replication fork were transient blockages, whereas those oriented head-on were severe, stable blockages. On the other hand, replisomes easily bypassed R-loops on either template strand. Replication encounters with R-loops on the leading-strand template (co-directional) resulted in gaps in the nascent leading strand, whereas lagging-strand template R-loops (head-on) had little impact on replication fork progression. We conclude that whereas R-loops alone can act as transient replication blocks, most genome-destabilizing replication fork stalling likely occurs because of proteins bound to the R-loops.

Linking high GC content to the repair of double strand breaks in prokaryotic genomes.

Genomic GC content varies widely among microbes for reasons unknown. While mutation bias partially explains this variation, prokaryotes near-universally have a higher GC content than predicted solely by this bias. Debate surrounds the relative importance of the remaining explanations of selection versus biased gene conversion favoring GC alleles. Some environments (e.g. soils) are associated with a high genomic GC content of their inhabitants, which implies that either high GC content is a selective adaptation to particular habitats, or that certain habitats favor increased rates of gene conversion. Here, we report a novel association between the presence of the non-homologous end joining DNA double-strand break repair pathway and GC content; this observation suggests that DNA damage may be a fundamental driver of GC content, leading in part to the many environmental patterns observed to-date. We discuss potential mechanisms accounting for the observed association, and provide preliminary evidence that sites experiencing higher rates of double-strand breaks are under selection for increased GC content relative to the genomic background.

Comparative analysis of human coronaviruses focusing on nucleotide variability and synonymous codon usage patterns.

The prevailing COVID-19 pandemic has drawn the attention of the scientific community to study the evolutionary origin of Severe Acute Respiratory Syndrome Corona Virus 2 (SARS-CoV-2). This study is a comprehensive quantitative analysis of the protein-coding sequences of seven human coronaviruses (HCoVs) to decipher the nucleotide sequence variability and codon usage patterns. It is essential to understand the survival ability of the viruses, their adaptation to hosts, and their evolution. The current analysis revealed a high abundance of the relative dinucleotide (odds ratio), GC and CT pairs in the first and last two codon positions, respectively, as well as a low abundance of the CG pair in the last two positions of the codon, which might be related to the evolution of the viruses. A remarkable level of variability of GC content in the third position of the codon among the seven coronaviruses was observed. Codons with high RSCU values are primarily from the aliphatic and hydroxyl amino acid groups, and codons with low RSCU values belong to the aliphatic, cyclic, positively charged, and sulfur-containing amino acid groups. In order to elucidate the evolutionary processes of the seven coronaviruses, a phylogenetic tree (dendrogram) was constructed based on the RSCU scores of the codons. The severe and mild categories CoVs were positioned in different clades. A comparative phylogenetic study with other coronaviruses depicted that SARS-CoV-2 is close to the CoV isolated from pangolins (Manis javanica, Pangolin-CoV) and cats (Felis catus, SARS(r)-CoV). Further analysis of the effective number of codon (ENC) usage bias showed a relatively higher bias for SARS-CoV and MERS-CoV compared to SARS-CoV-2. The ENC plot against GC3 suggested that the mutational bias might have a role in determining the codon usage variation among candidate viruses. A codon adaptability study on a few human host parasites (from different kingdoms), including CoVs, showed a diverse adaptability pattern. SARS-CoV-2 and SARS-CoV exhibit relatively lower but similar codon adaptability compared to MERS-CoV.

Genome-Wide Prediction of Transcription Start Sites in Conifers.

The identification of promoters is an essential step in the genome annotation process, providing a framework for gene regulatory networks and their role in transcription regulation. Despite considerable advances in the high-throughput determination of transcription start sites (TSSs) and transcription factor binding sites (TFBSs), experimental methods are still time-consuming and expensive. Instead, several computational approaches have been developed to provide fast and reliable means for predicting the location of TSSs and regulatory motifs on a genome-wide scale. Numerous studies have been carried out on the regulatory elements of mammalian genomes, but plant promoters, especially in gymnosperms, have been left out of the limelight and, therefore, have been poorly investigated. The aim of this study was to enhance and expand the existing genome annotations using computational approaches for genome-wide prediction of TSSs in the four conifer species: loblolly pine, white spruce, Norway spruce, and Siberian larch. Our pipeline will be useful for TSS predictions in other genomes, especially for draft assemblies, where reliable TSS predictions are not usually available. We also explored some of the features of the nucleotide composition of the predicted promoters and compared the GC properties of conifer genes with model monocot and dicot plants. Here, we demonstrate that even incomplete genome assemblies and partial annotations can be a reliable starting point for TSS annotation. The results of the TSS prediction in four conifer species have been deposited in the Persephone genome browser, which allows smooth visualization and is optimized for large data sets. This work provides the initial basis for future experimental validation and the study of the regulatory regions to understand gene regulation in gymnosperms.