The genomic origins of the world's first farmers.

The precise genetic origins of the first Neolithic farming populations in Europe and Southwest Asia, as well as the processes and the timing of their differentiation, remain largely unknown. Demogenomic modeling of high-quality ancient genomes reveals that the early farmers of Anatolia and Europe emerged from a multiphase mixing of a Southwest Asian population with a strongly bottlenecked western hunter-gatherer population after the last glacial maximum. Moreover, the ancestors of the first farmers of Europe and Anatolia went through a period of extreme genetic drift during their westward range expansion, contributing highly to their genetic distinctiveness. This modeling elucidates the demographic processes at the root of the Neolithic transition and leads to a spatial interpretation of the population history of Southwest Asia and Europe during the late Pleistocene and early Holocene.

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.

In embryonic stem cells, ZFP57/KAP1 recognize a methylated hexanucleotide to affect chromatin and DNA methylation of imprinting control regions.

The maintenance of H3K9 and DNA methylation at imprinting control regions (ICRs) during early embryogenesis is key to the regulation of imprinted genes. Here, we reveal that ZFP57, its cofactor KAP1, and associated effectors bind selectively to the H3K9me3-bearing, DNA-methylated allele of ICRs in ES cells. KAP1 deletion induces a loss of heterochromatin marks at ICRs, whereas deleting ZFP57 or DNMTs leads to ICR DNA demethylation. Accordingly, we find that ZFP57 and KAP1 associated with DNMTs and hemimethylated DNA-binding NP95. Finally, we identify the methylated TGCCGC hexanucleotide as the motif that is recognized by ZFP57 in all ICRs and in several tens of additional loci, several of which are at least ZFP57-dependently methylated in ES cells. These results significantly advance our understanding of imprinting and suggest a general mechanism for the protection of specific loci against the wave of DNA demethylation that affects the mammalian genome during early embryogenesis.

Interplay of TRIM28 and DNA methylation in controlling human endogenous retroelements.

Reverse transcription-derived sequences account for at least half of the human genome. Although these retroelements are formidable motors of evolution, they can occasionally cause disease, and accordingly are inactivated during early embryogenesis through epigenetic mechanisms. In the mouse, at least for endogenous retroviruses, important mediators of this process are the tetrapod-specific KRAB-containing zinc finger proteins (KRAB-ZFPs) and their cofactor TRIM28. The present study demonstrates that KRAB/TRIM28-mediated regulation is responsible for controlling a very broad range of human-specific endogenous retroelements (EREs) in human embryonic stem (ES) cells and that it exerts, as a consequence, a marked effect on the transcriptional dynamics of these cells. It further reveals reciprocal dependence between TRIM28 recruitment at specific families of EREs and DNA methylation. It finally points to the importance of persistent TRIM28-mediated control of ERE transcriptional impact beyond their presumed inactivation by DNA methylation.

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.

TRIM28 repression of retrotransposon-based enhancers is necessary to preserve transcriptional dynamics in embryonic stem cells.

TRIM28 is critical for the silencing of endogenous retroviruses (ERVs) in embryonic stem (ES) cells. Here, we reveal that an essential impact of this process is the protection of cellular gene expression in early embryos from perturbation by cis-acting activators contained within these retroelements. In TRIM28-depleted ES cells, repressive chromatin marks at ERVs are replaced by histone modifications typical of active enhancers, stimulating transcription of nearby cellular genes, notably those harboring bivalent promoters. Correspondingly, ERV-derived sequences can repress or enhance expression from an adjacent promoter in transgenic embryos depending on their TRIM28 sensitivity in ES cells. TRIM28-mediated control of ERVs is therefore crucial not just to prevent retrotransposition, but more broadly to safeguard the transcriptional dynamics of early embryos.

Probable zoonotic leprosy in the southern United States.

BACKGROUND: In the southern region of the United States, such as in Louisiana and Texas, there are autochthonous cases of leprosy among native-born Americans with no history of foreign exposure. In the same region, as well as in Mexico, wild armadillos are infected with Mycobacterium leprae. METHODS: Whole-genome resequencing of M. leprae from one wild armadillo and three U.S. patients with leprosy revealed that the infective strains were essentially identical. Comparative genomic analysis of these strains and M. leprae strains from Asia and Brazil identified 51 single-nucleotide polymorphisms and an 11-bp insertion-deletion. We genotyped these polymorphic sites, in combination with 10 variable-number tandem repeats, in M. leprae strains obtained from 33 wild armadillos from five southern states, 50 U.S. outpatients seen at a clinic in Louisiana, and 64 Venezuelan patients, as well as in four foreign reference strains. RESULTS: The M. leprae genotype of patients with foreign exposure generally reflected their country of origin or travel history. However, a unique M. leprae genotype (3I-2-v1) was found in 28 of the 33 wild armadillos and 25 of the 39 U.S. patients who resided in areas where exposure to armadillo-borne M. leprae was possible. This genotype has not been reported elsewhere in the world. CONCLUSIONS: Wild armadillos and many patients with leprosy in the southern United States are infected with the same strain of M. leprae. Armadillos are a large natural reservoir for M. leprae, and leprosy may be a zoonosis in the region. (Funded by the National Institute of Allergy and Infectious Diseases and others.).

KAP1 regulates gene networks controlling mouse B-lymphoid cell differentiation and function.

Chromatin remodeling is fundamental for B-cell differentiation. In the present study, we explored the role of KAP1, the cofactor of KRAB-ZFP transcriptional repressors, in this process. B-lymphoid-specific Kap1-KO mice displayed reduced numbers of mature B cells, lower steady-state levels of Abs, and accelerated rates of decay of neutralizing Abs after viral immunization. Transcriptome analyses of Kap1-deleted B splenocytes revealed an up-regulation of PTEN, the enzymatic counteractor of PIK3 signaling, and of genes encoding DNA-damage response factors, cell-cycle regulators, and chemokine receptors. ChIP/seq studies established that KAP1 bound at or close to several of these genes and controlled chromatin status at their promoters. Genome wide, KAP1 binding sites lacked active B cell-specific enhancers and were enriched in repressive histone marks, further supporting a role for this molecule in gene silencing in vivo. Likely responsible for tethering KAP1 to at least some of these targets, a discrete subset of KRAB-ZFPs is enriched in B lymphocytes. Our results therefore reveal the role of KRAB/KAP1-mediated epigenetic regulation in B-cell development and homeostasis.

Demogenomic inference from spatially and temporally heterogeneous samples.

Modern and ancient genomes are not necessarily drawn from homogeneous populations, as they may have been collected from different places and at different times. This heterogeneous sampling can be an issue for demographic inferences and results in biased demographic parameters and incorrect model choice if not properly considered. When explicitly accounted for, it can result in very complex models and high data dimensionality that are difficult to analyse. In this paper, we formally study the impact of such spatial and temporal sampling heterogeneity on demographic inference, and we introduce a way to circumvent this problem. To deal with structured samples without increasing the dimensionality of the site frequency spectrum (SFS), we introduce a new structured approach to the existing program fastsimcoal2. We assess the efficiency and relevance of this methodological update with simulated and modern human genomic data. We particularly focus on spatial and temporal heterogeneities to evidence the interest of this new SFS-based approach, which can be especially useful when handling scattered and ancient DNA samples, as in conservation genetics or archaeogenetics.

Liver-specific ablation of Krüppel-associated box-associated protein 1 in mice leads to male-predominant hepatosteatosis and development of liver adenoma.

UNLABELLED: The liver is characterized by sexually dimorphic gene expression translating into sex-specific differences in lipid, drug, steroid hormone, and xenobiotic metabolism, with distinct responses of males and females to environmental challenges. Here, we investigated the role of the Krüppel-associated box (KRAB)-associated protein 1 (KAP1) epigenetic regulator in this process. Liver-specific KAP1 knockout (KO) led to strikingly sexually dimorphic phenotypic disturbances, including male-predominant steatosis and hepatic tumors with up-regulation of protein kinase B and extracellular signal-related kinases 1/2 mitogen-activated protein kinase signaling. This correlated with the sex-specific transcriptional dysregulation of a wide range of metabolic genes, notably those involved in retinol and sex hormone processing as well as in detoxification. Furthermore, chromatin immunoprecipitation followed by deep sequencing indicated that a number of dysregulated genes are direct targets of the KRAB/KAP1 repression system. Those genes include sexually dimorphic cytochrome P 450 Cyp2d9, glutathione S-transferase π, Cyp2a, Cyp2b, and Cyp3a gene clusters. Additionally, we identified a male-restricted KAP1-binding site in the fat-specific protein 27 gene, correlating with its male-predominant up-regulation upon Kap1 deletion, suggesting that the latter might be an important trigger in the development of male-specific hepatosteatosis and secondary tumorigenesis. CONCLUSION: This work reveals KRAB/KAP1-mediated transcriptional regulation as a central event in metabolic control hormones, drugs, and xenobiotics in the liver and further links disturbances in these processes with hepatic carcinogenesis.

KAP1 regulates gene networks controlling T-cell development and responsiveness.

Chromatin remodeling at specific genomic loci controls lymphoid differentiation. Here, we investigated the role played in this process by Kruppel-associated box (KRAB)-associated protein 1 (KAP1), the universal cofactor of KRAB-zinc finger proteins (ZFPs), a tetrapod-restricted family of transcriptional repressors. T-cell-specific Kap1-deleted mice displayed a significant expansion of immature thymocytes, imbalances in CD4(+)/CD8(+) cell ratios, and altered responses to TCR and TGFß stimulation when compared to littermate KAP1 control mice. Transcriptome and chromatin studies revealed that KAP1 binds T-cell-specific cis-acting regulatory elements marked by the H3K9me3 repressive mark and enriched in Ikaros/NuRD complexes. Also, KAP1 directly controls the expression of several genes involved in TCR and cytokine signaling. Among these, regulation of FoxO1 seems to play a major role in this system. Likely responsible for tethering KAP1 to at least part of its genomic targets, a small number of KRAB-ZFPs are selectively expressed in T-lymphoid cells. These results reveal the so far unsuspected yet important role of KAP1-mediated epigenetic regulation in T-lymphocyte differentiation and activation.

Metagenomics and metabarcoding experimental choices and their impact on microbial community characterization in freshwater recirculating aquaculture systems.

BACKGROUND: Microbial communities in recirculating aquaculture systems (RAS) play a role in system success, nutrient cycling, and water quality. Considering the increasing socio-economic role of fish farming, e.g., regarding food security, an in-depth understanding of aquaculture microbial communities is also relevant from a management perspective, especially regarding the growth, development, and welfare of the farmed animal. However, the current data on the composition of microbial communities within RAS is patchy, which is partly attributable to diverging method choices that render comparative analyses challenging. Therefore, there is a need for accurate, standardized, and user-friendly methods to study microbial communities in aquaculture systems. RESULTS: We compared sequencing approach performances (3 types of 16S short amplicon sequencing, PacBio long-read amplicon sequencing, and amplification-free shotgun metagenomics) in the characterization of microbial communities in two commercial RAS fish farms. Results showed that 16S primer choice and amplicon length affect some values (e.g., diversity measures, number of assigned taxa or distinguishing ASVs) but have no impact on spatio-temporal patterns between sample types, farms and time points. This implies that 16S rRNA approaches are adequate for community studies. The long-read amplicons underperformed regarding the quantitative resolution of spatio-temporal patterns but were suited to identify functional services, e.g., nitrification cycling and the detection of pathogens. Finally, shotgun metagenomics extended the picture to fungi, viruses, and bacteriophages, opening avenues for exploring inter-domain interactions. All sequencing datasets agreed on major prokaryotic players, such as Actinobacteriota, Bacteroidota, Nitrospirota, and Proteobacteria. CONCLUSION: The different sequencing approaches yielded overlapping and highly complementary results, with each contributing unique data not obtainable with the other approaches. We conclude that a tiered approach constitutes a strategy for obtaining the maximum amount of information on aquaculture microbial communities and can inform basic research on community evolution dynamics. For specific and/or applied questions, single-method approaches are more practical and cost-effective and could lead to better farm management practices.

A gene-rich, transcriptionally active environment and the pre-deposition of repressive marks are predictive of susceptibility to KRAB/KAP1-mediated silencing.

BACKGROUND: KRAB-ZFPs (Krüppel-associated box domain-zinc finger proteins) are vertebrate-restricted transcriptional repressors encoded in the hundreds by the mouse and human genomes. They act via an essential cofactor, KAP1, which recruits effectors responsible for the formation of facultative heterochromatin. We have recently shown that KRAB/KAP1 can mediate long-range transcriptional repression through heterochromatin spreading, but also demonstrated that this process is at times countered by endogenous influences. METHOD: To investigate this issue further we used an ectopic KRAB-based repressor. This system allowed us to tether KRAB/KAP1 to hundreds of euchromatic sites within genes, and to record its impact on gene expression. We then correlated this KRAB/KAP1-mediated transcriptional effect to pre-existing genomic and chromatin structures to identify specific characteristics making a gene susceptible to repression. RESULTS: We found that genes that were susceptible to KRAB/KAP1-mediated silencing carried higher levels of repressive histone marks both at the promoter and over the transcribed region than genes that were insensitive. In parallel, we found a high enrichment in euchromatic marks within both the close and more distant environment of these genes. CONCLUSION: Together, these data indicate that high levels of gene activity in the genomic environment and the pre-deposition of repressive histone marks within a gene increase its susceptibility to KRAB/KAP1-mediated repression.

Impact of Genetic Variation in Gene Regulatory Sequences: A Population Genomics Perspective.

The unprecedented rise of high-throughput sequencing and assay technologies has provided a detailed insight into the non-coding sequences and their potential role as gene expression regulators. These regulatory non-coding sequences are also referred to as cis-regulatory elements (CREs). Genetic variants occurring within CREs have been shown to be associated with altered gene expression and phenotypic changes. Such variants are known to occur spontaneously and ultimately get fixed, due to selection and genetic drift, in natural populations and, in some cases, pave the way for speciation. Hence, the study of genetic variation at CREs has improved our overall understanding of the processes of local adaptation and evolution. Recent advances in high-throughput sequencing and better annotations of CREs have enabled the evaluation of the impact of such variation on gene expression, phenotypic alteration and fitness. Here, we review recent research on the evolution of CREs and concentrate on studies that have investigated genetic variation occurring in these regulatory sequences within the context of population genetics.

Demographic analyses of a new sample of haploid genomes from a Swedish population of Drosophila melanogaster.

European and African natural populations of Drosophila melanogaster have been the focus of several studies aiming at inferring demographic and adaptive processes based on genetic variation data. However, in these analyses little attention has been given to gene flow between African and European samples. Here we present a dataset consisting of 14 fully sequenced haploid genomes sampled from a natural population from the northern species range (Umeå, Sweden). We co-analyzed this new data with an African population to compare the likelihood of several competing demographic scenarios for European and African populations and show that gene flow improves the fit of demographic models to data.

Mitochondrial haplotypes affect metabolic phenotypes in the Drosophila Genetic Reference Panel.

The nature and extent of mitochondrial DNA variation in a population and how it affects traits is poorly understood. Here we resequence the mitochondrial genomes of 169 Drosophila Genetic Reference Panel lines, identifying 231 variants that stratify along 12 mitochondrial haplotypes. We identify 1,845 cases of mitonuclear allelic imbalances, thus implying that mitochondrial haplotypes are reflected in the nuclear genome. However, no major fitness effects are associated with mitonuclear imbalance, suggesting that such imbalances reflect population structure at the mitochondrial level rather than genomic incompatibilities. Although mitochondrial haplotypes have no direct impact on mitochondrial respiration, some haplotypes are associated with stress- and metabolism-related phenotypes, including food intake in males. Finally, through reciprocal swapping of mitochondrial genomes, we demonstrate that a mitochondrial haplotype associated with high food intake can rescue a low food intake phenotype. Together, our findings provide new insight into population structure at the mitochondrial level and point to the importance of incorporating mitochondrial haplotypes in genotype-phenotype relationship studies.

The Demographic History of African Drosophila melanogaster.

As one of the most commonly utilized organisms in the study of local adaptation, an accurate characterization of the demographic history of Drosophila melanogaster remains as an important research question. This owes both to the inherent interest in characterizing the population history of this model organism, as well as to the well-established importance of an accurate null demographic model for increasing power and decreasing false positive rates in genomic scans for positive selection. Although considerable attention has been afforded to this issue in non-African populations, less is known about the demographic history of African populations, including from the ancestral range of the species. While qualitative predictions and hypotheses have previously been forwarded, we here present a quantitative model fitting of the population history characterizing both the ancestral Zambian population range as well as the subsequently colonized west African populations, which themselves served as the source of multiple non-African colonization events. We here report the split time of the West African population at 72 kya, a date corresponding to human migration into this region as well as a period of climatic changes in the African continent. Furthermore, we have estimated population sizes at this split time. These parameter estimates thus represent an important null model for future investigations in to African and non-African D. melanogaster populations alike.

The evolution of gene expression and binding specificity of the largest transcription factor family in primates.

The KRAB-containing zinc finger (KRAB-ZF) proteins represent the largest family of transcription factors (TFs) in humans, yet for the great majority, their function and specific genomic target remain unknown. However, it has been shown that a large fraction of these genes arose from segmental duplications, and that they have expanded in gene and zinc finger number throughout vertebrate evolution. To determine whether this expansion is linked to selective pressures acting on different domains, we have manually curated all KRAB-ZF genes present in the human genome together with their orthologous genes in three closely related species and assessed the evolutionary forces acting at the sequence level as well as on their expression profiles. We provide evidence that KRAB-ZFs can be separated into two categories according to the polymorphism present in their DNA-contacting residues. Those carrying a nonsynonymous single nucleotide polymorphism (SNP) in their DNA-contacting amino acids exhibit significantly reduced expression in all tissues, have emerged in a recent lineage, and seem to be less strongly constrained evolutionarily than those without such a polymorphism. This work provides evidence for a link between age of the TF, as well as polymorphism in their DNA-contacting residues and expression levels-both of which may be jointly affected by selection.

Release of human cytomegalovirus from latency by a KAP1/TRIM28 phosphorylation switch.

Human cytomegalovirus (HCMV) is a highly prevalent pathogen that induces life-long infections notably through the establishment of latency in hematopoietic stem cells (HSC). Bouts of reactivation are normally controlled by the immune system, but can be fatal in immuno-compromised individuals such as organ transplant recipients. Here, we reveal that HCMV latency in human CD34(+) HSC reflects the recruitment on the viral genome of KAP1, a master co-repressor, together with HP1 and the SETDB1 histone methyltransferase, which results in transcriptional silencing. During lytic infection, KAP1 is still associated with the viral genome, but its heterochromatin-inducing activity is suppressed by mTOR-mediated phosphorylation. Correspondingly, HCMV can be forced out of latency by KAP1 knockdown or pharmacological induction of KAP1 phosphorylation, and this process can be potentiated by activating NFkB with TNF-α. These results suggest new approaches both to curtail CMV infection and to purge the virus from organ transplants.