Genetic determinants of blood pressure and heart rate identified through ENU-induced mutagenesis with automated meiotic mapping.

We used N-ethyl-N-nitrosurea-induced germline mutagenesis combined with automated meiotic mapping to identify specific systolic blood pressure (SBP) and heart rate (HR) determinant loci. We analyzed 43,627 third-generation (G3) mice from 841 pedigrees to assess the effects of 45,378 variant alleles within 15,760 genes, in both heterozygous and homozygous states. We comprehensively tested 23% of all protein-encoding autosomal genes and found 87 SBP and 144 HR (with 7 affecting both) candidates exhibiting detectable hypomorphic characteristics. Unexpectedly, only 18 of the 87 SBP genes were previously known, while 26 of the 144 genes linked to HR were previously identified. Furthermore, we confirmed the influence of two genes on SBP regulation and three genes on HR control through reverse genetics. This underscores the importance of our research in uncovering genes associated with these critical cardiovascular risk factors and illustrate the effectiveness of germline mutagenesis for defining key determinants of polygenic phenotypes that must be studied in an intact organism.

Obesity caused by an OVOL2 mutation reveals dual roles of OVOL2 in promoting thermogenesis and limiting white adipogenesis.

Using random germline mutagenesis in mice, we identified a viable hypomorphic allele (boh) of the transcription-factor-encoding gene Ovol2 that resulted in obesity, which initially developed with normal food intake and physical activity but decreased energy expenditure. Fat weight was dramatically increased, while lean weight was reduced in 12-week-old boh homozygous mice, culminating by 24 weeks in massive obesity, hepatosteatosis, insulin resistance, and diabetes. The Ovol2boh/boh genotype augmented obesity in Lepob/ob mice, and pair-feeding failed to normalize obesity in Ovol2boh/boh mice. OVOL2-deficient mice were extremely cold intolerant. OVOL2 is essential for brown/beige adipose tissue-mediated thermogenesis. In white adipose tissues, OVOL2 limited adipogenesis by blocking C/EBPα engagement of its transcriptional targets. Overexpression of OVOL2 in adipocytes of mice fed with a high-fat diet reduced total body and liver fat and improved insulin sensitivity. Our data reveal that OVOL2 plays dual functions in thermogenesis and adipogenesis to maintain energy balance.

De novo germline mutation in the dual specificity phosphatase 10 gene accelerates autoimmune diabetes.

Insulin-dependent or type 1 diabetes (T1D) is a polygenic autoimmune disease. In humans, more than 60 loci carrying common variants that confer disease susceptibility have been identified by genome-wide association studies, with a low individual risk contribution for most variants excepting those of the major histocompatibility complex (MHC) region (40 to 50% of risk); hence the importance of missing heritability due in part to rare variants. Nonobese diabetic (NOD) mice recapitulate major features of the human disease including genetic aspects with a key role for the MHC haplotype and a series of Idd loci. Here we mapped in NOD mice rare variants arising from genetic drift and significantly impacting disease risk. To that aim we established by selective breeding two sublines of NOD mice from our inbred NOD/Nck colony exhibiting a significant difference in T1D incidence. Whole-genome sequencing of high (H)- and low (L)-incidence sublines (NOD/NckH and NOD/NckL) revealed a limited number of subline-specific variants. Treating age of diabetes onset as a quantitative trait in automated meiotic mapping (AMM), enhanced susceptibility in NOD/NckH mice was unambiguously attributed to a recessive missense mutation of Dusp10, which encodes a dual specificity phosphatase. The causative effect of the mutation was verified by targeting Dusp10 with CRISPR-Cas9 in NOD/NckL mice, a manipulation that significantly increased disease incidence. The Dusp10 mutation resulted in islet cell down-regulation of type I interferon signature genes, which may exert protective effects against autoimmune aggression. De novo mutations akin to rare human susceptibility variants can alter the T1D phenotype.

Thousands of induced germline mutations affecting immune cells identified by automated meiotic mapping coupled with machine learning.

Forward genetic studies use meiotic mapping to adduce evidence that a particular mutation, normally induced by a germline mutagen, is causative of a particular phenotype. Particularly in small pedigrees, cosegregation of multiple mutations, occasional unawareness of mutations, and paucity of homozygotes may lead to erroneous declarations of cause and effect. We sought to improve the identification of mutations causing immune phenotypes in mice by creating Candidate Explorer (CE), a machine-learning software program that integrates 67 features of genetic mapping data into a single numeric score, mathematically convertible to the probability of verification of any putative mutation-phenotype association. At this time, CE has evaluated putative mutation-phenotype associations arising from screening damaging mutations in ∼55% of mouse genes for effects on flow cytometry measurements of immune cells in the blood. CE has therefore identified more than half of genes within which mutations can be causative of flow cytometric phenovariation in Mus musculus The majority of these genes were not previously known to support immune function or homeostasis. Mouse geneticists will find CE data informative in identifying causative mutations within quantitative trait loci, while clinical geneticists may use CE to help connect causative variants with rare heritable diseases of immunity, even in the absence of linkage information. CE displays integrated mutation, phenotype, and linkage data, and is freely available for query online.

SLFN2 protection of tRNAs from stress-induced cleavage is essential for T cell-mediated immunity.

Reactive oxygen species (ROS) increase in activated T cells because of metabolic activity induced to support T cell proliferation and differentiation. We show that these ROS trigger an oxidative stress response that leads to translation repression. This response is countered by Schlafen 2 (SLFN2), which directly binds transfer RNAs (tRNAs) to protect them from cleavage by the ribonuclease angiogenin. T cell-specific SLFN2 deficiency results in the accumulation of tRNA fragments, which inhibit translation and promote stress-granule formation. Interleukin-2 receptor ß (IL-2Rß) and IL-2Rγ fail to be translationally up-regulated after T cell receptor stimulation, rendering SLFN2-deficient T cells insensitive to interleukin-2's mitogenic effects. SLFN2 confers resistance against the ROS-mediated translation-inhibitory effects of oxidative stress normally induced by T cell activation, permitting the robust protein synthesis necessary for T cell expansion and immunity.

Dominant atopy risk mutations identified by mouse forward genetic analysis.

BACKGROUND: Atopy, the overall tendency to become sensitized to an allergen, is heritable but seldom ascribed to mutations within specific genes. Atopic individuals develop abnormally elevated IgE responses to immunization with potential allergens. To gain insight into the genetic causes of atopy, we carried out a forward genetic screen for atopy in mice. METHODS: We screened mice carrying homozygous and heterozygous N-ethyl-N-nitrosourea (ENU)-induced germline mutations for aberrant antigen-specific IgE and IgG1 production in response to immunization with the model allergen papain. Candidate genes were validated by independent gene mutation. RESULTS: Of 31 candidate genes selected for investigation, the effects of mutations in 23 genes on papain-specific IgE or IgG1 were verified. Among the 20 verified genes influencing the IgE response, eight were necessary for the response, while 12 repressed IgE. Nine genes were not previously implicated in the IgE response. Fifteen genes encoded proteins contributing to IgE class switch recombination or B-cell receptor signaling. The precise roles of the five remaining genes (Flcn, Map1lc3b, Me2, Prkd2, and Scarb2) remain to be determined. Loss-of-function mutations in nine of the 12 genes limiting the IgE response were dominant or semi-dominant for the IgE phenotype but did not cause immunodeficiency in the heterozygous state. Using damaging allele frequencies for the corresponding human genes and in silico simulations (Monte Carlo) of undiscovered atopy mutations, we estimated the percentage of humans with heterozygous atopy risk mutations. CONCLUSIONS: Up to 37% of individuals may be heterozygous carriers for at least one dominant atopy risk mutation.

LocNES: a computational tool for locating classical NESs in CRM1 cargo proteins.

MOTIVATION: Classical nuclear export signals (NESs) are short cognate peptides that direct proteins out of the nucleus via the CRM1-mediated export pathway. CRM1 regulates the localization of hundreds of macromolecules involved in various cellular functions and diseases. Due to the diverse and complex nature of NESs, reliable prediction of the signal remains a challenge despite several attempts made in the last decade. RESULTS: We present a new NES predictor, LocNES. LocNES scans query proteins for NES consensus-fitting peptides and assigns these peptides probability scores using Support Vector Machine model, whose feature set includes amino acid sequence, disorder propensity, and the rank of position-specific scoring matrix score. LocNES demonstrates both higher sensitivity and precision over existing NES prediction tools upon comparative analysis using experimentally identified NESs. AVAILABILITY AND IMPLEMENTATION: LocNES is freely available at http://prodata.swmed.edu/LocNES CONTACT: yuhmin.chook@utsouthwestern.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

NESdb: a database of NES-containing CRM1 cargoes.

The leucine-rich nuclear export signal (NES) is the only known class of targeting signal that directs macromolecules out of the cell nucleus. NESs are short stretches of 8-15 amino acids with regularly spaced hydrophobic residues that bind the export karyopherin CRM1. NES-containing proteins are involved in numerous cellular and disease processes. We compiled a database named NESdb that contains 221 NES-containing CRM1 cargoes that were manually curated from the published literature. Each NESdb entry is annotated with information about sequence and structure of both the NES and the cargo protein, as well as information about experimental evidence of NES-mapping and CRM1-mediated nuclear export. NESdb will be updated regularly and will serve as an important resource for nuclear export signals. NESdb is freely available to nonprofit organizations at http://prodata.swmed.edu/LRNes.

Sequence and structural analyses of nuclear export signals in the NESdb database.

We compiled >200 nuclear export signal (NES)-containing CRM1 cargoes in a database named NESdb. We analyzed the sequences and three-dimensional structures of natural, experimentally identified NESs and of false-positive NESs that were generated from the database in order to identify properties that might distinguish the two groups of sequences. Analyses of amino acid frequencies, sequence logos, and agreement with existing NES consensus sequences revealed strong preferences for the Φ1-X(3)-Φ2-X(2)-Φ3-X-Φ4 pattern and for negatively charged amino acids in the nonhydrophobic positions of experimentally identified NESs but not of false positives. Strong preferences against certain hydrophobic amino acids in the hydrophobic positions were also revealed. These findings led to a new and more precise NES consensus. More important, three-dimensional structures are now available for 68 NESs within 56 different cargo proteins. Analyses of these structures showed that experimentally identified NESs are more likely than the false positives to adopt α-helical conformations that transition to loops at their C-termini and more likely to be surface accessible within their protein domains or be present in disordered or unobserved parts of the structures. Such distinguishing features for real NESs might be useful in future NES prediction efforts. Finally, we also tested CRM1-binding of 40 NESs that were found in the 56 structures. We found that 16 of the NES peptides did not bind CRM1, hence illustrating how NESs are easily misidentified.

Recognition of nuclear targeting signals by Karyopherin-ß proteins.

The Karyopherin-ß family of nuclear transport factors mediates the majority of nucleocytoplasmic transport. Although each of the 19 Karyopherin-ßs transports unique sets of cargos, only three classes of nuclear localization and export signals, or NLSs and NESs, have been characterized. The short basic classical-NLS was first discovered in the 1980s and their karyopherin-bound structures were first reported more than 10 years ago. More recently, structural and biophysical studies of Karyopherin-ß2-cargo complexes led to definition of the complex and diverse PY-NLS. Structural knowledge of the leucine-rich NES is finally available more than 10 years after the discovery of its recognition by the exportin CRM1. We review recent findings relating to how these three classes of nuclear targeting signals are recognized by their Karyopherin-ß nuclear transport factors.

The electrostatic characteristics of G.U wobble base pairs.

G.U wobble base pairs are the most common and highly conserved non-Watson-Crick base pairs in RNA. Previous surface maps imply uniformly negative electrostatic potential at the major groove of G.U wobble base pairs embedded in RNA helices, suitable for entrapment of cationic ligands. In this work, we have used a Poisson-Boltzmann approach to gain a more detailed and accurate characterization of the electrostatic profile. We found that the major groove edge of an isolated G.U wobble displays distinctly enhanced negativity compared with standard GC or AU base pairs; however, in the context of different helical motifs, the electrostatic pattern varies. G.U wobbles with distinct widening have similar major groove electrostatic potentials to their canonical counterparts, whereas those with minimal widening exhibit significantly enhanced electronegativity, ranging from 0.8 to 2.5 kT/e, depending upon structural features. We propose that the negativity at the major groove of G.U wobble base pairs is determined by the combined effect of the base atoms and the sugar-phosphate backbone, which is impacted by stacking pattern and groove width as a result of base sequence. These findings are significant in that they provide predictive power with respect to which G.U sites in RNA are most likely to bind cationic ligands.

Recognition of the spliceosomal branch site RNA helix on the basis of surface and electrostatic features.

We have investigated electrostatic and surface features of an essential region of the catalytic core of the spliceosome, the eukaryotic precursor messenger (pre-m)RNA splicing apparatus. The nucleophile for the first of two splicing reactions is the 2'-hydroxyl (OH) of the ribose of a specific adenosine within the intron. During assembly of the spliceosome's catalytic core, this adenosine is positioned by pairing with a short region of the U2 small nuclear (sn)RNA to form the pre-mRNA branch site helix. The solution structure of the spliceosomal pre-mRNA branch site [Newby,M.I. and Greenbaum,N.L. (2002) Nature Struct. Biol., 9, 958-965] showed that a phylogenetically conserved pseudouridine (psi) residue in the segment of U2 snRNA that pairs with the intron induces a markedly different structure compared with that of its unmodified counterpart. In order to achieve a more detailed understanding of the factors that contribute to recognition of the spliceosome's branch site helix and activation of the nucleophile for the first step of pre-mRNA splicing, we have calculated surface areas and electrostatic potentials of psi-modified and unmodified branch site duplexes. There was no significant difference between the total accessible area or ratio of total polar:nonpolar groups between modified and unmodified duplexes. However, there was substantially greater exposure of nonpolar area of the adenine base, and less exposure of the 2'-OH, in the psi-modified structure. Electrostatic potentials computed using a hybrid boundary element and finite difference nonlinear Poisson-Boltzmann approach [Boschitsch, A.H. and Fenley, M.O. (2004) J. Comput. Chem., 25, 935-955] revealed a region of exceptionally negative potential in the major groove surrounding the 2'-OH of the branch site adenosine. These surface and electrostatic features may contribute to the overall recognition of the pre-mRNA branch site region by other components of the splicing reaction.