Unfinished Business: Evolution of the MHC and the Adaptive Immune System of Jawed Vertebrates.

The major histocompatibility complex (MHC) is a large genetic region with many genes, including the highly polymorphic classical class I and II genes that play crucial roles in adaptive as well as innate immune responses. The organization of the MHC varies enormously among jawed vertebrates, but class I and II genes have not been found in other animals. How did the MHC arise, and are there underlying principles that can help us to understand the evolution of the MHC? This review considers what it means to be an MHC and the potential importance of genome-wide duplication, gene linkage, and gene coevolution for the emergence and evolution of an adaptive immune system. Then it considers what the original antigen-specific receptor and MHC molecule might have looked like, how peptide binding might have evolved, and finally the importance of adaptive immunity in general.

Genetics of Infectious and Inflammatory Diseases: Overlapping Discoveries from Association and Exome-Sequencing Studies.

Genome technologies have defined a complex genetic architecture in major infectious, inflammatory, and autoimmune disorders. High density marker arrays and Immunochips have powered genome-wide association studies (GWAS) that have mapped nearly 450 genetic risk loci in 22 major inflammatory diseases, including a core of common genes that play a central role in pathological inflammation. Whole-exome and whole-genome sequencing have identified more than 265 genes in which mutations cause primary immunodeficiencies and rare forms of severe inflammatory bowel disease. Combined analysis of inflammatory disease GWAS and primary immunodeficiencies point to shared proteins and pathways that are required for immune cell development and protection against infections and are also associated with pathological inflammation. Finally, sequencing of chromatin immunoprecipitates containing specific transcription factors, with parallel RNA sequencing, has charted epigenetic regulation of gene expression by proinflammatory transcription factors in immune cells, providing complementary information to characterize morbid genes at infectious and inflammatory disease loci.

A multi-cohort genome-wide association study in African ancestry individuals reveals risk loci for primary open-angle glaucoma.

Primary open-angle glaucoma (POAG), the leading cause of irreversible blindness worldwide, disproportionately affects individuals of African ancestry. We conducted a genome-wide association study (GWAS) for POAG in 11,275 individuals of African ancestry (6,003 cases; 5,272 controls). We detected 46 risk loci associated with POAG at genome-wide significance. Replication and post-GWAS analyses, including functionally informed fine-mapping, multiple trait co-localization, and in silico validation, implicated two previously undescribed variants (rs1666698 mapping to DBF4P2; rs34957764 mapping to ROCK1P1) and one previously associated variant (rs11824032 mapping to ARHGEF12) as likely causal. For individuals of African ancestry, a polygenic risk score (PRS) for POAG from our mega-analysis (African ancestry individuals) outperformed a PRS from summary statistics of a much larger GWAS derived from European ancestry individuals. This study quantifies the genetic architecture similarities and differences between African and non-African ancestry populations for this blinding disease.

Complex effects of sequence variants on lipid levels and coronary artery disease.

Many sequence variants have additive effects on blood lipid levels and, through that, on the risk of coronary artery disease (CAD). We show that variants also have non-additive effects and interact to affect lipid levels as well as affecting variance and correlations. Variance and correlation effects are often signatures of epistasis or gene-environmental interactions. These complex effects can translate into CAD risk. For example, Trp154Ter in FUT2 protects against CAD among subjects with the A1 blood group, whereas it associates with greater risk of CAD in others. His48Arg in ADH1B interacts with alcohol consumption to affect lipid levels and CAD. The effect of variants in TM6SF2 on blood lipids is greatest among those who never eat oily fish but absent from those who often do. This work demonstrates that variants that affect variance of quantitative traits can allow for the discovery of epistasis and interactions of variants with the environment.

An evolved AAV variant enables efficient genetic engineering of murine T cells.

Precise targeting of large transgenes to T cells using homology-directed repair has been transformative for adoptive cell therapies and T cell biology. Delivery of DNA templates via adeno-associated virus (AAV) has greatly improved knockin efficiencies, but the tropism of current AAV serotypes restricts their use to human T cells employed in immunodeficient mouse models. To enable targeted knockins in murine T cells, we evolved Ark313, a synthetic AAV that exhibits high transduction efficiency in murine T cells. We performed a genome-wide knockout screen and identified QA2 as an essential factor for Ark313 infection. We demonstrate that Ark313 can be used for nucleofection-free DNA delivery, CRISPR-Cas9-mediated knockouts, and targeted integration of large transgenes. Ark313 enables preclinical modeling of Trac-targeted CAR-T and transgenic TCR-T cells in immunocompetent models. Efficient gene targeting in murine T cells holds great potential for improved cell therapies and opens avenues in experimental T cell immunology.

Limb development genes underlie variation in human fingerprint patterns.

Fingerprints are of long-standing practical and cultural interest, but little is known about the mechanisms that underlie their variation. Using genome-wide scans in Han Chinese cohorts, we identified 18 loci associated with fingerprint type across the digits, including a genetic basis for the long-recognized "pattern-block" correlations among the middle three digits. In particular, we identified a variant near EVI1 that alters regulatory activity and established a role for EVI1 in dermatoglyph patterning in mice. Dynamic EVI1 expression during human development supports its role in shaping the limbs and digits, rather than influencing skin patterning directly. Trans-ethnic meta-analysis identified 43 fingerprint-associated loci, with nearby genes being strongly enriched for general limb development pathways. We also found that fingerprint patterns were genetically correlated with hand proportions. Taken together, these findings support the key role of limb development genes in influencing the outcome of fingerprint patterning.

Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations.

Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation.

Identification of Required Host Factors for SARS-CoV-2 Infection in Human Cells.

To better understand host-virus genetic dependencies and find potential therapeutic targets for COVID-19, we performed a genome-scale CRISPR loss-of-function screen to identify host factors required for SARS-CoV-2 viral infection of human alveolar epithelial cells. Top-ranked genes cluster into distinct pathways, including the vacuolar ATPase proton pump, Retromer, and Commander complexes. We validate these gene targets using several orthogonal methods such as CRISPR knockout, RNA interference knockdown, and small-molecule inhibitors. Using single-cell RNA-sequencing, we identify shared transcriptional changes in cholesterol biosynthesis upon loss of top-ranked genes. In addition, given the key role of the ACE2 receptor in the early stages of viral entry, we show that loss of RAB7A reduces viral entry by sequestering the ACE2 receptor inside cells. Overall, this work provides a genome-scale, quantitative resource of the impact of the loss of each host gene on fitness/response to viral infection.

A Genome-wide ER-phagy Screen Highlights Key Roles of Mitochondrial Metabolism and ER-Resident UFMylation.

Selective autophagy of organelles is critical for cellular differentiation, homeostasis, and organismal health. Autophagy of the ER (ER-phagy) is implicated in human neuropathy but is poorly understood beyond a few autophagosomal receptors and remodelers. By using an ER-phagy reporter and genome-wide CRISPRi screening, we identified 200 high-confidence human ER-phagy factors. Two pathways were unexpectedly required for ER-phagy. First, reduced mitochondrial metabolism represses ER-phagy, which is opposite of general autophagy and is independent of AMPK. Second, ER-localized UFMylation is required for ER-phagy to repress the unfolded protein response via IRE1α. The UFL1 ligase is brought to the ER surface by DDRGK1 to UFMylate RPN1 and RPL26 and preferentially targets ER sheets for degradation, analogous to PINK1-Parkin regulation during mitophagy. Our data provide insight into the cellular logic of ER-phagy, reveal parallels between organelle autophagies, and provide an entry point to the relatively unexplored process of degrading the ER network.

Genomic History of Neolithic to Bronze Age Anatolia, Northern Levant, and Southern Caucasus.

Here, we report genome-wide data analyses from 110 ancient Near Eastern individuals spanning the Late Neolithic to Late Bronze Age, a period characterized by intense interregional interactions for the Near East. We find that 6th millennium BCE populations of North/Central Anatolia and the Southern Caucasus shared mixed ancestry on a genetic cline that formed during the Neolithic between Western Anatolia and regions in today's Southern Caucasus/Zagros. During the Late Chalcolithic and/or the Early Bronze Age, more than half of the Northern Levantine gene pool was replaced, while in the rest of Anatolia and the Southern Caucasus, we document genetic continuity with only transient gene flow. Additionally, we reveal a genetically distinct individual within the Late Bronze Age Northern Levant. Overall, our study uncovers multiple scales of population dynamics through time, from extensive admixture during the Neolithic period to long-distance mobility within the globalized societies of the Late Bronze Age. VIDEO ABSTRACT.

Per-Nucleus Crossover Covariation and Implications for Evolution.

Crossing over is a nearly universal feature of sexual reproduction. Here, analysis of crossover numbers on a per-chromosome and per-nucleus basis reveals a fundamental, evolutionarily conserved feature of meiosis: within individual nuclei, crossover frequencies covary across different chromosomes. This effect results from per-nucleus covariation of chromosome axis lengths. Crossovers can promote evolutionary adaptation. However, the benefit of creating favorable new allelic combinations must outweigh the cost of disrupting existing favorable combinations. Covariation concomitantly increases the frequencies of gametes with especially high, or especially low, numbers of crossovers, and thus might concomitantly enhance the benefits of crossing over while reducing its costs. A four-locus population genetic model suggests that such an effect can pertain in situations where the environment fluctuates: hyper-crossover gametes are advantageous when the environment changes while hypo-crossover gametes are advantageous in periods of environmental stasis. These findings reveal a new feature of the basic meiotic program and suggest a possible adaptive advantage.

Root System Depth in Arabidopsis Is Shaped by EXOCYST70A3 via the Dynamic Modulation of Auxin Transport.

Root system architecture (RSA), the distribution of roots in soil, plays a major role in plant survival. RSA is shaped by multiple developmental processes that are largely governed by the phytohormone auxin, suggesting that auxin regulates responses of roots that are important for local adaptation. However, auxin has a central role in numerous processes, and it is unclear which molecular mechanisms contribute to the variation in RSA for environmental adaptation. Using natural variation in Arabidopsis, we identify EXOCYST70A3 as a modulator of the auxin system that causes variation in RSA by acting on PIN4 protein distribution. Allelic variation and genetic perturbation of EXOCYST70A3 lead to alteration of root gravitropic responses, resulting in a different RSA depth profile and drought resistance. Overall our findings suggest that the local modulation of the pleiotropic auxin pathway can gives rise to distinct RSAs that can be adaptive in specific environments.

Intra- and Inter-cellular Rewiring of the Human Colon during Ulcerative Colitis.

Genome-wide association studies (GWAS) have revealed risk alleles for ulcerative colitis (UC). To understand their cell type specificities and pathways of action, we generate an atlas of 366,650 cells from the colon mucosa of 18 UC patients and 12 healthy individuals, revealing 51 epithelial, stromal, and immune cell subsets, including BEST4+ enterocytes, microfold-like cells, and IL13RA2+IL11+ inflammatory fibroblasts, which we associate with resistance to anti-TNF treatment. Inflammatory fibroblasts, inflammatory monocytes, microfold-like cells, and T cells that co-express CD8 and IL-17 expand with disease, forming intercellular interaction hubs. Many UC risk genes are cell type specific and co-regulated within relatively few gene modules, suggesting convergence onto limited sets of cell types and pathways. Using this observation, we nominate and infer functions for specific risk genes across GWAS loci. Our work provides a framework for interrogating complex human diseases and mapping risk variants to cell types and pathways.

Genomic Analyses from Non-invasive Prenatal Testing Reveal Genetic Associations, Patterns of Viral Infections, and Chinese Population History.

We analyze whole-genome sequencing data from 141,431 Chinese women generated for non-invasive prenatal testing (NIPT). We use these data to characterize the population genetic structure and to investigate genetic associations with maternal and infectious traits. We show that the present day distribution of alleles is a function of both ancient migration and very recent population movements. We reveal novel phenotype-genotype associations, including several replicated associations with height and BMI, an association between maternal age and EMB, and between twin pregnancy and NRG1. Finally, we identify a unique pattern of circulating viral DNA in plasma with high prevalence of hepatitis B and other clinically relevant maternal infections. A GWAS for viral infections identifies an exceptionally strong association between integrated herpesvirus 6 and MOV10L1, which affects piwi-interacting RNA (piRNA) processing and PIWI protein function. These findings demonstrate the great value and potential of accumulating NIPT data for worldwide medical and genetic analyses.

Genome-wide CRISPR Screens in Primary Human T Cells Reveal Key Regulators of Immune Function.

Human T cells are central effectors of immunity and cancer immunotherapy. CRISPR-based functional studies in T cells could prioritize novel targets for drug development and improve the design of genetically reprogrammed cell-based therapies. However, large-scale CRISPR screens have been challenging in primary human cells. We developed a new method, single guide RNA (sgRNA) lentiviral infection with Cas9 protein electroporation (SLICE), to identify regulators of stimulation responses in primary human T cells. Genome-wide loss-of-function screens identified essential T cell receptor signaling components and genes that negatively tune proliferation following stimulation. Targeted ablation of individual candidate genes characterized hits and identified perturbations that enhanced cancer cell killing. SLICE coupled with single-cell RNA sequencing (RNA-seq) revealed signature stimulation-response gene programs altered by key genetic perturbations. SLICE genome-wide screening was also adaptable to identify mediators of immunosuppression, revealing genes controlling responses to adenosine signaling. The SLICE platform enables unbiased discovery and characterization of functional gene targets in primary cells.

A Genetic Variant Associated with Five Vascular Diseases Is a Distal Regulator of Endothelin-1 Gene Expression.

Genome-wide association studies (GWASs) implicate the PHACTR1 locus (6p24) in risk for five vascular diseases, including coronary artery disease, migraine headache, cervical artery dissection, fibromuscular dysplasia, and hypertension. Through genetic fine mapping, we prioritized rs9349379, a common SNP in the third intron of the PHACTR1 gene, as the putative causal variant. Epigenomic data from human tissue revealed an enhancer signature at rs9349379 exclusively in aorta, suggesting a regulatory function for this SNP in the vasculature. CRISPR-edited stem cell-derived endothelial cells demonstrate rs9349379 regulates expression of endothelin 1 (EDN1), a gene located 600 kb upstream of PHACTR1. The known physiologic effects of EDN1 on the vasculature may explain the pattern of risk for the five associated diseases. Overall, these data illustrate the integration of genetic, phenotypic, and epigenetic analysis to identify the biologic mechanism by which a common, non-coding variant can distally regulate a gene and contribute to the pathogenesis of multiple vascular diseases.

A Natural Allele of a Transcription Factor in Rice Confers Broad-Spectrum Blast Resistance.

Rice feeds half the world's population, and rice blast is often a destructive disease that results in significant crop loss. Non-race-specific resistance has been more effective in controlling crop diseases than race-specific resistance because of its broad spectrum and durability. Through a genome-wide association study, we report the identification of a natural allele of a C2H2-type transcription factor in rice that confers non-race-specific resistance to blast. A survey of 3,000 sequenced rice genomes reveals that this allele exists in 10% of rice, suggesting that this favorable trait has been selected through breeding. This allele causes a single nucleotide change in the promoter of the bsr-d1 gene, which results in reduced expression of the gene through the binding of the repressive MYB transcription factor and, consequently, an inhibition of H2O2 degradation and enhanced disease resistance. Our discovery highlights this novel allele as a strategy for breeding durable resistance in rice.

A Genome-wide CRISPR Screen in Toxoplasma Identifies Essential Apicomplexan Genes.

Apicomplexan parasites are leading causes of human and livestock diseases such as malaria and toxoplasmosis, yet most of their genes remain uncharacterized. Here, we present the first genome-wide genetic screen of an apicomplexan. We adapted CRISPR/Cas9 to assess the contribution of each gene from the parasite Toxoplasma gondii during infection of human fibroblasts. Our analysis defines ∼200 previously uncharacterized, fitness-conferring genes unique to the phylum, from which 16 were investigated, revealing essential functions during infection of human cells. Secondary screens identify as an invasion factor the claudin-like apicomplexan microneme protein (CLAMP), which resembles mammalian tight-junction proteins and localizes to secretory organelles, making it critical to the initiation of infection. CLAMP is present throughout sequenced apicomplexan genomes and is essential during the asexual stages of the malaria parasite Plasmodium falciparum. These results provide broad-based functional information on T. gondii genes and will facilitate future approaches to expand the horizon of antiparasitic interventions.

Complementary strategies for directing in vivo transcription factor binding through DNA binding domains and intrinsically disordered regions.

DNA binding domains (DBDs) of transcription factors (TFs) recognize DNA sequence motifs that are highly abundant in genomes. Within cells, TFs bind a subset of motif-containing sites as directed by either their DBDs or DBD-external (nonDBD) sequences. To define the relative roles of DBDs and nonDBDs in directing binding preferences, we compared the genome-wide binding of 48 (∼30%) budding yeast TFs with their DBD-only, nonDBD-truncated, and nonDBD-only mutants. With a few exceptions, binding locations differed between DBDs and TFs, resulting from the cumulative action of multiple determinants mapped mostly to disordered nonDBD regions. Furthermore, TFs' preferences for promoters of the fuzzy nucleosome architecture were lost in DBD-only mutants, whose binding spread across promoters, implicating nonDBDs' preferences in this hallmark of budding yeast regulatory design. We conclude that DBDs and nonDBDs employ complementary DNA-targeting strategies, whose balance defines TF binding specificity along genomes.

Negative DNA supercoiling induces genome-wide Cas9 off-target activity.

CRISPR-Cas9 is a powerful gene-editing technology; however, off-target activity remains an important consideration for therapeutic applications. We have previously shown that force-stretching DNA induces off-target activity and hypothesized that distortions of the DNA topology in vivo, such as negative DNA supercoiling, could reduce Cas9 specificity. Using single-molecule optical-tweezers, we demonstrate that negative supercoiling λ-DNA induces sequence-specific Cas9 off-target binding at multiple sites, even at low forces. Using an adapted CIRCLE-seq approach, we detect over 10,000 negative-supercoiling-induced Cas9 off-target double-strand breaks genome-wide caused by increased mismatch tolerance. We further demonstrate in vivo that directed local DNA distortion increases off-target activity in cells and that induced off-target events can be detected during Cas9 genome editing. These data demonstrate that Cas9 off-target activity is regulated by DNA topology in vitro and in vivo, suggesting that cellular processes, such as transcription and replication, could induce off-target activity at previously overlooked sites.