Maize HapMap2 identifies extant variation from a genome in flux

Whereas breeders have exploited diversity in maize for yield improvements, there has been limited progress in using beneficial alleles in undomesticated varieties. Characterizing standing variation in this complex genome has been challenging, with only a small fraction of it described to date. Using a population genetics scoring model, we identified 55 million SNPs in 103 lines across pre-domestication and domesticated Zea mays varieties, including a representative from the sister genus Tripsacum. We find that structural variations are pervasive in the Z. mays genome and are enriched at loci associated with important traits. By investigating the drivers of genome size variation, we find that the larger Tripsacum genome can be explained by transposable element abundance rather than an allopolyploid origin. In contrast, intraspecies genome size variation seems to be controlled by chromosomal knob content. There is tremendous overlap in key gene content in maize and Tripsacum, suggesting that adaptations from Tripsacum (for example, perennialism and frost and drought tolerance) can likely be integrated into maize.     

Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm

A high-density haplotype map recently enabled a genome-wide association study (GWAS) in a population of indica subspecies of Chinese rice landraces. Here we extend this methodology to a larger and more diverse sample of 950 worldwide rice varieties, including the Oryza sativa indica and Oryza sativa japonica subspecies, to perform an additional GWAS. We identified a total of 32 new loci associated with flowering time and with ten grain-related traits, indicating that the larger sample increased the power to detect trait-associated variants using GWAS. To characterize various alleles and complex genetic variation, we developed an analytical framework for haplotype-based de novo assembly of the low-coverage sequencing data in rice. We identified candidate genes for 18 associated loci through detailed annotation. This study shows that the integrated approach of sequence-based GWAS and functional genome annotation has the potential to match complex traits to their causal polymorphisms in rice.     

A genome-wide association study in Han Chinese identifies new susceptibility loci for ankylosing spondylitis

To identify susceptibility loci for ankylosing spondylitis, we performed a two-stage genome-wide association study in Han Chinese. In the discovery stage, we analyzed 1,356,350 autosomal SNPs in 1,837 individuals with ankylosing spondylitis and 4,231 controls; in the validation stage, we analyzed 30 suggestive SNPs in an additional 2,100 affected individuals and 3,496 controls. We identified two new susceptibility loci between EDIL3 and HAPLN1 at 5q14.3 (rs4552569; P = 8.77 × 10(-10)) and within ANO6 at 12q12 (rs17095830; P = 1.63 × 10(-8)). We also confirmed previously reported associations in Europeans within the major histocompatibility complex (MHC) region (top SNP, rs13202464; P < 5 × 10(-324)) and at 2p15 (rs10865331; P = 1.98 × 10(-8)). We show that rs13202464 within the MHC region mainly represents the risk effect of HLA-B*27 variants (including HLA-B*2704, HLA-B*2705 and HLA-B*2715) in Chinese. The two newly discovered loci implicate genes related to bone formation and cartilage development, suggesting their potential involvement in the etiology of ankylosing spondylitis.     

A genome-wide association study in Chinese men identifies three risk loci for non-obstructive azoospermia

Non-obstructive azoospermia (NOA) is one of the most severe forms of male infertility. Its pathophysiology is largely unknown, and few genetic influences have been defined. To identify common variants contributing to NOA in Han Chinese men, we performed a three-stage genome-wide association study of 2,927 individuals with NOA and 5,734 controls. The combined analyses identified significant (P < 5.0 × 10(-8)) associations between NOA risk and common variants near PRMT6 (rs12097821 at 1p13.3: odds ratio (OR) = 1.25, P = 5.7 × 10(-10)), PEX10 (rs2477686 at 1p36.32: OR = 1.39, P = 5.7 × 10(-12)) and SOX5 (rs10842262 at 12p12.1: OR = 1.23, P = 2.3 × 10(-9)). These findings implicate genetic variants at 1p13.3, 1p36.32 and 12p12.1 in the etiology of NOA in Han Chinese men.     

Extremely low-coverage sequencing and imputation increases power for genome-wide association studies

Genome-wide association studies (GWAS) have proven to be a powerful method to identify common genetic variants contributing to susceptibility to common diseases. Here, we show that extremely low-coverage sequencing (0.1-0.5×) captures almost as much of the common (>5%) and low-frequency (1-5%) variation across the genome as SNP arrays. As an empirical demonstration, we show that genome-wide SNP genotypes can be inferred at a mean r(2) of 0.71 using off-target data (0.24× average coverage) in a whole-exome study of 909 samples. Using both simulated and real exome-sequencing data sets, we show that association statistics obtained using extremely low-coverage sequencing data attain similar P values at known associated variants as data from genotyping arrays, without an excess of false positives. Within the context of reductions in sample preparation and sequencing costs, funds invested in extremely low-coverage sequencing can yield several times the effective sample size of GWAS based on SNP array data and a commensurate increase in statistical power.     

Mutations in NMNAT1 cause Leber congenital amaurosis and identify a new disease pathway for retinal degeneration

Leber congenital amaurosis (LCA) is a blinding retinal disease that presents within the first year after birth. Using exome sequencing, we identified mutations in the nicotinamide adenine dinucleotide (NAD) synthase gene NMNAT1 encoding nicotinamide mononucleotide adenylyltransferase 1 in eight families with LCA, including the family in which LCA was originally linked to the LCA9 locus. Notably, all individuals with NMNAT1 mutations also have macular colobomas, which are severe degenerative entities of the central retina (fovea) devoid of tissue and photoreceptors. Functional assays of the proteins encoded by the mutant alleles identified in our study showed that the mutations reduce the enzymatic activity of NMNAT1 in NAD biosynthesis and affect protein folding. Of note, recent characterization of the slow Wallerian degeneration (Wld(s)) mouse model, in which prolonged axonal survival after injury is observed, identified NMNAT1 as a neuroprotective protein when ectopically expressed. Our findings identify a new disease mechanism underlying LCA and provide the first link between endogenous NMNAT1 dysfunction and a human nervous system disorder.     

Calcium-induced cleavage of DNA topoisomerase I involves the cytoplasmic-nuclear shuttling of calpain 2

Important to the function of calpains is temporal and spatial regulation of their proteolytic activity. Here, we demonstrate that cytoplasm-resident calpain 2 cleaves human nuclear topoisomerase I (hTOP1) via Ca²⁺-activated proteolysis and nucleoplasmic shuttling of proteases. This proteolysis of hTOP1 was induced by either ionomycin-caused Ca²⁺ influx or addition of Ca²⁺ in cellular extracts. Ca²⁺ failed to induce hTOP1 proteolysis in calpain 2-knockdown cells. Moreover, calpain 2 cleaved hTOP1 in vitro. Furthermore, calpain 2 entered the nucleus upon Ca²⁺ influx, and calpastatin interfered with this process. Calpain 2 cleavage sites were mapped at K¹⁵⁸ and K¹⁸³ of hTOP1. Calpain 2-truncated hTOP1 exhibited greater relaxation activity but remained able to interact with nucleolin and to form cleavable complexes. Interestingly, calpain 2 appears to be involved in ionomycin-induced protection from camptothecin-induced cytotoxicity. Thus, our data suggest that nucleocytoplasmic shuttling may serve as a novel type of regulation for calpain 2-mediated nuclear proteolysis.     

CEP152 is a genome maintenance protein disrupted in Seckel syndrome

Functional impairment of DNA damage response pathways leads to increased genomic instability. Here we describe the centrosomal protein CEP152 as a new regulator of genomic integrity and cellular response to DNA damage. Using homozygosity mapping and exome sequencing, we identified CEP152 mutations in Seckel syndrome and showed that impaired CEP152 function leads to accumulation of genomic defects resulting from replicative stress through enhanced activation of ATM signaling and increased H2AX phosphorylation.     

Genome-wide association and linkage identify modifier loci of lung disease severity in cystic fibrosis at 11p13 and 20q13.2

A combined genome-wide association and linkage study was used to identify loci causing variation in cystic fibrosis lung disease severity. We identified a significant association (P = 3.34 × 10(-8)) near EHF and APIP (chr11p13) in p.Phe508del homozygotes (n = 1,978). The association replicated in p.Phe508del homozygotes (P = 0.006) from a separate family based study (n = 557), with P = 1.49 × 10(-9) for the three-study joint meta-analysis. Linkage analysis of 486 sibling pairs from the family based study identified a significant quantitative trait locus on chromosome 20q13.2 (log(10) odds = 5.03). Our findings provide insight into the causes of variation in lung disease severity in cystic fibrosis and suggest new therapeutic targets for this life-limiting disorder.