Discovery and fine-mapping of height loci via high-density imputation of GWASs in individuals of African ancestry.

Although many loci have been associated with height in European ancestry populations, very few have been identified in African ancestry individuals. Furthermore, many of the known loci have yet to be generalized to and fine-mapped within a large-scale African ancestry sample. We performed sex-combined and sex-stratified meta-analyses in up to 52,764 individuals with height and genome-wide genotyping data from the African Ancestry Anthropometry Genetics Consortium (AAAGC). We additionally combined our African ancestry meta-analysis results with published European genome-wide association study (GWAS) data. In the African ancestry analyses, we identified three novel loci (SLC4A3, NCOA2, ECD/FAM149B1) in sex-combined results and two loci (CRB1, KLF6) in women only. In the African plus European sex-combined GWAS, we identified an additional three novel loci (RCCD1, G6PC3, CEP95) which were equally driven by AAAGC and European results. Among 39 genome-wide significant signals at known loci, conditioning index SNPs from European studies identified 20 secondary signals. Two of the 20 new secondary signals and none of the 8 novel loci had minor allele frequencies (MAF) < 5%. Of 802 known European height signals, 643 displayed directionally consistent associations with height, of which 205 were nominally significant (p < 0.05) in the African ancestry sex-combined sample. Furthermore, 148 of 241 loci contained ≤20 variants in the credible sets that jointly account for 99% of the posterior probability of driving the associations. In summary, trans-ethnic meta-analyses revealed novel signals and further improved fine-mapping of putative causal variants in loci shared between African and European ancestry populations.

Structural basis for substrate binding and catalysis by a self-alkylating ribozyme.

Ribozymes that react with small-molecule probes have important applications in transcriptomics and chemical biology, such as RNA labeling and imaging. Understanding the structural basis for these RNA-modifying reactions will enable the development of better tools for studying RNA. Nevertheless, high-resolution structures and underlying catalytic mechanisms for members of this ribozyme class remain elusive. Here, we focus on a self-alkylating ribozyme that catalyzes nitrogen-carbon bond formation between a specific guanine and a 2,3-disubstituted epoxide substrate and report the crystal structures of a self-alkylating ribozyme, including both alkylated and apo forms, at 1.71-Å and 2.49-Å resolution, respectively. The ribozyme assumes an elongated hairpin-like architecture preorganized to accommodate the epoxide substrate in a hook-shaped conformation. Observed reactivity of substrate analogs together with an inverse, log-linear pH dependence of the reaction rate suggests a requirement for epoxide protonation, possibly assisted by the ether oxygens within the substrate.

Structure of an RNA silencing complex of the CRISPR-Cas immune system.

Bacterial and archaeal clustered regularly interspaced short palindromic repeat (CRISPR) loci capture virus and plasmid sequences and use them to recognize and eliminate these invaders. CRISPR RNAs (crRNAs) containing the acquired sequences are incorporated into effector complexes that destroy matching invader nucleic acids. The multicomponent Cmr effector complex cleaves RNA targets complementary to the crRNAs. Here, we report cryoelectron microscopy reconstruction of a functional Cmr complex bound with a target RNA at ~12 Å. Pairs of the Cmr4 and Cmr5 proteins form a helical core that is asymmetrically capped on each end by distinct pairs of the four remaining subunits: Cmr2 and Cmr3 at the conserved 5' crRNA tag sequence and Cmr1 and Cmr6 near the 3' end of the crRNA. The shape and organization of the RNA-targeting Cmr complex is strikingly similar to the DNA-targeting Cascade complex. Our results reveal a remarkably conserved architecture among very distantly related CRISPR-Cas complexes.

A Crystal Structure of a Functional RNA Molecule Containing an Artificial Nucleobase Pair.

As one of its goals, synthetic biology seeks to increase the number of building blocks in nucleic acids. While efforts towards this goal are well advanced for DNA, they have hardly begun for RNA. Herein, we present a crystal structure for an RNA riboswitch where a stem C:G pair has been replaced by a pair between two components of an artificially expanded genetic-information system (AEGIS), Z and P, (6-amino-5-nitro-2(1H)-pyridone and 2-amino-imidazo[1,2-a]-1,3,5-triazin-4-(8H)-one). The structure shows that the Z:P pair does not greatly change the conformation of the RNA molecule nor the details of its interaction with a hypoxanthine ligand. This was confirmed in solution by in-line probing, which also measured a 3.7 nM affinity of the riboswitch for guanine. These data show that the Z:P pair mimics the natural Watson-Crick geometry in RNA in the first example of a crystal structure of an RNA molecule that contains an orthogonal added nucleobase pair.

Pleiotropic genes for metabolic syndrome and inflammation.

Metabolic syndrome (MetS) has become a health and financial burden worldwide. The MetS definition captures clustering of risk factors that predict higher risk for diabetes mellitus and cardiovascular disease. Our study hypothesis is that additional to genes influencing individual MetS risk factors, genetic variants exist that influence MetS and inflammatory markers forming a predisposing MetS genetic network. To test this hypothesis a staged approach was undertaken. (a) We analyzed 17 metabolic and inflammatory traits in more than 85,500 participants from 14 large epidemiological studies within the Cross Consortia Pleiotropy Group. Individuals classified with MetS (NCEP definition), versus those without, showed on average significantly different levels for most inflammatory markers studied. (b) Paired average correlations between 8 metabolic traits and 9 inflammatory markers from the same studies as above, estimated with two methods, and factor analyses on large simulated data, helped in identifying 8 combinations of traits for follow-up in meta-analyses, out of 130,305 possible combinations between metabolic traits and inflammatory markers studied. (c) We performed correlated meta-analyses for 8 metabolic traits and 6 inflammatory markers by using existing GWAS published genetic summary results, with about 2.5 million SNPs from twelve predominantly largest GWAS consortia. These analyses yielded 130 unique SNPs/genes with pleiotropic associations (a SNP/gene associating at least one metabolic trait and one inflammatory marker). Of them twenty-five variants (seven loci newly reported) are proposed as MetS candidates. They map to genes MACF1, KIAA0754, GCKR, GRB14, COBLL1, LOC646736-IRS1, SLC39A8, NELFE, SKIV2L, STK19, TFAP2B, BAZ1B, BCL7B, TBL2, MLXIPL, LPL, TRIB1, ATXN2, HECTD4, PTPN11, ZNF664, PDXDC1, FTO, MC4R and TOMM40. Based on large data evidence, we conclude that inflammation is a feature of MetS and several gene variants show pleiotropic genetic associations across phenotypes and might explain a part of MetS correlated genetic architecture. These findings warrant further functional investigation.

Structures of artificially designed discrete RNA nanoarchitectures at near-atomic resolution.

Although advances in nanotechnology have enabled the construction of complex and functional synthetic nucleic acid­based nanoarchitectures, high-resolution discrete structures are lacking because of the difficulty in obtaining good diffracting crystals. Here, we report the design and construction of RNA nanostructures based on homooligomerizable one-stranded tiles for x-ray crystallographic determination. We solved three structures to near-atomic resolution: a 2D parallelogram, a 3D nanobracelet unexpectedly formed from an RNA designed for a nanocage, and, eventually, a bona fide 3D nanocage designed with the guidance of the two previous structures. Structural details of their constituent motifs, such as kissing loops, branched kissing loops, and T-junctions, that resemble natural RNA motifs and resisted x-ray determination are revealed, providing insights into those natural motifs. This work unveils the largely unexplored potential of crystallography in gaining high-resolution feedback for nanoarchitectural design and suggests a route to investigate RNA motif structures by configuring them into nanoarchitectures.

Branched kissing loops for the construction of diverse RNA homooligomeric nanostructures.

In biological systems, large and complex structures are often assembled from multiple simpler identical subunits. This strategy-homooligomerization-allows efficient genetic encoding of structures and avoids the need to control the stoichiometry of multiple distinct units. It also allows the minimal number of distinct subunits when designing artificial nucleic acid structures. Here, we present a robust self-assembly system in which homooligomerizable tiles are formed from intramolecularly folded RNA single strands. Tiles are linked through an artificially designed branched kissing-loop motif, involving Watson-Crick base pairing between the single-stranded regions of a bulged helix and a hairpin loop. By adjusting the tile geometry to gain control over the curvature, torsion and the number of helices, we have constructed 16 different linear and circular structures, including a finite-sized three-dimensional cage. We further demonstrate cotranscriptional self-assembly of tiles based on branched kissing loops, and show that tiles inserted into a transfer RNA scaffold can be overexpressed in bacterial cells.

A conserved RNA structural motif for organizing topology within picornaviral internal ribosome entry sites.

Picornaviral IRES elements are essential for initiating the cap-independent viral translation. However, three-dimensional structures of these elements remain elusive. Here, we report a 2.84-Å resolution crystal structure of hepatitis A virus IRES domain V (dV) in complex with a synthetic antibody fragment-a crystallization chaperone. The RNA adopts a three-way junction structure, topologically organized by an adenine-rich stem-loop motif. Despite no obvious sequence homology, the dV architecture shows a striking similarity to a circularly permuted form of encephalomyocarditis virus J-K domain, suggesting a conserved strategy for organizing the domain architecture. Recurrence of the motif led us to use homology modeling tools to compute a 3-dimensional structure of the corresponding domain of foot-and-mouth disease virus, revealing an analogous domain organizing motif. The topological conservation observed among these IRESs and other viral domains implicates a structured three-way junction as an architectural scaffold to pre-organize helical domains for recruiting the translation initiation machinery.

Exome Chip Meta-analysis Fine Maps Causal Variants and Elucidates the Genetic Architecture of Rare Coding Variants in Smoking and Alcohol Use.

BACKGROUND: Smoking and alcohol use have been associated with common genetic variants in multiple loci. Rare variants within these loci hold promise in the identification of biological mechanisms in substance use. Exome arrays and genotype imputation can now efficiently genotype rare nonsynonymous and loss of function variants. Such variants are expected to have deleterious functional consequences and to contribute to disease risk. METHODS: We analyzed ∼250,000 rare variants from 16 independent studies genotyped with exome arrays and augmented this dataset with imputed data from the UK Biobank. Associations were tested for five phenotypes: cigarettes per day, pack-years, smoking initiation, age of smoking initiation, and alcoholic drinks per week. We conducted stratified heritability analyses, single-variant tests, and gene-based burden tests of nonsynonymous/loss-of-function coding variants. We performed a novel fine-mapping analysis to winnow the number of putative causal variants within associated loci. RESULTS: Meta-analytic sample sizes ranged from 152,348 to 433,216, depending on the phenotype. Rare coding variation explained 1.1% to 2.2% of phenotypic variance, reflecting 11% to 18% of the total single nucleotide polymorphism heritability of these phenotypes. We identified 171 genome-wide associated loci across all phenotypes. Fine mapping identified putative causal variants with double base-pair resolution at 24 of these loci, and between three and 10 variants for 65 loci. Twenty loci contained rare coding variants in the 95% credible intervals. CONCLUSIONS: Rare coding variation significantly contributes to the heritability of smoking and alcohol use. Fine-mapping genome-wide association study loci identifies specific variants contributing to the biological etiology of substance use behavior.

Structural basis for activation of fluorogenic dyes by an RNA aptamer lacking a G-quadruplex motif.

The DIR2s RNA aptamer, a second-generation, in-vitro selected binder to dimethylindole red (DIR), activates the fluorescence of cyanine dyes, DIR and oxazole thiazole blue (OTB), allowing detection of two well-resolved emission colors. Using Fab BL3-6 and its cognate hairpin as a crystallization module, we solved the crystal structures of both the apo and OTB-SO3 bound forms of DIR2s at 2.0 Å and 1.8 Å resolution, respectively. DIR2s adopts a compact, tuning fork-like architecture comprised of a helix and two short stem-loops oriented in parallel to create the ligand binding site through tertiary interactions. The OTB-SO3 fluorophore binds in a planar conformation to a claw-like structure formed by a purine base-triple, which provides a stacking platform for OTB-SO3, and an unpaired nucleotide, which partially caps the binding site from the top. The absence of a G-quartet or base tetrad makes the DIR2s aptamer unique among fluorogenic RNAs with known 3D structure.

The nicotinic acetylcholine receptor gene family of the silkworm, Bombyx mori.

BACKGROUND: Nicotinic acetylcholine receptors (nAChRs) mediate fast synaptic cholinergic transmission in the insect central nervous system. The insect nAChR is the molecular target of a class of insecticides, neonicotinoids. Like mammalian nAChRs, insect nAChRs are considered to be made up of five subunits, coded by homologous genes belonging to the same family. The nAChR subunit genes of Drosophila melanogaster, Apis mellifera and Anopheles gambiae have been cloned previously based on their genome sequences. The silkworm Bombyx mori is a model insect of Lepidoptera, among which are many agricultural pests. Identification and characterization of B. mori nAChR genes could provide valuable basic information for this important family of receptor genes and for the study of the molecular mechanisms of neonicotinoid action and resistance. RESULTS: We searched the genome sequence database of B. mori with the fruit fly and honeybee nAChRs by tBlastn and cloned all putative silkworm nAChR cDNAs by reverse transcriptase-polymerase chain reaction (RT-PCR) and rapid amplification of cDNA ends (RACE) methods. B. mori appears to have the largest known insect nAChR gene family to date, including nine alpha-type subunits and three beta-type subunits. The silkworm possesses three genes having low identity with others, including one alpha and two beta subunits, alpha 9, beta2 and beta 3. Like the fruit fly and honeybee counterparts, silkworm nAChR gene alpha 6 has RNA-editing sites, and alpha 4, alpha 6 and alpha 8 undergo alternative splicing. In particular, alternative exon 7 of Bm alpha 8 may have arisen from a recent duplication event. Truncated transcripts were found for Bm alpha 4 and Bm alpha 5. CONCLUSION: B. mori possesses a largest known insect nAChR gene family characterized to date, including nine alpha-type subunits and three beta-type subunits. RNA-editing, alternative splicing and truncated transcripts were found in several subunit genes, which might enhance the diversity of the gene family.

Synthesizing topological structures containing RNA.

Though knotting and entanglement have been observed in DNA and proteins, their existence in RNA remains an enigma. Synthetic RNA topological structures are significant for understanding the physical and biological properties pertaining to RNA topology, and these properties in turn could facilitate identifying naturally occurring topologically nontrivial RNA molecules. Here we show that topological structures containing single-stranded RNA (ssRNA) free of strong base pairing interactions can be created either by configuring RNA-DNA hybrid four-way junctions or by template-directed synthesis with a single-stranded DNA (ssDNA) topological structure. By using a constructed ssRNA knot as a highly sensitive topological probe, we find that Escherichia coli DNA topoisomerase I has low RNA topoisomerase activity and that the R173A point mutation abolishes the unknotting activity for ssRNA, but not for ssDNA. Furthermore, we discover the topological inhibition of reverse transcription (RT) and obtain different RT-PCR patterns for an ssRNA knot and circle of the same sequence.

A Non-Stem-Loop CRISPR RNA Is Processed by Dual Binding Cas6.

A subclass of recently discovered CRISPR repeat RNA in bacteria contains minimally recognizable structural features that facilitate an unknown mechanism of recognition and processing by the Cas6 family of endoribonucleases. Cocrystal structures of Cas6 from Methanococcus maripaludis (MmCas6b) bound with its repeat RNA revealed a dual site binding structure and a cleavage site conformation poised for phosphodiester bond breakage. Two non-interacting MmCas6b bind to two separate AAYAA motifs within the same repeat, one distal and one adjacent to the cleavage site. This bound structure potentially competes with a stable but non-productive RNA structure. At the cleavage site, MmCas6b supplies a base pair mimic to stabilize a short 2 base pair stem immediately upstream of the scissile phosphate. Complementary biochemical analyses support the dual-AAYAA binding model and a critical role of the protein-RNA base pair mimic. Our results reveal a previously unknown method of processing non-stem-loop CRISPR RNA by Cas6.

Specific Recognition of a Single-Stranded RNA Sequence by a Synthetic Antibody Fragment.

Antibodies that bind RNA represent an unrealized source of reagents for synthetic biology and for characterizing cellular transcriptomes. However, facile access to RNA-binding antibodies requires the engineering of effective Fab libraries guided by the knowledge of the principles that govern RNA recognition. Here, we describe a Fab identified from a minimalist synthetic library during phage display against a branched RNA target. The Fab (BRG) binds with 20nM dissociation constant to a single-stranded RNA (ssRNA) sequence adjacent to the branch site and can block the action of debranchase enzyme. We report the crystal structure in complex with RNA target at 2.38Å. The Fab traps the RNA in a hairpin conformation that contains a 2-bp duplex capped by a tetraloop. The paratope surface consists of residues located in four complementarity-determining regions including a major contribution from H3, which adopts a helical structure that projects into a deep, wide groove formed by the RNA. The amino acid composition of the paratope reflects the library diversity, consisting mostly of tyrosine and serine residues and a small but significant contribution from a single arginine residue. This structure, involving the recognition of ssRNA via a stem-loop conformation, together with our two previous structures involving the recognition of an RNA hairpin loop and an RNA tertiary structure, reveals the capacity of minimalist libraries biased with tyrosine, serine, glycine, and arginine to form binding surfaces for specific RNA conformations and distinct levels of RNA structural hierarchy.

Inhibition of cGAS DNA Sensing by a Herpesvirus Virion Protein.

Invading viral DNA can be recognized by the host cytosolic DNA sensor, cyclic GMP-AMP (cGAMP) synthase (cGAS), resulting in production of the second messenger cGAMP, which directs the adaptor protein STING to stimulate production of type I interferons (IFNs). Although several DNA viruses are sensed by cGAS, viral strategies targeting cGAS are virtually unknown. We report here that Kaposi's sarcoma-associated herpesvirus (KSHV) ORF52, an abundant gammaherpesvirus-specific tegument protein, subverts cytosolic DNA sensing by directly inhibiting cGAS enzymatic activity through a mechanism involving both cGAS binding and DNA binding. Moreover, ORF52 homologs in other gammaherpesviruses also inhibit cGAS activity and similarly bind cGAS and DNA, suggesting conserved inhibitory mechanisms. Furthermore, KSHV infection evokes cGAS-dependent responses that can limit the infection, and an ORF52 null mutant exhibits increased cGAS signaling. Our findings reveal a mechanism through which gammaherpesviruses antagonize host cGAS DNA sensing.

Essential structural and functional roles of the Cmr4 subunit in RNA cleavage by the Cmr CRISPR-Cas complex.

The Cmr complex is the multisubunit effector complex of the type III-B clustered regularly interspaced short palindromic repeats (CRISPR)-Cas immune system. The Cmr complex recognizes a target RNA through base pairing with the integral CRISPR RNA (crRNA) and cleaves the target at multiple regularly spaced locations within the complementary region. To understand the molecular basis of the function of this complex, we have assembled information from electron microscopic and X-ray crystallographic structural studies and mutagenesis of a complete Pyrococcus furiosus Cmr complex. Our findings reveal that four helically packed Cmr4 subunits, which make up the backbone of the Cmr complex, act as a platform to support crRNA binding and target RNA cleavage. Interestingly, we found a hook-like structural feature associated with Cmr4 that is likely the site of target RNA binding and cleavage. Our results also elucidate analogies in the mechanisms of crRNA and target molecule binding by the distinct Cmr type III-A and Cascade type I-E complexes.

Sequence variation in TMEM18 in association with body mass index: Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium Targeted Sequencing Study.

BACKGROUND: Genome-wide association studies for body mass index (BMI) previously identified a locus near TMEM18. We conducted targeted sequencing of this region to investigate the role of common, low-frequency, and rare variants influencing BMI. METHODS AND RESULTS: We sequenced TMEM18 and regions downstream of TMEM18 on chromosome 2 in 3976 individuals of European ancestry from 3 community-based cohorts (Atherosclerosis Risk in Communities, Cardiovascular Health Study, and Framingham Heart Study), including 200 adults selected for high BMI. We examined the association between BMI and variants identified in the region from nucleotide position 586 432 to 677 539 (hg18). Rare variants (minor allele frequency, <1%) were analyzed using a burden test and the sequence kernel association test. Results from the 3 cohort studies were meta-analyzed. We estimate that mean BMI is 0.43 kg/m(2) higher for each copy of the G allele of single-nucleotide polymorphism rs7596758 (minor allele frequency, 29%; P=3.46×10(-4)) using a Bonferroni threshold of P<4.6×10(-4). Analyses conditional on previous genome-wide association study single-nucleotide polymorphisms associated with BMI in the region led to attenuation of this signal and uncovered another independent (r(2)<0.2), statistically significant association, rs186019316 (P=2.11×10(-4)). Both rs186019316 and rs7596758 or proxies are located in transcription factor binding regions. No significant association with rare variants was found in either the exons of TMEM18 or the 3' genome-wide association study region. CONCLUSIONS: Targeted sequencing around TMEM18 identified 2 novel BMI variants with possible regulatory function.

Recognition and cleavage of a nonstructured CRISPR RNA by its processing endoribonuclease Cas6.

Clustered regularly interspaced short palindromic repeats (CRISPRs) confer adaptive immunity to prokaryotes through a small RNA-mediated mechanism. Specific endoribonucleases are required by all CRISPR-bearing organisms to process CRISPR RNAs into small RNA that serve as guides for defensive effector complexes. The molecular mechanism of how the endoribonucleases process the class of CRISPR RNA containing no predicted secondary structural features remains largely elusive. Here, we report cocrystal structures of a processing endoribonuclease bound with a noncleavable RNA substrate and its product-like fragment derived from a nonpalindramic repeat. The enzyme stabilizes a short RNA stem-loop structure near the cleavage site and cleaves the phosphodiester bond using an active site comprised of arginine and lysine residues. The distinct RNA binding and cleavage mechanisms underline the diversity in CRISPR RNA processing.

Structure of the Cmr2-Cmr3 subcomplex of the Cmr RNA silencing complex.

The Cmr complex is an RNA-guided effector complex that cleaves invader RNA in the prokaryotic immune response mediated by the CRISPR (Clustered Regularly Interspaced Short Palindromic Repeat)-Cas system. Here, we report the crystal structure of a Cmr subcomplex containing Cmr2 (Cas10) and Cmr3 subunits at 2.8 Å resolution. The structure revealed a dual ferredoxin fold and glycine-rich loops characteristic of previously known repeat-associated mysterious proteins and two unique insertion elements in Cmr3 that mediate its interaction with Cmr2. Surprisingly, while mutation of both insertion elements significantly weakened Cmr3-Cmr2 interaction, they exhibit differential effects on Cmr-mediated RNA cleavage by the Cmr complex, suggesting stabilization of Cmr2-Cmr3 interactions by other subunits. Further mutational analysis of the two conserved (but non-Cmr2-binding) glycine-rich loops of Cmr3 identified a region that is likely involved in assembly or the RNA cleavage function of the Cmr complex.