Causes and Consequences of Purifying Selection on SARS-CoV-2.

Owing to a lag between a deleterious mutation's appearance and its selective removal, gold-standard methods for mutation rate estimation assume no meaningful loss of mutations between parents and offspring. Indeed, from analysis of closely related lineages, in SARS-CoV-2, the Ka/Ks ratio was previously estimated as 1.008, suggesting no within-host selection. By contrast, we find a higher number of observed SNPs at 4-fold degenerate sites than elsewhere and, allowing for the virus's complex mutational and compositional biases, estimate that the mutation rate is at least 49-67% higher than would be estimated based on the rate of appearance of variants in sampled genomes. Given the high Ka/Ks one might assume that the majority of such intrahost selection is the purging of nonsense mutations. However, we estimate that selection against nonsense mutations accounts for only ∼10% of all the "missing" mutations. Instead, classical protein-level selective filters (against chemically disparate amino acids and those predicted to disrupt protein functionality) account for many missing mutations. It is less obvious why for an intracellular parasite, amino acid cost parameters, notably amino acid decay rate, is also significant. Perhaps most surprisingly, we also find evidence for real-time selection against synonymous mutations that move codon usage away from that of humans. We conclude that there is common intrahost selection on SARS-CoV-2 that acts on nonsense, missense, and possibly synonymous mutations. This has implications for methods of mutation rate estimation, for determining times to common ancestry and the potential for intrahost evolution including vaccine escape.

Transcription, mRNA Export, and Immune Evasion Shape the Codon Usage of Viruses.

The nucleotide composition, dinucleotide composition, and codon usage of many viruses differ from their hosts. These differences arise because viruses are subject to unique mutation and selection pressures that do not apply to host genomes; however, the molecular mechanisms that underlie these evolutionary forces are unclear. Here, we analyzed the patterns of codon usage in 1,520 vertebrate-infecting viruses, focusing on parameters known to be under selection and associated with gene regulation. We find that GC content, dinucleotide content, and splicing and m6A modification-related sequence motifs are associated with the type of genetic material (DNA or RNA), strandedness, and replication compartment of viruses. In an experimental follow-up, we find that the effects of GC content on gene expression depend on whether the genetic material is delivered to the cell as DNA or mRNA, whether it is transcribed by endogenous or exogenous RNA polymerase, and whether transcription takes place in the nucleus or cytoplasm. Our results suggest that viral codon usage cannot be explained by a simple adaptation to the codon usage of the host-instead, it reflects the combination of multiple selective and mutational pressures, including the need for efficient transcription, export, and immune evasion.

Evidence for Strong Mutation Bias toward, and Selection against, U Content in SARS-CoV-2: Implications for Vaccine Design.

Large-scale re-engineering of synonymous sites is a promising strategy to generate vaccines either through synthesis of attenuated viruses or via codon-optimized genes in DNA vaccines. Attenuation typically relies on deoptimization of codon pairs and maximization of CpG dinucleotide frequencies. So as to formulate evolutionarily informed attenuation strategies that aim to force nucleotide usage against the direction favored by selection, here, we examine available whole-genome sequences of SARS-CoV-2 to infer patterns of mutation and selection on synonymous sites. Analysis of mutational profiles indicates a strong mutation bias toward U. In turn, analysis of observed synonymous site composition implicates selection against U. Accounting for dinucleotide effects reinforces this conclusion, observed UU content being a quarter of that expected under neutrality. Possible mechanisms of selection against U mutations include selection for higher expression, for high mRNA stability or lower immunogenicity of viral genes. Consistent with gene-specific selection against CpG dinucleotides, we observe systematic differences of CpG content between SARS-CoV-2 genes. We propose an evolutionarily informed approach to attenuation that, unusually, seeks to increase usage of the already most common synonymous codons. Comparable analysis of H1N1 and Ebola finds that GC3 deviated from neutral equilibrium is not a universal feature, cautioning against generalization of results.

Deficient pregnenolone synthesis associated with congenital adrenal hyperplasia and organelle dysfunction.

Steroid hormones are essential for the survival of all mammals. In adrenal glands and gonads, cytochrome P450 side chain cleavage enzyme (SCC or CYP11A1), catalyzes conversion of cholesterol to pregnenolone. We studied a patient with ambiguous genitalia by the absence of Müllerian ducts and the presence of an incompletely formed vagina, who had extremely high adrenocorticotropic hormone (ACTH) and reduced pregnenolone levels with enlarged adrenal glands. The testes revealed seminiferous tubules, stroma, rete testis with interstitial fibrosis and reduced number of germ cells. Electron microscopy showed that the patient's testicular mitochondrial size was small with little SCC expression within the mitochondria. The mitochondria were not close to the mitochondria-associated ER membrane (MAM), and cells were filled with the microfilaments. Our result revealed that absence of pregnenolone is associated with organelle stress, leading to altered protein organization that likely created steric hindrance in testicular cells. Learning points: Testes revealed seminiferous tubules, stroma, rete testis with interstitial fibrosis and reduced number of germ cells; Testicular mitochondrial size was small with little SCC expression within the mitochondria; Absence of pregnenolone is associated with organelle stress.

A Family of Vertebrate-Specific Polycombs Encoded by the LCOR/LCORL Genes Balance PRC2 Subtype Activities.

The polycomb repressive complex 2 (PRC2) consists of core subunits SUZ12, EED, RBBP4/7, and EZH1/2 and is responsible for mono-, di-, and tri-methylation of lysine 27 on histone H3. Whereas two distinct forms exist, PRC2.1 (containing one polycomb-like protein) and PRC2.2 (containing AEBP2 and JARID2), little is known about their differential functions. Here, we report the discovery of a family of vertebrate-specific PRC2.1 proteins, "PRC2 associated LCOR isoform 1" (PALI1) and PALI2, encoded by the LCOR and LCORL gene loci, respectively. PALI1 promotes PRC2 methyltransferase activity in vitro and in vivo and is essential for mouse development. Pali1 and Aebp2 define mutually exclusive, antagonistic PRC2 subtypes that exhibit divergent H3K27-tri-methylation activities. The balance of these PRC2.1/PRC2.2 activities is required for the appropriate regulation of polycomb target genes during differentiation. PALI1/2 potentially link polycombs with transcriptional co-repressors in the regulation of cellular identity during development and in cancer.

Dosage-sensitive genes in evolution and disease.

For a subset of genes in our genome a change in gene dosage, by duplication or deletion, causes a phenotypic effect. These dosage-sensitive genes may confer an advantage upon copy number change, but more typically they are associated with disease, including heart disease, cancers and neuropsychiatric disorders. This gene copy number sensitivity creates characteristic evolutionary constraints that can serve as a diagnostic to identify dosage-sensitive genes. Though the link between copy number change and disease is well-established, the mechanism of pathogenicity is usually opaque. We propose that gene expression level may provide a common basis for the pathogenic effects of many copy number variants.

De novo disruption of promoter and exon 1 of STAR gene reveals essential role for gonadal development.

SUMMARY: Cholesterol transport into the mitochondria is required for synthesis of the first steroid, pregnenolone. Cholesterol is transported by the steroidogenic acute regulatory protein (STAR), which acts at the outer mitochondrial membrane prior to its import. Mutations in the STAR protein result in lipoid congenital adrenal hyperplasia (CAH). Although the STAR protein consists of seven exons, biochemical analysis in nonsteroidogenic COS-1 cells showed that the first two were not essential for pregnenolone synthesis. Here, we present a patient with ambiguous genitalia, salt-lossing crisis within two weeks after birth and low cortisol levels. Sequence analysis of the STAR, including the exon-intron boundaries, showed the complete deletion of exon 1 as well as more than 50 nucleotides upstream of STAR promoter. Mitochondrial protein import with the translated protein through synthesis cassette of the mutant STAR lacking exon 1 showed protein translation, but it is less likely to have synthesized without a promoter in our patient. Thus, a full-length STAR gene is necessary for physiological mitochondrial cholesterol transport in vivo. LEARNING POINTS: STAR exon 1 deletion caused lipoid CAH.Exon 1 substitution does not affect biochemical activity.StAR promoter is responsible for gonadal development.

Dosage sensitivity is a major determinant of human copy number variant pathogenicity.

Human copy number variants (CNVs) account for genome variation an order of magnitude larger than single-nucleotide polymorphisms. Although much of this variation has no phenotypic consequences, some variants have been associated with disease, in particular neurodevelopmental disorders. Pathogenic CNVs are typically very large and contain multiple genes, and understanding the cause of the pathogenicity remains a major challenge. Here we show that pathogenic CNVs are significantly enriched for genes involved in development and genes that have greater evolutionary copy number conservation across mammals, indicative of functional constraints. Conversely, genes found in benign CNV regions have more variable copy number. These evolutionary constraints are characteristic of genes in pathogenic CNVs and can only be explained by dosage sensitivity of those genes. These results implicate dosage sensitivity of individual genes as a common cause of CNV pathogenicity. These evolutionary metrics suggest a path to identifying disease genes in pathogenic CNVs.

Novel SCC mutation in a patient of Mexican descent with sex reversal, salt-losing crisis and adrenal failure.

Congenital adrenal hyperplasia (CAH) is caused by mutations in cytochrome P450 side chain cleavage enzyme (CYP11A1 and old name, SCC). Errors in cholesterol side chain cleavage by the mitochondrial resident CYP11A1 results in an inadequate amount of pregnenolone production. This study was performed to evaluate the cause of salt-losing crisis and possible adrenal failure in a pediatric patient whose mother had a history of two previous stillbirths and loss of another baby within a week of birth. CAH can appear in any population in any region of the world. The study was conducted at Memorial University Medical Center and Mercer University School of Medicine. The patient was admitted to Pediatric Endocrinology Clinic due to salt-losing crisis and possible adrenal failure. The patient had CAH, an autosomal recessive disease, due to a novel mutation in exon 5 of the CYP11A1 gene, which generated a truncated protein of 286 amino acids compared with wild-type protein that has 521 amino acids (W286X). Although unrelated, both parents are carriers. Mitochondrial protein import analysis of the mutant CYP11A1 in steroidogenic MA-10 cells showed that the protein is imported in a similar fashion as observed for the wild-type protein and was cleaved to a shorter fragment. However, mutant's activity was 10% of that obtained for the wild-type protein in non-steroidogenic COS-1 cells. In a patient of Mexican descent, a homozygous CYP11A1 mutation caused CAH, suggesting that this disease is not geographically restricted even in a homogeneous population. LEARNING POINTS: Novel mutation in CYP11A1 causes CAH;This is a pure population from Central Mexico;Novel mutation created early truncated protein.

A chromatin-independent role of Polycomb-like 1 to stabilize p53 and promote cellular quiescence.

Polycomb-like proteins 1-3 (PCL1-3) are substoichiometric components of the Polycomb-repressive complex 2 (PRC2) that are essential for association of the complex with chromatin. However, it remains unclear why three proteins with such apparent functional redundancy exist in mammals. Here we characterize their divergent roles in both positively and negatively regulating cellular proliferation. We show that while PCL2 and PCL3 are E2F-regulated genes expressed in proliferating cells, PCL1 is a p53 target gene predominantly expressed in quiescent cells. Ectopic expression of any PCL protein recruits PRC2 to repress the INK4A gene; however, only PCL2 and PCL3 confer an INK4A-dependent proliferative advantage. Remarkably, PCL1 has evolved a PRC2- and chromatin-independent function to negatively regulate proliferation. We show that PCL1 binds to and stabilizes p53 to induce cellular quiescence. Moreover, depletion of PCL1 phenocopies the defects in maintaining cellular quiescence associated with p53 loss. This newly evolved function is achieved by the binding of the PCL1 N-terminal PHD domain to the C-terminal domain of p53 through two unique serine residues, which were acquired during recent vertebrate evolution. This study illustrates the functional bifurcation of PCL proteins, which act in both a chromatin-dependent and a chromatin-independent manner to regulate the INK4A and p53 pathways.