A defective viral genome strategy elicits broad protective immunity against respiratory viruses.

RNA viruses generate defective viral genomes (DVGs) that can interfere with replication of the parental wild-type virus. To examine their therapeutic potential, we created a DVG by deleting the capsid-coding region of poliovirus. Strikingly, intraperitoneal or intranasal administration of this genome, which we termed eTIP1, elicits an antiviral response, inhibits replication, and protects mice from several RNA viruses, including enteroviruses, influenza, and SARS-CoV-2. While eTIP1 replication following intranasal administration is limited to the nasal cavity, its antiviral action extends non-cell-autonomously to the lungs. eTIP1 broad-spectrum antiviral effects are mediated by both local and distal type I interferon responses. Importantly, while a single eTIP1 dose protects animals from SARS-CoV-2 infection, it also stimulates production of SARS-CoV-2 neutralizing antibodies that afford long-lasting protection from SARS-CoV-2 reinfection. Thus, eTIP1 is a safe and effective broad-spectrum antiviral generating short- and long-term protection against SARS-CoV-2 and other respiratory infections in animal models.

NAD metabolism modulates inflammation and mitochondria function in diabetic kidney disease.

Diabetes mellitus is the leading cause of cardiovascular and renal disease in the United -States. Despite the beneficial interventions available for patients with diabetes, there remains a need for additional therapeutic targets and therapies in diabetic kidney disease (DKD). Inflammation and oxidative stress are increasingly recognized as important causes of renal diseases. Inflammation is closely associated with mitochondrial damage. The molecular connection between inflammation and mitochondrial metabolism remains to be elucidated. Recently, nicotinamide adenine nucleotide (NAD+) metabolism has been found to regulate immune function and inflammation. In the present studies, we tested the hypothesis that enhancing NAD metabolism could prevent inflammation in and progression of DKD. We found that treatment of db/db mice with type 2 diabetes with nicotinamide riboside (NR) prevented several manifestations of kidney dysfunction (i.e., albuminuria, increased urinary kidney injury marker-1 (KIM1) excretion, and pathologic changes). These effects were associated with decreased inflammation, at least in part via inhibiting the activation of the cyclic GMP-AMP synthase-stimulator of interferon genes (cGAS-STING) signaling pathway. An antagonist of the serum stimulator of interferon genes (STING) and whole-body STING deletion in diabetic mice showed similar renoprotection. Further analysis found that NR increased SIRT3 activity and improved mitochondrial function, which led to decreased mitochondrial DNA damage, a trigger for mitochondrial DNA leakage which activates the cGAS-STING pathway. Overall, these data show that NR supplementation boosted NAD metabolism to augment mitochondrial function, reducing inflammation and thereby preventing the progression of diabetic kidney disease.

Estrogen-Related Receptor Agonism Reverses Mitochondrial Dysfunction and Inflammation in the Aging Kidney.

A gradual decline in renal function occurs even in healthy aging individuals. In addition to aging, per se, concurrent metabolic syndrome and hypertension, which are common in the aging population, can induce mitochondrial dysfunction and inflammation, which collectively contribute to age-related kidney dysfunction and disease. This study examined the role of the nuclear hormone receptors, the estrogen-related receptors (ERRs), in regulation of age-related mitochondrial dysfunction and inflammation. The ERRs were decreased in both aging human and mouse kidneys and were preserved in aging mice with lifelong caloric restriction (CR). A pan-ERR agonist, SLU-PP-332, was used to treat 21-month-old mice for 8 weeks. In addition, 21-month-old mice were treated with a stimulator of interferon genes (STING) inhibitor, C-176, for 3 weeks. Remarkably, similar to CR, an 8-week treatment with a pan-ERR agonist reversed the age-related increases in albuminuria, podocyte loss, mitochondrial dysfunction, and inflammatory cytokines, via the cyclic GMP-AMP synthase-STING and STAT3 signaling pathways. A 3-week treatment of 21-month-old mice with a STING inhibitor reversed the increases in inflammatory cytokines and the senescence marker, p21/cyclin dependent kinase inhibitor 1A (Cdkn1a), but also unexpectedly reversed the age-related decreases in PPARG coactivator (PGC)-1α, ERRα, mitochondrial complexes, and medium chain acyl coenzyme A dehydrogenase (MCAD) expression. These studies identified ERRs as CR mimetics and as important modulators of age-related mitochondrial dysfunction and inflammation. These findings highlight novel druggable pathways that can be further evaluated to prevent progression of age-related kidney disease.

Transcriptome landscapes that signify Botrylloides leachi (Ascidiacea) torpor states.

Harsh environments enforce the expression of behavioural, morphological, physiological, and reproductive rejoinders, including torpor. Here we study the morphological, cellular, and molecular alterations in torpor architype in the colonial urochordate Botrylloides aff. leachii by employing whole organism Transmission electron (TEM) and light microscope observations, RNA sequencing, real-time polymerase chain reaction (qPCR) quantification of selected genes, and immunolocalization of WNT, SMAD and SOX2 gene expressions. On the morphological level, torpor starts with gradual regression of all zooids and buds which leaves the colony surviving as condensed vasculature remnants that may be 'aroused' to regenerate fully functional colonies upon changes in the environment. Simultaneously, we observed altered distributions of hemolymph cell types. Phagocytes doubled in number, while the number of morula cells declined by half. In addition, two new circulating cell types were observed, multi-nucleated and bacteria-bearing cells. RNA sequencing technology revealed marked differences in gene expression between different organism compartments and states: active zooids and ampullae, and between mid-torpor and naive colonies, or naive and torpid colonies. Gene Ontology term enrichment analyses further showed disparate biological processes. In torpid colonies, we observed overall 233 up regulated genes. These genes included NR4A2, EGR1, MUC5AC, HMCN2 and. Also, 27 transcription factors were upregulated in torpid colonies including ELK1, HDAC3, RBMX, MAZ, STAT1, STAT4 and STAT6. Interestingly, genes involved in developmental processes such as SPIRE1, RHOA, SOX11, WNT5A and SNX18 were also upregulated in torpid colonies. We further validated the dysregulation of 22 genes during torpor by utilizing qPCR. Immunohistochemistry of representative genes from three signaling pathways revealed high expression of these genes in circulated cells along torpor. WNT agonist administration resulted in early arousal from torpor in 80% of the torpid colonies while in active colonies WNT agonist triggered the torpor state. Abovementioned results thus connote unique transcriptome landscapes associated with Botrylloides leachii torpor.

Genome evolution of blind subterranean mole rats: Adaptive peripatric versus sympatric speciation.

Speciation mechanisms remain controversial. Two speciation models occur in Israeli subterranean mole rats, genus Spalax: a regional speciation cline southward of four peripatric climatic chromosomal species and a local, geologic-edaphic, genic, and sympatric speciation. Here we highlight their genome evolution. The five species were separated into five genetic clusters by single nucleotide polymorphisms, copy number variations (CNVs), repeatome, and methylome in sympatry. The regional interspecific divergence correspond to Pleistocene climatic cycles. Climate warmings caused chromosomal speciation. Triple effective population size, Ne , declines match glacial cold cycles. Adaptive genes evolved under positive selection to underground stresses and to divergent climates, involving interspecies reproductive isolation. Genomic islands evolved mainly due to adaptive evolution involving ancient polymorphisms. Repeatome, including both CNV and LINE1 repetitive elements, separated the five species. Methylation in sympatry identified geologically chalk-basalt species that differentially affect thermoregulation, hypoxia, DNA repair, P53, and other pathways. Genome adaptive evolution highlights climatic and geologic-edaphic stress evolution and the two speciation models, peripatric and sympatric.

Intra- and Inter-cellular Modeling of Dynamic Interaction between Zika Virus and Its Naturally Occurring Defective Viral Genomes.

Here, we examine in silico the infection dynamics and interactions of two Zika virus (ZIKV) genomes: one is the full-length ZIKV genome (wild type [WT]), and the other is one of the naturally occurring defective viral genomes (DVGs), which can replicate in the presence of the WT genome, appears under high-MOI (multiplicity of infection) passaging conditions, and carries a deletion encompassing part of the structural and NS1 protein-coding region. Ordinary differential equations (ODEs) were used to simulate the infection of cells by virus particles and the intracellular replication of the WT and DVG genomes that produce these particles. For each virus passage in Vero and C6/36 cell cultures, the rates of the simulated processes were fitted to two types of observations: virus titer data and the assembled haplotypes of the replicate passage samples. We studied the consistency of the model with the experimental data across all passages of infection in each cell type separately as well as the sensitivity of the model's parameters. We also determined which simulated processes of virus evolution are the most important for the adaptation of the WT and DVG interplay in these two disparate cell culture environments. Our results demonstrate that in the majority of passages, the rates of DVG production are higher inC6/36 cells than in Vero cells, which might result in tolerance and therefore drive the persistence of the mosquito vector in the context of ZIKV infection. Additionally, the model simulations showed a slower accumulation of infected cells under higher activation of the DVG-associated processes, which indicates a potential role of DVGs in virus attenuation. IMPORTANCE One of the ideas for lessening Zika pathogenicity is the addition of its natural or engineered defective virus genomes (DVGs) (have no pathogenicity) to the infection pool: a DVG is redirecting the wild-type (WT)-associated virus development resources toward its own maturation. The mathematical model presented here, attuned to the data from interplays between WT Zika viruses and their natural DVGs in mammalian and mosquito cells, provides evidence that the loss of uninfected cells is attenuated by the DVG development processes. This model enabled us to estimate the rates of virus development processes in the WT/DVG interplay, determine the key processes, and show that the key processes are faster in mosquito cells than in mammalian ones. In general, the presented model and its detailed study suggest in what important virus development processes the therapeutically efficient DVG might compete with the WT; this may help in assembling engineered DVGs for ZIKV and other flaviviruses.

Extensive editing of cellular and viral double-stranded RNA structures accounts for innate immunity suppression and the proviral activity of ADAR1p150.

The interferon (IFN)-mediated innate immune response is the first line of defense against viruses. However, an IFN-stimulated gene, the adenosine deaminase acting on RNA 1 (ADAR1), favors the replication of several viruses. ADAR1 binds double-stranded RNA and converts adenosine to inosine by deamination. This form of editing makes duplex RNA unstable, thereby preventing IFN induction. To better understand how ADAR1 works at the cellular level, we generated cell lines that express exclusively either the IFN-inducible, cytoplasmic isoform ADAR1p150, the constitutively expressed nuclear isoform ADAR1p110, or no isoform. By comparing the transcriptome of these cell lines, we identified more than 150 polymerase II transcripts that are extensively edited, and we attributed most editing events to ADAR1p150. Editing is focused on inverted transposable elements, located mainly within introns and untranslated regions, and predicted to form duplex RNA structures. Editing of these elements occurs also in primary human samples, and there is evidence for cross-species evolutionary conservation of editing patterns in primates and, to a lesser extent, in rodents. Whereas ADAR1p150 rarely edits tightly encapsidated standard measles virus (MeV) genomes, it efficiently edits genomes with inverted repeats accidentally generated by a mutant MeV. We also show that immune activation occurs in fully ADAR1-deficient (ADAR1KO) cells, restricting virus growth, and that complementation of these cells with ADAR1p150 rescues virus growth and suppresses innate immunity activation. Finally, by knocking out either protein kinase R (PKR) or mitochondrial antiviral signaling protein (MAVS)-another protein controlling the response to duplex RNA-in ADAR1KO cells, we show that PKR activation elicits a stronger antiviral response. Thus, ADAR1 prevents innate immunity activation by cellular transcripts that include extensive duplex RNA structures. The trade-off is that viruses take advantage of ADAR1 to elude innate immunity control.

Mutational and fitness landscapes of an RNA virus revealed through population sequencing.

RNA viruses exist as genetically diverse populations. It is thought that diversity and genetic structure of viral populations determine the rapid adaptation observed in RNA viruses and hence their pathogenesis. However, our understanding of the mechanisms underlying virus evolution has been limited by the inability to accurately describe the genetic structure of virus populations. Next-generation sequencing technologies generate data of sufficient depth to characterize virus populations, but are limited in their utility because most variants are present at very low frequencies and are thus indistinguishable from next-generation sequencing errors. Here we present an approach that reduces next-generation sequencing errors and allows the description of virus populations with unprecedented accuracy. Using this approach, we define the mutation rates of poliovirus and uncover the mutation landscape of the population. Furthermore, by monitoring changes in variant frequencies on serially passaged populations, we determined fitness values for thousands of mutations across the viral genome. Mapping of these fitness values onto three-dimensional structures of viral proteins offers a powerful approach for exploring structure-function relationships and potentially uncovering new functions. To our knowledge, our study provides the first single-nucleotide fitness landscape of an evolving RNA virus and establishes a general experimental platform for studying the genetic changes underlying the evolution of virus populations.

FACT is a sensor of DNA torsional stress in eukaryotic cells.

Transitions of B-DNA to alternative DNA structures (ADS) can be triggered by negative torsional strain, which occurs during replication and transcription, and may lead to genomic instability. However, how ADS are recognized in cells is unclear. We found that the binding of candidate anticancer drug, curaxin, to cellular DNA results in uncoiling of nucleosomal DNA, accumulation of negative supercoiling and conversion of multiple regions of genomic DNA into left-handed Z-form. Histone chaperone FACT binds rapidly to the same regions via the SSRP1 subunit in curaxin-treated cells. In vitro binding of purified SSRP1 or its isolated CID domain to a methylated DNA fragment containing alternating purine/pyrimidines, which is prone to Z-DNA transition, is much stronger than to other types of DNA. We propose that FACT can recognize and bind Z-DNA or DNA in transition from a B to Z form. Binding of FACT to these genomic regions triggers a p53 response. Furthermore, FACT has been shown to bind to other types of ADS through a different structural domain, which also leads to p53 activation. Thus, we propose that FACT acts as a sensor of ADS formation in cells. Recognition of ADS by FACT followed by a p53 response may explain the role of FACT in DNA damage prevention.

The evidence for a microRNA product of human DROSHA gene.

MicroRNAs are short RNA molecules that regulate function and stability of a large subset of eukaryotic mRNAs. In the main pathway of microRNA biogenesis, a short "hairpin" is excised from a primary transcript by ribonuclease DROSHA, followed by additional nucleolytic processing by DICER and inclusion of the mature microRNA into the RNA-induced silencing complex. We report that a microRNA-like molecule is encoded by human DROSHA gene within a predicted stem-loop element of the respective transcript. This putative mature microRNA is complementary to DROSHA transcript variant 1 and can attenuate expression of the corresponding protein. The findings suggest a possibility for a negative feedback loop, wherein DROSHA processes its own transcript and produces an inhibitor of its own biosynthesis.

High-resolution metabolic mapping of cell types in plant roots.

Metabolite composition offers a powerful tool for understanding gene function and regulatory processes. However, metabolomics studies on multicellular organisms have thus far been performed primarily on whole organisms, organs, or cell lines, losing information about individual cell types within a tissue. With the goal of profiling metabolite content in different cell populations within an organ, we used FACS to dissect GFP-marked cells from Arabidopsis roots for metabolomics analysis. Here, we present the metabolic profiles obtained from five GFP-tagged lines representing core cell types in the root. Fifty metabolites were putatively identified, with the most prominent groups being glucosinolates, phenylpropanoids, and dipeptides, the latter of which is not yet explored in roots. The mRNA expression of enzymes or regulators in the corresponding biosynthetic pathways was compared with the relative metabolite abundance. Positive correlations suggest that the rate-limiting steps in biosynthesis of glucosinolates in the root are oxidative modifications of side chains. The current study presents a work flow for metabolomics analyses of cell-type populations.

p53 cooperates with DNA methylation and a suicidal interferon response to maintain epigenetic silencing of repeats and noncoding RNAs.

Large parts of mammalian genomes are transcriptionally inactive and enriched with various classes of interspersed and tandem repeats. Here we show that the tumor suppressor protein p53 cooperates with DNA methylation to maintain silencing of a large portion of the mouse genome. Massive transcription of major classes of short, interspersed nuclear elements (SINEs) B1 and B2, both strands of near-centromeric satellite DNAs consisting of tandem repeats, and multiple species of noncoding RNAs was observed in p53-deficient but not in p53 wild-type mouse fibroblasts treated with the DNA demethylating agent 5-aza-2'-deoxycytidine. The abundance of these transcripts exceeded the level of ß-actin mRNA by more than 150-fold. Accumulation of these transcripts, which are capable of forming double-stranded RNA (dsRNA), was accompanied by a strong, endogenous, apoptosis-inducing type I IFN response. This phenomenon, which we named "TRAIN" (for "transcription of repeats activates interferon"), was observed in spontaneous tumors in two models of cancer-prone mice, presumably reflecting naturally occurring DNA hypomethylation and p53 inactivation in cancer. These observations suggest that p53 and IFN cooperate to prevent accumulation of cells with activated repeats and provide a plausible explanation for the deregulation of IFN function frequently seen in tumors. Overall, this work reveals roles for p53 and IFN that are key for genetic stability and therefore relevant to both tumorigenesis and the evolution of species.

Identification and characterization of a highly motile and antibiotic refractory subpopulation involved in the expansion of swarming colonies of Paenibacillus vortex.

Bacteria often use sophisticated cooperative behaviours, such as the development of complex colonies, elaborate biofilms and advanced dispersal strategies, to cope with the harsh and variable conditions of natural habitats, including the presence of antibiotics. Paenibacillus vortex uses swarming motility and cell-to-cell communication to form complex, structured colonies. The modular organization of P. vortex colony has been found to facilitate its dispersal on agar surfaces. The current study reveals that the complex structure of the colony is generated by the coexistence and transition between two morphotypes--'builders' and 'explorers'--with distinct functions in colony formation. Here, we focused on the explorers, which are highly motile and spearhead colonial expansion. Explorers are characterized by high expression levels of flagellar genes, such as flagellin (hag), motA, fliI, flgK and sigD, hyperflagellation, decrease in ATP (adenosine-5'-triphosphate) levels, and increased resistance to antibiotics. Their tolerance to many antibiotics gives them the advantage of translocation through antibiotics-containing areas. This work gives new insights on the importance of cell differentiation and task distribution in colony morphogenesis and adaptation to antibiotics.

A binary search approach to whole-genome data analysis.

A sequence analysis-oriented binary search-like algorithm was transformed to a sensitive and accurate analysis tool for processing whole-genome data. The advantage of the algorithm over previous methods is its ability to detect the margins of both short and long genome fragments, enriched by up-regulated signals, at equal accuracy. The score of an enriched genome fragment reflects the difference between the actual concentration of up-regulated signals in the fragment and the chromosome signal baseline. The "divide-and-conquer"-type algorithm detects a series of nonintersecting fragments of various lengths with locally optimal scores. The procedure is applied to detected fragments in a nested manner by recalculating the lower-than-baseline signals in the chromosome. The algorithm was applied to simulated whole-genome data, and its sensitivity/specificity were compared with those of several alternative algorithms. The algorithm was also tested with four biological tiling array datasets comprising Arabidopsis (i) expression and (ii) histone 3 lysine 27 trimethylation CHIP-on-chip datasets; Saccharomyces cerevisiae (iii) spliced intron data and (iv) chromatin remodeling factor binding sites. The analyses' results demonstrate the power of the algorithm in identifying both the short up-regulated fragments (such as exons and transcription factor binding sites) and the long--even moderately up-regulated zones--at their precise genome margins. The algorithm generates an accurate whole-genome landscape that could be used for cross-comparison of signals across the same genome in evolutionary and general genomic studies.

Physiological changes during torpor favor association with Endozoicomonas endosymbionts in the urochordate Botrylloides leachii.

Environmental perturbations evoke down-regulation of metabolism in some multicellular organisms, leading to dormancy, or torpor. Colonies of the urochordate Botrylloides leachii enter torpor in response to changes in seawater temperature and may survive for months as small vasculature remnants that lack feeding and reproductive organs but possess torpor-specific microbiota. Upon returning to milder conditions, the colonies rapidly restore their original morphology, cytology and functionality while harboring re-occurring microbiota, a phenomenon that has not been described in detail to date. Here we investigated the stability of B. leachii microbiome and its functionality in active and dormant colonies, using microscopy, qPCR, in situ hybridization, genomics and transcriptomics. A novel lineage of Endozoicomonas, proposed here as Candidatus Endozoicomonas endoleachii, was dominant in torpor animals (53-79% read abundance), and potentially occupied specific hemocytes found only in torpid animals. Functional analysis of the metagenome-assembled genome and genome-targeted transcriptomics revealed that Endozoicomonas can use various cellular substrates, like amino acids and sugars, potentially producing biotin and thiamine, but also expressing various features involved in autocatalytic symbiosis. Our study suggests that the microbiome can be linked to the metabolic and physiological states of the host, B. leachii, introducing a model organism for the study of symbioses during drastic physiological changes, such as torpor.

Genome sequence of the pattern-forming social bacterium Paenibacillus dendritiformis C454 chiral morphotype.

Paenibacillus dendritiformis is a Gram-positive, soil-dwelling, spore-forming social microorganism. An intriguing collective faculty of this strain is manifested by its ability to switch between different morphotypes, such as the branching (T) and the chiral (C) morphotypes. Here we report the 6.3-Mb draft genome sequence of the P. dendritiformis C454 chiral morphotype.

Tandem repeats on an eco-geographical scale: outcomes from the genome of Aegilops speltoides.

The chromosomal pattern of tandem repeat fractions of repetitive DNA is one of the most important characteristics of a species. In the present research, we aimed to detect and evaluate the level of intraspecific variability in the chromosomal distribution of species-specific Spelt 1 and Aegilops-Triticum-specific Spelt 52 tandem repeats in Aegilops speltoides and in closely related diploid and polyploid species. There is a distinct eco-geographical gradient in Spelt 1 and Spelt 52 blocks abundance in Ae. speltoides. In marginal populations, the number of Spelt 1 chromosomal blocks could be 12-14 times lower than in the center of the species distribution. Also, in related diploid species, the abundance of Spelt 52 correlates with evolutionary proximity to Ae. speltoides. Finally, the B- and G-genomes of allopolyploid wheats have Spelt 1 chromosomal distribution patterns similar to those of the types of Ae. speltoides with poor and rich contents of Spelt 1, respectively. The observed changes in numbers of blocks of Spelt 1 and Spelt 52 tandem repeats along the eco-geographical gradient may due to their depletion in the marginal populations as a result of increased recombination frequency under stressful conditions. Alternatively, it may be accumulation of tandem repeats in conducive climatic/edaphic environments in the center of the species' geographical distribution. Anyway, we observe a bidirectional shift of repetitive DNA genomic patterns on the population level leading to the formation of population-specific chromosomal patterns of tandem repeats. The appearance of a new chromosomal pattern is considered an important factor in promoting the emergence of interbreeding barriers.

Extensive metabolic cross-talk in melon fruit revealed by spatial and developmental combinatorial metabolomics.

• Variations in tissue development and spatial composition have a major impact on the nutritional and organoleptic qualities of ripe fleshy fruit, including melon (Cucumis melo). To gain a deeper insight into the mechanisms involved in these changes, we identified key metabolites for rational food quality design. • The metabolome, volatiles and mineral elements were profiled employing an unprecedented range of complementary analytical technologies. Fruits were followed at a number of time points during the final ripening process and tissues were collected across the fruit flesh from rind to seed cavity. Approximately 2000 metabolite signatures and 15 mineral elements were determined in an assessment of temporal and spatial melon fruit development. • This study design enabled the identification of: coregulated hubs (including aspartic acid, 2-isopropylmalic acid, ß-carotene, phytoene and dihydropseudoionone) in metabolic association networks; global patterns of coordinated compositional changes; and links of primary and secondary metabolism to key mineral and volatile fruit complements. • The results reveal the extent of metabolic interactions relevant to ripe fruit quality and thus have enabled the identification of essential candidate metabolites for the high-throughput screening of melon breeding populations for targeted breeding programmes aimed at nutrition and flavour improvement.

Development of infrastructure for a systemic multidisciplinary approach to study aging in retired sled dogs.

Canines represent a valuable model for mammalian aging studies as large animals with short lifespans, allowing longitudinal analyses within a reasonable time frame. Moreover, they develop a spectrum of aging-related diseases resembling that of humans, are exposed to similar environments, and have been reasonably well studied in terms of physiology and genetics. To overcome substantial variables that complicate studies of privately-owned household dogs, we have focused on a more uniform population composed of retired Alaskan sled dogs that shared similar lifestyles, including exposure to natural stresses, and are less prone to breed-specific biases than a pure breed population. To reduce variability even further, we have collected a population of 103 retired (8-11 years-old) sled dogs from multiple North American kennels in a specialized research facility named Vaika. Vaika dogs are maintained under standardized conditions with professional veterinary care and participate in a multidisciplinary program to assess the longitudinal dynamics of aging. The established Vaika infrastructure enables periodic gathering of quantitative data reflecting physical, physiological, immunological, neurological, and cognitive decline, as well as monitoring of aging-associated genetic and epigenetic alterations occurring in somatic cells. In addition, we assess the development of age-related diseases such as arthritis and cancer. In-depth data analysis, including artificial intelligence-based approaches, will build a comprehensive, integrated model of canine aging and potentially identify aging biomarkers that will allow use of this model for future testing of antiaging therapies.

Comparative Genomics Analysis of Repetitive Elements in Ten Gymnosperm Species: "Dark Repeatome" and Its Abundance in Conifer and Gnetum Species.

Repetitive elements (RE) and transposons (TE) can comprise up to 80% of some plant genomes and may be essential for regulating their evolution and adaptation. The "repeatome" information is often unavailable in assembled genomes because genomic areas of repeats are challenging to assemble and are often missing from final assembly. However, raw genomic sequencing data contain rich information about RE/TEs. Here, raw genomic NGS reads of 10 gymnosperm species were studied for the content and abundance patterns of their "repeatome". We utilized a combination of alignment on databases of repetitive elements and de novo assembly of highly repetitive sequences from genomic sequencing reads to characterize and calculate the abundance of known and putative repetitive elements in the genomes of 10 conifer plants: Pinus taeda, Pinus sylvestris, Pinus sibirica, Picea glauca, Picea abies, Abies sibirica, Larix sibirica, Juniperus communis, Taxus baccata, and Gnetum gnemon. We found that genome abundances of known and newly discovered putative repeats are specific to phylogenetically close groups of species and match biological taxa. The grouping of species based on abundances of known repeats closely matches the grouping based on abundances of newly discovered putative repeats (kChains) and matches the known taxonomic relations.