Genome-wide genetic changes during modern breeding of maize

The success of modern maize breeding has been demonstrated by remarkable increases in productivity over the last four decades. However, the underlying genetic changes correlated with these gains remain largely unknown. We report here the sequencing of 278 temperate maize inbred lines from different stages of breeding history, including deep resequencing of 4 lines with known pedigree information. The results show that modern breeding has introduced highly dynamic genetic changes into the maize genome. Artificial selection has affected thousands of targets, including genes and non-genic regions, leading to a reduction in nucleotide diversity and an increase in the proportion of rare alleles. Genetic changes during breeding happen rapidly, with extensive variation (SNPs, indels and copy-number variants (CNVs)) occurring, even within identity-by-descent regions. Our genome-wide assessment of genetic changes during modern maize breeding provides new strategies as well as practical targets for future crop breeding and biotechnology.     

Automating resequencing-based detection of insertion-deletion polymorphisms

Structural and insertion-deletion (indel) variants have received considerable recent attention, partly because of their phenotypic consequences. Among these variants, the most common are small indels ( approximately 1-30 bp). Identifying and genotyping indels using sequence traces obtained from diploid samples requires extensive manual review, which makes large-scale studies inconvenient. We report a new algorithm, implemented in available software (PolyPhred version 6.0), to help automate detection and genotyping of indels from sequence traces. The algorithm identifies heterozygous individuals, which permits the discovery of low-frequency indels. It finds 80% of all indel polymorphisms with almost no false positives and finds 97% with a false discovery rate of 10%. Additionally, genotyping accuracy exceeds 99%, and it correctly infers indel length in 96% of the cases. Using this approach, we identify indels in the HapMap ENCODE regions, providing the first report of these polymorphisms in this data set.     

Estimating coverage and power for genetic association studies using near-complete variation data

Although studies suggest that SNPs derived from HapMap provide promising coverage and power for association studies, the lack of alternative variation datasets limits independent analysis. Using near-complete variation data for 76 genes resequenced in HapMap samples, we find that coverage of common variation by commercial genotyping arrays is substantially lower compared to the HapMap-based estimates. We quantify the power offered by these arrays for a range of disease models.     

Dwarf8 polymorphisms associate with variation in flowering time.

Historically, association tests have been used extensively in medical genetics, but have had virtually no application in plant genetics. One obstacle to their application is the structured populations often found in crop plants, which may lead to nonfunctional, spurious associations. In this study, statistical methods to account for population structure were extended for use with quantitative variation and applied to our evaluation of maize flowering time. Mutagenesis and quantitative trait locus (QTL) studies suggested that the maize gene Dwarf8 might affect the quantitative variation of maize flowering time and plant height. The wheat orthologs of this gene contributed to the increased yields seen in the 'Green Revolution' varieties. We used association approaches to evaluate Dwarf8 sequence polymorphisms from 92 maize inbred lines. Population structure was estimated using a Bayesian analysis of 141 simple sequence repeat (SSR) loci. Our results indicate that a suite of polymorphisms associate with differences in flowering time, which include a deletion that may alter a key domain in the coding region. The distribution of nonsynonymous polymorphisms suggests that Dwarf8 has been a target of selection.     

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 quantitative resistance to southern leaf blight in the maize nested association mapping population

Nested association mapping (NAM) offers power to resolve complex, quantitative traits to their causal loci. The maize NAM population, consisting of 5,000 recombinant inbred lines (RILs) from 25 families representing the global diversity of maize, was evaluated for resistance to southern leaf blight (SLB) disease. Joint-linkage analysis identified 32 quantitative trait loci (QTLs) with predominantly small, additive effects on SLB resistance. Genome-wide association tests of maize HapMap SNPs were conducted by imputing founder SNP genotypes onto the NAM RILs. SNPs both within and outside of QTL intervals were associated with variation for SLB resistance. Many of these SNPs were within or near sequences homologous to genes previously shown to be involved in plant disease resistance. Limited linkage disequilibrium was observed around some SNPs associated with SLB resistance, indicating that the maize NAM population enables high-resolution mapping of some genome regions.     

PNPLA1 mutations cause autosomal recessive congenital ichthyosis in golden retriever dogs and humans

Ichthyoses comprise a heterogeneous group of genodermatoses characterized by abnormal desquamation over the whole body, for which the genetic causes of several human forms remain unknown. We used a spontaneous dog model in the golden retriever breed, which is affected by a lamellar ichthyosis resembling human autosomal recessive congenital ichthyoses (ARCI), to carry out a genome-wide association study. We identified a homozygous insertion-deletion (indel) mutation in PNPLA1 that leads to a premature stop codon in all affected golden retriever dogs. We subsequently found one missense and one nonsense mutation in the catalytic domain of human PNPLA1 in six individuals with ARCI from two families. Further experiments highlighted the importance of PNPLA1 in the formation of the epidermal lipid barrier. This study identifies a new gene involved in human ichthyoses and provides insights into the localization and function of this yet uncharacterized member of the PNPLA protein family.     

A unified mixed-model method for association mapping that accounts for multiple levels of relatedness

As population structure can result in spurious associations, it has constrained the use of association studies in human and plant genetics. Association mapping, however, holds great promise if true signals of functional association can be separated from the vast number of false signals generated by population structure. We have developed a unified mixed-model approach to account for multiple levels of relatedness simultaneously as detected by random genetic markers. We applied this new approach to two samples: a family-based sample of 14 human families, for quantitative gene expression dissection, and a sample of 277 diverse maize inbred lines with complex familial relationships and population structure, for quantitative trait dissection. Our method demonstrates improved control of both type I and type II error rates over other methods. As this new method crosses the boundary between family-based and structured association samples, it provides a powerful complement to currently available methods for association mapping.     

Mutations in LGI1 cause autosomal-dominant partial epilepsy with auditory features.

The epilepsies are a common, clinically heterogeneous group of disorders defined by recurrent unprovoked seizures. Here we describe identification of the causative gene in autosomal-dominant partial epilepsy with auditory features (ADPEAF, MIM 600512), a rare form of idiopathic lateral temporal lobe epilepsy characterized by partial seizures with auditory disturbances. We constructed a complete, 4.2-Mb physical map across the genetically implicated disease-gene region, identified 28 putative genes (Fig. 1) and resequenced all or part of 21 genes before identifying presumptive mutations in one copy of the leucine-rich, glioma-inactivated 1 gene (LGI1) in each of five families with ADPEAF. Previous studies have indicated that loss of both copies of LGI1 promotes glial tumor progression. We show that the expression pattern of mouse Lgi1 is predominantly neuronal and is consistent with the anatomic regions involved in temporal lobe epilepsy. Discovery of LGI1 as a cause of ADPEAF suggests new avenues for research on pathogenic mechanisms of idiopathic epilepsies.     

Quality and completeness of SNP databases

To address the quality and completeness of single-nucleotide polymorphism (SNP) databases, we resequenced 173 kb (spanning 17 loci) in 150 chromosomes of west African and European ancestry. Over 88% of SNPs in the public (TSC and BAC overlap) and Celera databases were confirmed in independent resequencing. Approximately 45% of all human heterozygosity is attributable to SNPs already available from the two databases, and of SNPs with minor-allele frequencies >10%, more than half are represented.     

The heterochronic maize mutant Corngrass1 results from overexpression of a tandem microRNA

Retention of juvenile traits in the adult reproductive phase characterizes a process known as neoteny, and speculation exists over whether it has contributed to the evolution of new species. The dominant Corngrass1 (Cg1) mutant of maize is a neotenic mutation that results in phenotypes that may be present in the grass-like ancestors of maize. We cloned Cg1 and found that it encodes two tandem miR156 genes that are overexpressed in the meristem and lateral organs. Furthermore, a target of Cg1 is teosinte glume architecture1 (tga1), a gene known to have had a role in the domestication of maize from teosinte. Cg1 mutant plants overexpressing miR156 have lower levels of mir172, a microRNA that targets genes controlling juvenile development. By altering the relative levels of both microRNAs, it is possible to either prolong or shorten juvenile development in maize, thus providing a mechanism for how species-level heterochronic changes can occur in nature.     

Additional SNPs and linkage-disequilibrium analyses are necessary for whole-genome association studies in humans

More than 5 million single-nucleotide polymorphisms (SNPs) with minor-allele frequency greater than 10% are expected to exist in the human genome. Some of these SNPs may be associated with risk of developing common diseases. To assess the power of currently available SNPs to detect such associations, we resequenced 50 genes in two ethnic samples and measured patterns of linkage disequilibrium between the subset of SNPs reported in dbSNP and the complete set of common SNPs. Our results suggest that using all 2.7 million SNPs currently in the database would detect nearly 80% of all common SNPs in European populations but only 50% of those common in the African American population and that efficient selection of a minimal subset of SNPs for use in association studies requires measurement of allele frequency and linkage disequilibrium relationships for all SNPs in dbSNP.     

Mutations in the chromatin modifier gene KANSL1 cause the 17q21.31 microdeletion syndrome

We show that haploinsufficiency of KANSL1 is sufficient to cause the 17q21.31 microdeletion syndrome, a multisystem disorder characterized by intellectual disability, hypotonia and distinctive facial features. The KANSL1 protein is an evolutionarily conserved regulator of the chromatin modifier KAT8, which influences gene expression through histone H4 lysine 16 (H4K16) acetylation. RNA sequencing studies in cell lines derived from affected individuals and the presence of learning deficits in Drosophila melanogaster mutants suggest a role for KANSL1 in neuronal processes.     

MicroRNAs can generate thresholds in target gene expression

MicroRNAs (miRNAs) are short, highly conserved noncoding RNA molecules that repress gene expression in a sequence-dependent manner. We performed single-cell measurements using quantitative fluorescence microscopy and flow cytometry to monitor a target gene's protein expression in the presence and absence of regulation by miRNA. We find that although the average level of repression is modest, in agreement with previous population-based measurements, the repression among individual cells varies dramatically. In particular, we show that regulation by miRNAs establishes a threshold level of target mRNA below which protein production is highly repressed. Near this threshold, protein expression responds sensitively to target mRNA input, consistent with a mathematical model of molecular titration. These results show that miRNAs can act both as a switch and as a fine-tuner of gene expression.     

Human cells are protected from mitochondrial dysfunction by complementation of DNA products in fused mitochondria.

Extensive complementation between fused mitochondria is indicated by recombination of 'parental' mitochondrial (mt) DNA (ref. 1,2) of yeast and plant cells. It has been difficult, however, to demonstrate the occurrence of complementation between fused mitochondria in mammalian species through the presence of recombinant mtDNA molecules, because sequence of mtDNA throughout an individual tends to be uniform owing to its strictly maternal inheritance. We isolated two types of respiration-deficient cell lines, with pathogenic mutations in mitochondrial tRNAIle or tRNALeu(UUR) genes from patients with mitochondrial diseases. The coexistence of their mitochondria within hybrids restored their normal morphology and respiratory enzyme activity by 10-14 days after fusion, indicating the presence of an extensive and continuous exchange of genetic contents between the mitochondria. This complementation between fused mitochondria may represent a defence of highly oxidative organelles against mitochondrial dysfunction caused by the accumulation of mtDNA lesions with age.