Genomic epidemiology of SARS-CoV-2 in the UAE reveals novel virus mutation, patterns of co-infection and tissue specific host immune response.

To unravel the source of SARS-CoV-2 introduction and the pattern of its spreading and evolution in the United Arab Emirates, we conducted meta-transcriptome sequencing of 1067 nasopharyngeal swab samples collected between May 9th and Jun 29th, 2020 during the first peak of the local COVID-19 epidemic. We identified global clade distribution and eleven novel genetic variants that were almost absent in the rest of the world and that defined five subclades specific to the UAE viral population. Cross-settlement human-to-human transmission was related to the local business activity. Perhaps surprisingly, at least 5% of the population were co-infected by SARS-CoV-2 of multiple clades within the same host. We also discovered an enrichment of cytosine-to-uracil mutation among the viral population collected from the nasopharynx, that is different from the adenosine-to-inosine change previously reported in the bronchoalveolar lavage fluid samples and a previously unidentified upregulation of APOBEC4 expression in nasopharynx among infected patients, indicating the innate immune host response mediated by ADAR and APOBEC gene families could be tissue-specific. The genomic epidemiological and molecular biological knowledge reported here provides new insights for the SARS-CoV-2 evolution and transmission and points out future direction on host-pathogen interaction investigation.

Author Correction: Elephant shark genome provides unique insights into gnathostome evolution.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

C/EBPα mediates the growth inhibitory effect of progestins on breast cancer cells.

Steroid hormones are key gene regulators in breast cancer cells. While estrogens stimulate cell proliferation, progestins activate a single cell cycle followed by proliferation arrest. Here, we use biochemical and genome-wide approaches to show that progestins achieve this effect via a functional crosstalk with C/EBPα. Using ChIP-seq, we identify around 1,000 sites where C/EBPα binding precedes and helps binding of progesterone receptor (PR) in response to hormone. These regions exhibit epigenetic marks of active enhancers, and C/EBPα maintains an open chromatin conformation that facilitates loading of ligand-activated PR. Prior to hormone exposure, C/EBPα favors promoter-enhancer contacts that assure hormonal regulation of key genes involved in cell proliferation by facilitating binding of RAD21, YY1, and the Mediator complex. Knockdown of C/EBPα disrupts enhancer-promoter contacts and decreases the presence of these architectural proteins, highlighting its key role in 3D chromatin looping. Thus, C/EBPα fulfills a previously unknown function as a potential growth modulator in hormone-dependent breast cancer.

Spanish list of potentially inappropriate drugs in the elderly (ES-PIA project).

PURPOSE: In the last decades, different criteria have been developed for detecting inappropriate prescription in older patients. In Spain, translations and adaptations of international lists are available but it would be necessary a national list which could cope with the peculiarities of our health system, existing pharmaceutical market, and prescription habits. We propose in this project the creation of a Spanish potentially inappropriate drugs list which could be applicable in our clinical scenario. METHODS: We use a Delphi method involving 25 experts from different backgrounds (Clinical Pharmacology, Geriatrics, Rational Use of Drugs and Pharmacy, Primary Care and Pharmacoepidemiology, and Pharmacovigilance) that were asked to participate in two-round questionnaires. For analysis, current recommendations of Worth and Pigni were applied, and every statement was classified into one of three groups: strong, moderate, or low agreement. Statements with strong agreement were accepted to be part of the inadequate prescription list. Moderate agreement statements were selected to enter the second questionnaire, and statements with low agreement were further analyzed to determine if it was due to heterogeneity or due to dispersion in the answers. RESULTS: The first questionnaire consisted of 160 proposed sentences, of which 106 reached a high agreement, 32 a moderate agreement, and 22 a low agreement. All sentences proposed in the second questionnaire reached a strong agreement. The total accepted sentences were 138. CONCLUSIONS: We offer a list of inadequate prescription in older patients adapted to the Spanish pharmacopeia and according to the prescription habits in our environment.

Rapid reversible changes in compartments and local chromatin organization revealed by hyperosmotic shock.

Nuclear architecture is decisive for the assembly of transcriptional responses. However, how chromosome organization is dynamically modulated to permit rapid and transient transcriptional changes in response to environmental challenges remains unclear. Here we show that hyperosmotic stress disrupts different levels of chromosome organization, ranging from A/B compartment changes to reduction in the number and insulation of topologically associating domains (TADs). Concomitantly, transcription is greatly affected, TAD borders weaken, and RNA Polymerase II runs off from hundreds of transcription end sites. Stress alters the binding profiles of architectural proteins, which explains the disappearance of local chromatin organization. These processes are dynamic, and cells rapidly reconstitute their default chromatin conformation after stress removal, uncovering an intrinsic organization. Transcription is not required for local chromatin reorganization, while compartment recovery is partially transcription-dependent. Thus, nuclear organization in mammalian cells can be rapidly modulated by environmental changes in a reversible manner.

Arginine Citrullination at the C-Terminal Domain Controls RNA Polymerase II Transcription.

The post-translational modification of key residues at the C-terminal domain of RNA polymerase II (RNAP2-CTD) coordinates transcription, splicing, and RNA processing by modulating its capacity to act as a landing platform for a variety of protein complexes. Here, we identify a new modification at the CTD, the deimination of arginine and its conversion to citrulline by peptidyl arginine deiminase 2 (PADI2), an enzyme that has been associated with several diseases, including cancer. We show that, among PADI family members, only PADI2 citrullinates R1810 (Cit1810) at repeat 31 of the CTD. Depletion of PADI2 or loss of R1810 results in accumulation of RNAP2 at transcription start sites, reduced gene expression, and inhibition of cell proliferation. Cit1810 is needed for interaction with the P-TEFb (positive transcription elongation factor b) kinase complex and for its recruitment to chromatin. In this way, CTD-Cit1810 favors RNAP2 pause release and efficient transcription in breast cancer cells.

Hormone-control regions mediate steroid receptor-dependent genome organization.

In breast cancer cells, some topologically associating domains (TADs) behave as hormonal gene regulation units, within which gene transcription is coordinately regulated in response to steroid hormones. Here we further describe that responsive TADs contain 20- to 100-kb-long clusters of intermingled estrogen receptor (ESR1) and progesterone receptor (PGR) binding sites, hereafter called hormone-control regions (HCRs). In T47D cells, we identified more than 200 HCRs, which are frequently bound by unliganded ESR1 and PGR. These HCRs establish steady long-distance inter-TAD interactions between them and organize characteristic looping structures with promoters in their TADs even in the absence of hormones in ESR1+-PGR+ cells. This organization is dependent on the expression of the receptors and is further dynamically modulated in response to steroid hormones. HCRs function as platforms that integrate different signals, resulting in some cases in opposite transcriptional responses to estrogens or progestins. Altogether, these results suggest that steroid hormone receptors act not only as hormone-regulated sequence-specific transcription factors but also as local and global genome organizers.

Unliganded Progesterone Receptor Governs Estrogen Receptor Gene Expression by Regulating DNA Methylation in Breast Cancer Cells.

Breast cancer prognosis and response to endocrine therapy strongly depends on the expression of the estrogen and progesterone receptors (ER and PR, respectively). Although much is known about ERα gene (ESR1) regulation after hormonal stimulation, how it is regulated in hormone-free condition is not fully understood. We used ER-/PR-positive breast cancer cells to investigate the role of PR in ESR1 regulation in the absence of hormones. We show that PR binds to the low-methylated ESR1 promoter and maintains both gene expression and DNA methylation of the ESR1 locus in hormone-deprived breast cancer cells. Depletion of PR reduces ESR1 expression, with a concomitant increase in gene promoter methylation. The high amount of methylation in the ESR1 promoter of PR-depleted cells persists after the stable re-expression of PR and inhibits PR binding to this genomic region. As a consequence, the rescue of PR expression in PR-depleted cells is insufficient to restore ESR1 expression. Consistently, DNA methylation impedes PR binding to consensus progesterone responsive elements. These findings contribute to understanding the complex crosstalk between PR and ER and suggest that the analysis of ESR1 promoter methylation in breast cancer cells can help to design more appropriate targeted therapies for breast cancer patients.

OneD: increasing reproducibility of Hi-C samples with abnormal karyotypes.

The three-dimensional conformation of genomes is an essential component of their biological activity. The advent of the Hi-C technology enabled an unprecedented progress in our understanding of genome structures. However, Hi-C is subject to systematic biases that can compromise downstream analyses. Several strategies have been proposed to remove those biases, but the issue of abnormal karyotypes received little attention. Many experiments are performed in cancer cell lines, which typically harbor large-scale copy number variations that create visible defects on the raw Hi-C maps. The consequences of these widespread artifacts on the normalized maps are mostly unexplored. We observed that current normalization methods are not robust to the presence of large-scale copy number variations, potentially obscuring biological differences and enhancing batch effects. To address this issue, we developed an alternative approach designed to take into account chromosomal abnormalities. The method, called OneD, increases reproducibility among replicates of Hi-C samples with abnormal karyotype, outperforming previous methods significantly. On normal karyotypes, OneD fared equally well as state-of-the-art methods, making it a safe choice for Hi-C normalization. OneD is fast and scales well in terms of computing resources for resolutions up to 5 kb.

Transcription factors orchestrate dynamic interplay between genome topology and gene regulation during cell reprogramming.

Chromosomal architecture is known to influence gene expression, yet its role in controlling cell fate remains poorly understood. Reprogramming of somatic cells into pluripotent stem cells (PSCs) by the transcription factors (TFs) OCT4, SOX2, KLF4 and MYC offers an opportunity to address this question but is severely limited by the low proportion of responding cells. We have recently developed a highly efficient reprogramming protocol that synchronously converts somatic into pluripotent stem cells. Here, we used this system to integrate time-resolved changes in genome topology with gene expression, TF binding and chromatin-state dynamics. The results showed that TFs drive topological genome reorganization at multiple architectural levels, often before changes in gene expression. Removal of locus-specific topological barriers can explain why pluripotency genes are activated sequentially, instead of simultaneously, during reprogramming. Together, our results implicate genome topology as an instructive force for implementing transcriptional programs and cell fate in mammals.

Similar genomic proportions of copy number variation within gray wolves and modern dog breeds inferred from whole genome sequencing.

BACKGROUND: Whole genome re-sequencing data from dogs and wolves are now commonly used to study how natural and artificial selection have shaped the patterns of genetic diversity. Single nucleotide polymorphisms, microsatellites and variants in mitochondrial DNA have been interrogated for links to specific phenotypes or signals of domestication. However, copy number variation (CNV), despite its increasingly recognized importance as a contributor to phenotypic diversity, has not been extensively explored in canids. RESULTS: Here, we develop a new accurate probabilistic framework to create fine-scale genomic maps of segmental duplications (SDs), compare patterns of CNV across groups and investigate their role in the evolution of the domestic dog by using information from 34 canine genomes. Our analyses show that duplicated regions are enriched in genes and hence likely possess functional importance. We identify 86 loci with large CNV differences between dogs and wolves, enriched in genes responsible for sensory perception, immune response, metabolic processes, etc. In striking contrast to the observed loss of nucleotide diversity in domestic dogs following the population bottlenecks that occurred during domestication and breed creation, we find a similar proportion of CNV loci in dogs and wolves, suggesting that other dynamics are acting to particularly select for CNVs with potentially functional impacts. CONCLUSIONS: This work is the first comparison of genome wide CNV patterns in domestic and wild canids using whole-genome sequencing data and our findings contribute to study the impact of novel kinds of genetic changes on the evolution of the domestic dog.

Parallel sequencing lives, or what makes large sequencing projects successful.

T47D_rep2 and b1913e6c1_51720e9cf were 2 Hi-C samples. They were born and processed at the same time, yet their fates were very different. The life of b1913e6c1_51720e9cf was simple and fruitful, while that of T47D_rep2 was full of accidents and sorrow. At the heart of these differences lies the fact that b1913e6c1_51720e9cf was born under a lab culture of Documentation, Automation, Traceability, and Autonomy and compliance with the FAIR Principles. Their lives are a lesson for those who wish to embark on the journey of managing high-throughput sequencing data.

Extreme genomic erosion after recurrent demographic bottlenecks in the highly endangered Iberian lynx.

BACKGROUND: Genomic studies of endangered species provide insights into their evolution and demographic history, reveal patterns of genomic erosion that might limit their viability, and offer tools for their effective conservation. The Iberian lynx (Lynx pardinus) is the most endangered felid and a unique example of a species on the brink of extinction. RESULTS: We generate the first annotated draft of the Iberian lynx genome and carry out genome-based analyses of lynx demography, evolution, and population genetics. We identify a series of severe population bottlenecks in the history of the Iberian lynx that predate its known demographic decline during the 20th century and have greatly impacted its genome evolution. We observe drastically reduced rates of weak-to-strong substitutions associated with GC-biased gene conversion and increased rates of fixation of transposable elements. We also find multiple signatures of genetic erosion in the two remnant Iberian lynx populations, including a high frequency of potentially deleterious variants and substitutions, as well as the lowest genome-wide genetic diversity reported so far in any species. CONCLUSIONS: The genomic features observed in the Iberian lynx genome may hamper short- and long-term viability through reduced fitness and adaptive potential. The knowledge and resources developed in this study will boost the research on felid evolution and conservation genomics and will benefit the ongoing conservation and management of this emblematic species.

Hormone-induced repression of genes requires BRG1-mediated H1.2 deposition at target promoters.

Eukaryotic gene regulation is associated with changes in chromatin compaction that modulate access to DNA regulatory sequences relevant for transcriptional activation or repression. Although much is known about the mechanism of chromatin remodeling in hormonal gene activation, how repression is accomplished is much less understood. Here we report that in breast cancer cells, ligand-activated progesterone receptor (PR) is directly recruited to transcriptionally repressed genes involved in cell proliferation along with the kinases ERK1/2 and MSK1. PR recruits BRG1 associated with the HP1γ-LSD1 complex repressor complex, which is further anchored via binding of HP1γ to the H3K9me3 signal deposited by SUV39H2. In contrast to what is observed during gene activation, only BRG1 and not the BAF complex is recruited to repressed promoters, likely due to local enrichment of the pioneer factor FOXA1. BRG1 participates in gene repression by interacting with H1.2, facilitating its deposition and stabilizing nucleosome positioning around the transcription start site. Our results uncover a mechanism of hormone-dependent transcriptional repression and a novel role for BRG1 in progestin regulation of breast cancer cell growth.

Polymorphic tandem repeats within gene promoters act as modifiers of gene expression and DNA methylation in humans.

Despite representing an important source of genetic variation, tandem repeats (TRs) remain poorly studied due to technical difficulties. We hypothesized that TRs can operate as expression (eQTLs) and methylation (mQTLs) quantitative trait loci. To test this we analyzed the effect of variation at 4849 promoter-associated TRs, genotyped in 120 individuals, on neighboring gene expression and DNA methylation. Polymorphic promoter TRs were associated with increased variance in local gene expression and DNA methylation, suggesting functional consequences related to TR variation. We identified >100 TRs associated with expression/methylation levels of adjacent genes. These potential eQTL/mQTL TRs were enriched for overlaps with transcription factor binding and DNaseI hypersensitivity sites, providing a rationale for their effects. Moreover, we showed that most TR variants are poorly tagged by nearby single nucleotide polymorphisms (SNPs) markers, indicating that many functional TR variants are not effectively assayed by SNP-based approaches. Our study assigns biological significance to TR variations in the human genome, and suggests that a significant fraction of TR variations exert functional effects via alterations of local gene expression or epigenetics. We conclude that targeted studies that focus on genotyping TR variants are required to fully ascertain functional variation in the genome.

Germline methylation patterns determine the distribution of recombination events in the dog genome.

The positive-regulatory domain containing nine gene, PRDM9, which strongly associates with the location of recombination events in several vertebrates, is inferred to be inactive in the dog genome. Here, we address several questions regarding the control of recombination and its influence on genome evolution in dogs. First, we address whether the association between CpG islands (CGIs) and recombination hotspots is generated by lack of methylation, GC-biased gene conversion (gBGC), or both. Using a genome-wide dog single nucleotide polymorphism data set and comparisons of the dog genome with related species, we show that recombination-associated CGIs have low CpG mutation rates, and that CpG mutation rate is negatively correlated with recombination rate genome wide, indicating that nonmethylation attracts the recombination machinery. We next use a neighbor-dependent model of nucleotide substitution to disentangle the effects of CpG mutability and gBGC and analyze the effects that loss of PRDM9 has on these rates. We infer that methylation patterns have been stable during canid genome evolution, but that dog CGIs have experienced a drastic increase in substitution rate due to gBGC, consistent with increased levels of recombination in these regions. We also show that gBGC is likely to have generated many new CGIs in the dog genome, but these mostly occur away from genes, whereas the number of CGIs in gene promoter regions has not increased greatly in recent evolutionary history. Recombination has a major impact on the distribution of CGIs that are detected in the dog genome due to the interaction between methylation and gBGC. The results indicate that germline methylation patterns are the main determinant of recombination rates in the absence of PRDM9.

Comparative analysis of the domestic cat genome reveals genetic signatures underlying feline biology and domestication.

Little is known about the genetic changes that distinguish domestic cat populations from their wild progenitors. Here we describe a high-quality domestic cat reference genome assembly and comparative inferences made with other cat breeds, wildcats, and other mammals. Based upon these comparisons, we identified positively selected genes enriched for genes involved in lipid metabolism that underpin adaptations to a hypercarnivorous diet. We also found positive selection signals within genes underlying sensory processes, especially those affecting vision and hearing in the carnivore lineage. We observed an evolutionary tradeoff between functional olfactory and vomeronasal receptor gene repertoires in the cat and dog genomes, with an expansion of the feline chemosensory system for detecting pheromones at the expense of odorant detection. Genomic regions harboring signatures of natural selection that distinguish domestic cats from their wild congeners are enriched in neural crest-related genes associated with behavior and reward in mouse models, as predicted by the domestication syndrome hypothesis. Our description of a previously unidentified allele for the gloving pigmentation pattern found in the Birman breed supports the hypothesis that cat breeds experienced strong selection on specific mutations drawn from random bred populations. Collectively, these findings provide insight into how the process of domestication altered the ancestral wildcat genome and build a resource for future disease mapping and phylogenomic studies across all members of the Felidae.

Annotated features of domestic cat - Felis catus genome.

BACKGROUND: Domestic cats enjoy an extensive veterinary medical surveillance which has described nearly 250 genetic diseases analogous to human disorders. Feline infectious agents offer powerful natural models of deadly human diseases, which include feline immunodeficiency virus, feline sarcoma virus and feline leukemia virus. A rich veterinary literature of feline disease pathogenesis and the demonstration of a highly conserved ancestral mammal genome organization make the cat genome annotation a highly informative resource that facilitates multifaceted research endeavors. FINDINGS: Here we report a preliminary annotation of the whole genome sequence of Cinnamon, a domestic cat living in Columbia (MO, USA), bisulfite sequencing of Boris, a male cat from St. Petersburg (Russia), and light 30× sequencing of Sylvester, a European wildcat progenitor of cat domestication. The annotation includes 21,865 protein-coding genes identified by a comparative approach, 217 loci of endogenous retrovirus-like elements, repetitive elements which comprise about 55.7% of the whole genome, 99,494 new SNVs, 8,355 new indels, 743,326 evolutionary constrained elements, and 3,182 microRNA homologues. The methylation sites study shows that 10.5% of cat genome cytosines are methylated. An assisted assembly of a European wildcat, Felis silvestris silvestris, was performed; variants between F. silvestris and F. catus genomes were derived and compared to F. catus. CONCLUSIONS: The presented genome annotation extends beyond earlier ones by closing gaps of sequence that were unavoidable with previous low-coverage shotgun genome sequencing. The assembly and its annotation offer an important resource for connecting the rich veterinary and natural history of cats to genome discovery.

Digital genotyping of macrosatellites and multicopy genes reveals novel biological functions associated with copy number variation of large tandem repeats.

Tandem repeats are common in eukaryotic genomes, but due to difficulties in assaying them remain poorly studied. Here, we demonstrate the utility of Nanostring technology as a targeted approach to perform accurate measurement of tandem repeats even at extremely high copy number, and apply this technology to genotype 165 HapMap samples from three different populations and five species of non-human primates. We observed extreme variability in copy number of tandemly repeated genes, with many loci showing 5-10 fold variation in copy number among humans. Many of these loci show hallmarks of genome assembly errors, and the true copy number of many large tandem repeats is significantly under-represented even in the high quality 'finished' human reference assembly. Importantly, we demonstrate that most large tandem repeat variations are not tagged by nearby SNPs, and are therefore essentially invisible to SNP-based GWAS approaches. Using association analysis we identify many cis correlations of large tandem repeat variants with nearby gene expression and DNA methylation levels, indicating that variations of tandem repeat length are associated with functional effects on the local genomic environment. This includes an example where expansion of a macrosatellite repeat is associated with increased DNA methylation and suppression of nearby gene expression, suggesting a mechanism termed "repeat induced gene silencing", which has previously been observed only in transgenic organisms. We also observed multiple signatures consistent with altered selective pressures at tandemly repeated loci, suggesting important biological functions. Our studies show that tandemly repeated loci represent a highly variable fraction of the genome that have been systematically ignored by most previous studies, copy number variation of which can exert functionally significant effects. We suggest that future studies of tandem repeat loci will lead to many novel insights into their role in modulating both genomic and phenotypic diversity.