How much should we sequence? An analysis of the Swiss SARS-CoV-2 surveillance effort.

During the SARS-CoV-2 pandemic, many countries directed substantial resources toward genomic surveillance to detect and track viral variants. There is a debate over how much sequencing effort is necessary in national surveillance programs for SARS-CoV-2 and future pandemic threats. We aimed to investigate the effect of reduced sequencing on surveillance outcomes in a large genomic data set from Switzerland, comprising more than 143k sequences. We employed a uniform downsampling strategy using 100 iterations each to investigate the effects of fewer available sequences on the surveillance outcomes: (i) first detection of variants of concern (VOCs), (ii) speed of introduction of VOCs, (iii) diversity of lineages, (iv) first cluster detection of VOCs, (v) density of active clusters, and (vi) geographic spread of clusters. The impact of downsampling on VOC detection is disparate for the three VOC lineages, but many outcomes including introduction and cluster detection could be recapitulated even with only 35% of the original sequencing effort. The effect on the observed speed of introduction and first detection of clusters was more sensitive to reduced sequencing effort for some VOCs, in particular Omicron and Delta, respectively. A genomic surveillance program needs a balance between societal benefits and costs. While the overall national dynamics of the pandemic could be recapitulated by a reduced sequencing effort, the effect is strongly lineage-dependent-something that is unknown at the time of sequencing-and comes at the cost of accuracy, in particular for tracking the emergence of potential VOCs.IMPORTANCESwitzerland had one of the most comprehensive genomic surveillance systems during the COVID-19 pandemic. Such programs need to strike a balance between societal benefits and program costs. Our study aims to answer the question: How would surveillance outcomes have changed had we sequenced less? We find that some outcomes but also certain viral lineages are more affected than others by sequencing less. However, sequencing to around a third of the original effort still captured many important outcomes for the variants of concern such as their first detection but affected more strongly other measures like the detection of first transmission clusters for some lineages. Our work highlights the importance of setting predefined targets for a national genomic surveillance program based on which sequencing effort should be determined. Additionally, the use of a centralized surveillance platform facilitates aggregating data on a national level for rapid public health responses as well as post-analyses.

Pharmacological Induction of a Progenitor State for the Efficient Expansion of Primary Human Hepatocytes.

The liver is an organ with strong regenerative capacity, yet primary hepatocytes have a low amplification potential in vitro, a major limitation for the cell-based therapy of liver disorders and for ex vivo biological screens. Induced pluripotent stem cells (iPSCs) may help to circumvent this obstacle but often harbor genetic and epigenetic abnormalities, limiting their potential. Here, we describe the pharmacological induction of proliferative human hepatic progenitor cells (HPCs) through a cocktail of growth factors and small molecules mimicking the signaling events involved in liver regeneration. Human HPCs from healthy donors and pediatric patients proliferated vigorously while maintaining their genomic stability and could be redifferentiated in vitro into metabolically competent cells that supported the replication of hepatitis B and delta viruses. Redifferentiation efficiency was boosted by three-dimensional culture. Finally, transcriptome analysis showed that HPCs were more closely related to mature hepatocytes than iPSC-derived hepatocyte-like cells were. Conclusion: HPC induction holds promise for a variety of applications such as ex vivo disease modeling, personalized drug testing or metabolic studies, and development of a bioartificial liver.

Molecular Criteria for Defining the Naive Human Pluripotent State.

Recent studies have aimed to convert cultured human pluripotent cells to a naive state, but it remains unclear to what extent the resulting cells recapitulate in vivo naive pluripotency. Here we propose a set of molecular criteria for evaluating the naive human pluripotent state by comparing it to the human embryo. We show that transcription of transposable elements provides a sensitive measure of the concordance between pluripotent stem cells and early human development. We also show that induction of the naive state is accompanied by genome-wide DNA hypomethylation, which is reversible except at imprinted genes, and that the X chromosome status resembles that of the human preimplantation embryo. However, we did not see efficient incorporation of naive human cells into mouse embryos. Overall, the different naive conditions we tested showed varied relationships to human embryonic states based on molecular criteria, providing a backdrop for future analysis of naive human pluripotency.

Human Hepatocyte-Derived Induced Pluripotent Stem Cells: MYC Expression, Similarities to Human Germ Cell Tumors, and Safety Issues.

Induced pluripotent stem cells (iPSC) are a most promising approach to the development of a hepatocyte transplantable mass sufficient to induce long-term correction of inherited liver metabolic diseases, thus avoiding liver transplantation. Their intrinsic self-renewal ability and potential to differentiate into any of the three germ layers identify iPSC as the most promising cell-based therapeutics, but also as drivers of tumor development. Teratoma development currently represents the gold standard to assess iPSC pluripotency. We analyzed the tumorigenic potential of iPSC generated from human hepatocytes (HEP-iPSC) and compared their immunohistochemical profiles to that of tumors developed from fibroblast and hematopoietic stem cell-derived iPSC. HEP-iPSC generated tumors significantly presented more malignant morphological features than reprogrammed fibroblasts or CD34+ iPSC. Moreover, the protooncogene myc showed the strongest expression in HEP-iPSC, compared to only faint expression in the other cell subsets. Random integration of transgenes and the use of potent protooncogenes such as myc might be a risk factor for malignant tumor development if hepatocytes are used for reprogramming. Nonviral vector delivery systems or reprogramming of cells obtained from less invasive harvesting methods would represent interesting options for future developments in stem cell-based approaches for liver metabolic diseases.

Drawing a fine line on endogenous retroelement activity.

Endogenous retroelements (EREs) are essential motors of evolution yet require careful control to prevent genomic catastrophes, notably during the vulnerable phases of epigenetic reprogramming that occur immediately after fertilization and in germ cells. Accordingly, a variety of mechanisms restrict these mobile genetic units. Previous studies have revealed the importance of KRAB-containing zinc finger proteins (KRAB-ZFPs) and their cofactor, KAP1, in the early embryonic silencing of endogenous retroviruses and so-called SVAs, but the implication of this transcriptional repression system in the control of LINE-1, the only known active autonomous retrotransposon in the human genome, was thought to be marginal. Two recent studies straighten the record by revealing that the KRAB/KAP system is key to the control of L1 in embryonic stem (ES) cells, and go further in demonstrating that DNA methylation and KRAB/KAP1-induced repression contribute to this process in an evolutionally dynamic fashion. These results shed light on the delicate equilibrium between higher vertebrates and endogenous retroelements, which are not just genetic invaders calling for strict control but rather a constantly renewed and nicely exploitable source of evolutionary potential.

The developmental control of transposable elements and the evolution of higher species.

Transposable elements (TEs) account for at least 50% of the human genome. They constitute essential motors of evolution through their ability to modify genomic architecture, mutate genes and regulate gene expression. Accordingly, TEs are subject to tight epigenetic control during the earliest phases of embryonic development via histone and DNA methylation. Key to this process is recognition by sequence-specific RNA- and protein-based repressors. Collectively, these mediators are responsible for silencing a very broad range of TEs in an evolutionarily dynamic fashion. As a consequence, mobile elements and their controllers exert a marked influence on transcriptional networks in embryonic stem cells and a variety of adult tissues. The emerging picture is not that of a simple arms race but rather of a massive and sophisticated enterprise of TE domestication for the evolutionary benefit of the host.

Evolutionally dynamic L1 regulation in embryonic stem cells.

Mobile elements are important evolutionary forces that challenge genomic integrity. Long interspersed element-1 (L1, also known as LINE-1) is the only autonomous transposon still active in the human genome. It displays an unusual pattern of evolution, with, at any given time, a single active L1 lineage amplifying to thousands of copies before getting replaced by a new lineage, likely under pressure of host restriction factors, which act notably by silencing L1 expression during early embryogenesis. Here, we demonstrate that in human embryonic stem (hES) cells, KAP1 (KRAB [Krüppel-associated box domain]-associated protein 1), the master cofactor of KRAB-containing zinc finger proteins (KRAB-ZFPs) previously implicated in the restriction of endogenous retroviruses, represses a discrete subset of L1 lineages predicted to have entered the ancestral genome between 26.8 million and 7.6 million years ago. In mice, we documented a similar chronologically conditioned pattern, albeit with a much contracted time scale. We could further identify an L1-binding KRAB-ZFP, suggesting that this rapidly evolving protein family is more globally responsible for L1 recognition. KAP1 knockdown in hES cells induced the expression of KAP1-bound L1 elements, but their younger, human-specific counterparts (L1Hs) were unaffected. Instead, they were stimulated by depleting DNA methyltransferases, consistent with recent evidence demonstrating that the PIWI-piRNA (PIWI-interacting RNA) pathway regulates L1Hs in hES cells. Altogether, these data indicate that the early embryonic control of L1 is an evolutionarily dynamic process and support a model in which newly emerged lineages are first suppressed by DNA methylation-inducing small RNA-based mechanisms before KAP1-recruiting protein repressors are selected.

Loss of transcriptional control over endogenous retroelements during reprogramming to pluripotency.

Endogenous retroelements (EREs) account for about half of the mouse or human genome, and their potential as insertional mutagens and transcriptional perturbators is suppressed by early embryonic epigenetic silencing. Here, we asked how ERE control is maintained during the generation of induced pluripotent stem cells (iPSCs), as this procedure involves profound epigenetic remodeling. We found that all EREs tested were markedly up-regulated during the reprogramming of either mouse embryonic fibroblasts, human CD34(+) cells, or human primary hepatocytes. At the iPSC stage, EREs of some classes were repressed, whereas others remained highly expressed, yielding a pattern somewhat reminiscent of that recorded in embryonic stem cells. However, variability persisted between individual iPSC clones in the control of specific ERE integrants. Both during reprogramming and in iPS cells, the up-regulation of specific EREs significantly impacted on the transcription of nearby cellular genes. While transcription triggered by specific ERE integrants at highly precise developmental stages may be an essential step toward obtaining pluripotent cells, the broad and unspecific unleashing of the repetitive genome observed here may contribute to the inefficiency of the reprogramming process and to the phenotypic heterogeneity of iPSCs.

De novo DNA methylation of endogenous retroviruses is shaped by KRAB-ZFPs/KAP1 and ESET.

Endogenous retroviruses (ERVs) undergo de novo DNA methylation during the first few days of mammalian embryogenesis, although the factors that control the targeting of this process are largely unknown. We asked whether KAP1 (KRAB-associated protein 1) is involved in this mechanism because of its previously defined role in maintaining the silencing of ERVs through the histone methyltransferase ESET and histone H3 lysine 9 trimethylation. Here, we demonstrate that introduced ERV sequences are sufficient to direct rapid de novo methylation of a flanked promoter in embryonic stem (ES) cells. This mechanism requires the presence of an ERV sequence-recognizing KRAB zinc-finger protein (ZFP) and both KAP1 and ESET. Furthermore, this process can also take place on a strong cellular promoter and leads to methylation signatures that are subsequently maintained in vivo throughout embryogenesis. Finally, we show that methylation of ERVs residing in the genome is affected by knockout of KAP1 in early embryos. KRAB-ZFPs, KAP1 and ESET are thus likely to be responsible for the early embryonic instatement of stable epigenetic marks at ERV-containing loci.

Chromosome conformation capture uncovers potential genome-wide interactions between human conserved non-coding sequences.

Comparative analyses of various mammalian genomes have identified numerous conserved non-coding (CNC) DNA elements that display striking conservation among species, suggesting that they have maintained specific functions throughout evolution. CNC function remains poorly understood, although recent studies have identified a role in gene regulation. We hypothesized that the identification of genomic loci that interact physically with CNCs would provide information on their functions. We have used circular chromosome conformation capture (4C) to characterize interactions of 10 CNCs from human chromosome 21 in K562 cells. The data provide evidence that CNCs are capable of interacting with loci that are enriched for CNCs. The number of trans interactions varies among CNCs; some show interactions with many loci, while others interact with few. Some of the tested CNCs are capable of driving the expression of a reporter gene in the mouse embryo, and associate with the oligodendrocyte genes OLIG1 and OLIG2. Our results underscore the power of chromosome conformation capture for the identification of targets of functional DNA elements and raise the possibility that CNCs exert their functions by physical association with defined genomic regions enriched in CNCs. These CNC-CNC interactions may in part explain their stringent conservation as a group of regulatory sequences.

A systematic enhancer screen using lentivector transgenesis identifies conserved and non-conserved functional elements at the Olig1 and Olig2 locus.

Finding sequences that control expression of genes is central to understanding genome function. Previous studies have used evolutionary conservation as an indicator of regulatory potential. Here, we present a method for the unbiased in vivo screen of putative enhancers in large DNA regions, using the mouse as a model. We cloned a library of 142 overlapping fragments from a 200 kb-long murine BAC in a lentiviral vector expressing LacZ from a minimal promoter, and used the resulting vectors to infect fertilized murine oocytes. LacZ staining of E11 embryos obtained by first using the vectors in pools and then testing individual candidates led to the identification of 3 enhancers, only one of which shows significant evolutionary conservation. In situ hybridization and 3C/4C experiments suggest that this enhancer, which is active in the neural tube and posterior diencephalon, influences the expression of the Olig1 and/or Olig2 genes. This work provides a new approach for the large-scale in vivo screening of transcriptional regulatory sequences, and further demonstrates that evolutionary conservation alone seems too limiting a criterion for the identification of enhancers.

Disease-causing 7.4 kb cis-regulatory deletion disrupting conserved non-coding sequences and their interaction with the FOXL2 promotor: implications for mutation screening.

To date, the contribution of disrupted potentially cis-regulatory conserved non-coding sequences (CNCs) to human disease is most likely underestimated, as no systematic screens for putative deleterious variations in CNCs have been conducted. As a model for monogenic disease we studied the involvement of genetic changes of CNCs in the cis-regulatory domain of FOXL2 in blepharophimosis syndrome (BPES). Fifty-seven molecularly unsolved BPES patients underwent high-resolution copy number screening and targeted sequencing of CNCs. Apart from three larger distant deletions, a de novo deletion as small as 7.4 kb was found at 283 kb 5' to FOXL2. The deletion appeared to be triggered by an H-DNA-induced double-stranded break (DSB). In addition, it disrupts a novel long non-coding RNA (ncRNA) PISRT1 and 8 CNCs. The regulatory potential of the deleted CNCs was substantiated by in vitro luciferase assays. Interestingly, Chromosome Conformation Capture (3C) of a 625 kb region surrounding FOXL2 in expressing cellular systems revealed physical interactions of three upstream fragments and the FOXL2 core promoter. Importantly, one of these contains the 7.4 kb deleted fragment. Overall, this study revealed the smallest distant deletion causing monogenic disease and impacts upon the concept of mutation screening in human disease and developmental disorders in particular.

Human nocturnal frontal lobe epilepsy: pharmocogenomic profiles of pathogenic nicotinic acetylcholine receptor beta-subunit mutations outside the ion channel pore.

Certain mutations in specific parts of the neuronal nicotinic acetylcholine receptor (nAChR) subunit genes CHRNA4, CHRNB2, and probably CHRNA2, can cause autosomal dominant nocturnal frontal lobe epilepsy (ADNFLE). All but one of the known causative mutations are located in the second transmembrane region (TM2), which serves as the major ion poreforming domain of the receptor. Functional characterization of these ADNFLE mutations has shown that although each mutant exhibits specific properties, they all confer a gain of function with increased sensitivity to acetylcholine. In this work, we characterize the second and third ADNFLE-associated mutations that are external to TM2 but affect different amino acid residues within the third transmembrane region (TM3). The two new CHRNB2 mutations were identified in three families of Turkish Cypriot, Scottish, and English origin. These TM3 mutations elicit the same gain of function pathomechanism as observed for the TM2 mutations with enhanced acetylcholine sensitivity, despite their unusual localization within the gene. Electrophysiological experiments, including single channel measurements, revealed that incorporation of these new mutant subunits does not affect the conductance of the ionic pore but increases the probability of opening. Determination of the sensitivity to nicotine for nAChRs carrying mutations in TM2 and TM3 showed clear differences in the direction and the extent to which the window current for nicotine sensitivity was shifted for individual mutations, indicating differences in pharmacogenomic properties that are not readily correlated with increased ACh affinity.

Characterization of mouse Dactylaplasia mutations: a model for human ectrodactyly SHFM3.

SHFM3 is a limb malformation characterized by the absence of central digits. It has been shown that this condition is associated with tandem duplications of about 500 kb at 10q24. The Dactylaplasia mice display equivalent limb defects and the two corresponding alleles (Dac1j and Dac2j) map in the region syntenic with the duplications in SHFM3. Dac1j was shown to be associated with an insertion of an unspecified ETn-like mouse endogenous transposon upstream of the Fbxw4 gene. Dac2j was also thought to be an insertion or a small inversion in intron 5 of Fbxw4, but the breakpoints and the exact molecular lesion have not yet been characterized. Here we report precise mapping and characterization of these alleles. We failed to identify any copy number differences within the SHFM3 orthologous genomic locus between Dac mutant and wild-type littermates, showing that the Dactylaplasia alleles are not associated with duplications of the region, in contrast with the described human SHFM3 cases. We further show that both Dac1j and Dac2j are caused by insertions of MusD retroelements that share 98% sequence identity. The differences between the nature of the human and mouse genomic abnormalities argue against models proposed so far that either envisioned SHFM3 as a local trisomy or Dac as a mutant allele of Fbxw4. Instead, both genetic conditions might lead to complex alterations of gene regulation mechanisms that would impair limb morphogenesis. Interestingly, the Dac2j mutation occurs within a highly conserved element that may represent a regulatory sequence for a neighboring gene.

Identification of a novel member of the CLIC family, CLIC6, mapping to 21q22.12.

Five members of the newly identified chloride intracellular channel (CLIC) gene family of intracellular chloride channels (CLIC1-CLIC5) have previously been described in humans. Here we report the molecular cloning and initial characterisation of two splice forms of a novel member of this family, CLIC6, mapping to human chromosome 21. Two essential features distinguish CLIC6 from other members of the family. The CLIC6 protein is significantly longer and the CLIC6 gene contains a GC rich segment, which encodes a 10 amino acid motif repeated 14 times in the amino-terminus. Surprisingly, the repeat is conserved in the lagomorphs, but not in the rodents lineage. The putative bovine orthologue of CLIC5, p64, also exhibits a repeated motif, which is different from that of CLIC6. Attempts to functionally characterise CLIC6 by voltage clamp failed to show any chloride channel activity. Hence, the exact function of this protein remains unknown.

Broad antiretroviral defence by human APOBEC3G through lethal editing of nascent reverse transcripts.

Viral replication usually requires that innate intracellular lines of defence be overcome, a task usually accomplished by specialized viral gene products. The virion infectivity factor (Vif) protein of human immunodeficiency virus (HIV) is required during the late stages of viral production to counter the antiviral activity of APOBEC3G (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3G; also known as CEM15), a protein expressed notably in human T lymphocytes. When produced in the presence of APOBEC3G, vif-defective virus is non-infectious. APOBEC3G is closely related to APOBEC1, the central component of an RNA-editing complex that deaminates a cytosine residue in apoB messenger RNA. APOBEC family members also have potent DNA mutator activity through dC deamination; however, whether the editing potential of APOBEC3G has any relevance to HIV inhibition is unknown. Here, we demonstrate that it does, as APOBEC3G exerts its antiviral effect during reverse transcription to trigger G-to-A hypermutation in the nascent retroviral DNA. We also find that APOBEC3G can act on a broad range of retroviruses in addition to HIV, suggesting that hypermutation by editing is a general innate defence mechanism against this important group of pathogens.