Genome-wide analysis of deletions in maize population reveals abundant genetic diversity and functional impact.

KEY MESSAGE: Two read depth methods were jointly used in next-generation sequencing data to identify deletions in maize population. GWAS by deletions were analyzed for gene expression pattern and classical traits, respectively. Many studies have confirmed that structural variation (SV) is pervasive throughout the maize genome. Deletion is one type of SV that may impact gene expression and cause phenotypic changes in quantitative traits. In this study, two read count approaches were used to analyze the deletions in the whole-genome sequencing data of 270 maize inbred lines. A total of 19,754 deletion windows overlapped 12,751 genes, which were unevenly distributed across the genome. The deletions explained population structure well and correlated with genomic features. The deletion proportion of genes was determined to be negatively correlated with its expression. The detection of gene expression quantitative trait loci (eQTL) indicated that local eQTL were fewer but had larger effects than distant ones. The common associated genes were related to basic metabolic processes, whereas unique associated genes with eQTL played a role in the stress or stimulus responses in multiple tissues. Compared with the eQTL detected by SNPs derived from the same sequencing data, 89.4% of the associated genes could be detected by both markers. The effect of top eQTL detected by SNPs was usually larger than that detected by deletions for the same gene. A genome-wide association study (GWAS) on flowering time and plant height illustrated that only a few loci could be consistently captured by SNPs, suggesting that combining deletion and SNP for GWAS was an excellent strategy to dissect trait architecture. Our findings will provide insights into characteristic and biological function of genome-wide deletions in maize.

Coregulation of ribosomal RNA with hundreds of genes contributes to phenotypic variation.

Ribosomal repeats occupy 5% of a plant genome, yet there has been little study of their diversity in the modern age of genomics. Ribosomal copy number and expression variation present an opportunity to tap a novel source of diversity. In the present study, we estimated the ribosomal DNA (rDNA) copy number and ribosomal RNA (rRNA) expression for a population of maize inbred lines and investigated the potential role of rDNA and rRNA dosage in regulating global gene expression. Extensive variation was found in both ribosomal DNA copy number and ribosomal RNA expression among maize inbred lines. However, rRNA abundance was not consistent with the copy number of the rDNA. We have not found that the rDNA gene dosage has a regulatory role in gene expression; however, thousands of genes are identified to be coregulated with rRNA expression, including genes participating in ribosome biogenesis and other functionally relevant pathways. We further investigated the potential roles of copy number and the expression level of rDNA on agronomic traits and found that both correlated with flowering time but through different regulatory mechanisms. This comprehensive analysis suggested that rRNA expression variation is a valuable source of functional diversity that affects gene expression variation and field-based phenotypic changes.

The maize W22 genome provides a foundation for functional genomics and transposon biology.

The maize W22 inbred has served as a platform for maize genetics since the mid twentieth century. To streamline maize genome analyses, we have sequenced and de novo assembled a W22 reference genome using short-read sequencing technologies. We show that significant structural heterogeneity exists in comparison to the B73 reference genome at multiple scales, from transposon composition and copy number variation to single-nucleotide polymorphisms. The generation of this reference genome enables accurate placement of thousands of Mutator (Mu) and Dissociation (Ds) transposable element insertions for reverse and forward genetics studies. Annotation of the genome has been achieved using RNA-seq analysis, differential nuclease sensitivity profiling and bisulfite sequencing to map open reading frames, open chromatin sites and DNA methylation profiles, respectively. Collectively, the resources developed here integrate W22 as a community reference genome for functional genomics and provide a foundation for the maize pan-genome.

Construction of the third-generation Zea mays haplotype map.

Background: Characterization of genetic variations in maize has been challenging, mainly due to deterioration of collinearity between individual genomes in the species. An international consortium of maize research groups combined resources to develop the maize haplotype version 3 (HapMap 3), built from whole-genome sequencing data from 1218 maize lines, covering predomestication and domesticated Zea mays varieties across the world. Results: A new computational pipeline was set up to process more than 12 trillion bp of sequencing data, and a set of population genetics filters was applied to identify more than 83 million variant sites. Conclusions: We identified polymorphisms in regions where collinearity is largely preserved in the maize species. However, the fact that the B73 genome used as the reference only represents a fraction of all haplotypes is still an important limiting factor.

Genomic estimation of complex traits reveals ancient maize adaptation to temperate North America.

By 4000 years ago, people had introduced maize to the southwestern United States; full agriculture was established quickly in the lowland deserts but delayed in the temperate highlands for 2000 years. We test if the earliest upland maize was adapted for early flowering, a characteristic of modern temperate maize. We sequenced fifteen 1900-year-old maize cobs from Turkey Pen Shelter in the temperate Southwest. Indirectly validated genomic models predicted that Turkey Pen maize was marginally adapted with respect to flowering, as well as short, tillering, and segregating for yellow kernel color. Temperate adaptation drove modern population differentiation and was selected in situ from ancient standing variation. Validated prediction of polygenic traits improves our understanding of ancient phenotypes and the dynamics of environmental adaptation.

Host genetic variation impacts microbiome composition across human body sites.

BACKGROUND: The composition of bacteria in and on the human body varies widely across human individuals, and has been associated with multiple health conditions. While microbial communities are influenced by environmental factors, some degree of genetic influence of the host on the microbiome is also expected. This study is part of an expanding effort to comprehensively profile the interactions between human genetic variation and the composition of this microbial ecosystem on a genome- and microbiome-wide scale. RESULTS: Here, we jointly analyze the composition of the human microbiome and host genetic variation. By mining the shotgun metagenomic data from the Human Microbiome Project for host DNA reads, we gathered information on host genetic variation for 93 individuals for whom bacterial abundance data are also available. Using this dataset, we identify significant associations between host genetic variation and microbiome composition in 10 of the 15 body sites tested. These associations are driven by host genetic variation in immunity-related pathways, and are especially enriched in host genes that have been previously associated with microbiome-related complex diseases, such as inflammatory bowel disease and obesity-related disorders. Lastly, we show that host genomic regions associated with the microbiome have high levels of genetic differentiation among human populations, possibly indicating host genomic adaptation to environment-specific microbiomes. CONCLUSIONS: Our results highlight the role of host genetic variation in shaping the composition of the human microbiome, and provide a starting point toward understanding the complex interaction between human genetics and the microbiome in the context of human evolution and disease.

Ab initio water pair potential with flexible monomers.

A potential energy surface for the water dimer with explicit dependence on monomer coordinates is presented. The surface was fitted to a set of previously published interaction energies computed on a grid of over a quarter million points in the 12-dimensional configurational space using symmetry-adapted perturbation theory and coupled-cluster methods. The present fit removes small errors in published fits, and its accuracy is critically evaluated. The minimum and saddle-point structures of the potential surface were found to be very close to predictions from direct ab initio optimizations. The computed second virial coefficients agreed well with experimental values. At low temperatures, the effects of monomer flexibility in the virial coefficients were found to be much smaller than the quantum effects.

The role of cis regulatory evolution in maize domestication.

Gene expression differences between divergent lineages caused by modification of cis regulatory elements are thought to be important in evolution. We assayed genome-wide cis and trans regulatory differences between maize and its wild progenitor, teosinte, using deep RNA sequencing in F1 hybrid and parent inbred lines for three tissue types (ear, leaf and stem). Pervasive regulatory variation was observed with approximately 70% of ∼17,000 genes showing evidence of regulatory divergence between maize and teosinte. However, many fewer genes (1,079 genes) show consistent cis differences with all sampled maize and teosinte lines. For ∼70% of these 1,079 genes, the cis differences are specific to a single tissue. The number of genes with cis regulatory differences is greatest for ear tissue, which underwent a drastic transformation in form during domestication. As expected from the domestication bottleneck, maize possesses less cis regulatory variation than teosinte with this deficit greatest for genes showing maize-teosinte cis regulatory divergence, suggesting selection on cis regulatory differences during domestication. Consistent with selection on cis regulatory elements, genes with cis effects correlated strongly with genes under positive selection during maize domestication and improvement, while genes with trans regulatory effects did not. We observed a directional bias such that genes with cis differences showed higher expression of the maize allele more often than the teosinte allele, suggesting domestication favored up-regulation of gene expression. Finally, this work documents the cis and trans regulatory changes between maize and teosinte in over 17,000 genes for three tissues.

Defining the Role of prolamin-box binding factor1 Gene During Maize Domestication.

The prolamin-box binding factor1 (pbf1) gene encodes a transcription factor that controls the expression of seed storage protein (zein) genes in maize. Prior studies show that pbf1 underwent selection during maize domestication although how it affected trait change during domestication is unknown. To assay how pbf1 affects phenotypic differences between maize and teosinte, we compared nearly isogenic lines (NILs) that differ for a maize versus teosinte allele of pbf1 Kernel weight for the teosinte NIL (162mg) is slightly but significantly greater than that for the maize NIL (156mg). RNAseq data for developing kernels show that the teosinte allele of pbf1 is expressed at about twice the level of the maize allele. However, RNA and protein assays showed no difference in zein profile between the two NILs. The lower expression for the maize pbf1 allele suggests that selection may have favored this change; however, how reduced pbf1 expression alters phenotype remains unknown. One possibility is that pbf1 regulates genes other than zeins and thereby is a domestication trait. The observed drop in seed weight associated with the maize allele of pbf1 is counterintuitive but could represent a negative pleiotropic effect of selection on some other aspect of kernel composition.

Comprehensive high-resolution analysis of the role of an Arabidopsis gene family in RNA editing.

In flowering plants, mitochondrial and chloroplast mRNAs are edited by C-to-U base modification. In plant organelles, RNA editing appears to be generally a correcting mechanism that restores the proper function of the encoded product. Members of the Arabidopsis RNA editing-Interacting Protein (RIP) family have been recently shown to be essential components of the plant editing machinery. We report the use of a strand- and transcript-specific RNA-seq method (STS-PCRseq) to explore the effect of mutation or silencing of every RIP gene on plant organelle editing. We confirm RIP1 to be a major editing factor that controls the editing extent of 75% of the mitochondrial sites and 20% of the plastid C targets of editing. The quantitative nature of RNA sequencing allows the precise determination of overlapping effects of RIP factors on RNA editing. Over 85% of the sites under the influence of RIP3 and RIP8, two moderately important mitochondrial factors, are also controlled by RIP1. Previously uncharacterized RIP family members were found to have only a slight effect on RNA editing. The preferential location of editing sites controlled by RIP7 on some transcripts suggests an RNA metabolism function for this factor other than editing. In addition to a complete characterization of the RIP factors for their effect on RNA editing, our study highlights the potential of RNA-seq for studying plant organelle editing. Unlike previous attempts to use RNA-seq to analyze RNA editing extent, our methodology focuses on sequencing of organelle cDNAs corresponding to known transcripts. As a result, the depth of coverage of each editing site reaches unprecedented values, assuring a reliable measurement of editing extent and the detection of numerous new sites. This strategy can be applied to the study of RNA editing in any organism.

Comparison of public health impact of Listeria monocytogenes product-to-product and environment-to-product contamination of deli meats at retail.

This study compared the relative public health impact in deli meats at retail contaminated with Listeria monocytogenes by either (i) other products or (ii) the retail environment. Modeling was performed using the risk of listeriosis-associated deaths as a public health outcome of interest and using two deli meat products (i.e., ham and turkey, both formulated without growth inhibitors) as model systems. Based on reported data, deli meats coming to retail were assumed to be contaminated at a frequency of 0.4%. Three contamination scenarios were investigated: (i) a baseline scenario, in which no additional cross-contamination occurred at retail, (ii) a scenario in which an additional 2.3% of products were cross-contaminated at retail due to transfer of L. monocytogenes cells from already contaminated ready-to-eat deli meats, and (iii) a scenario in which an additional 2.3% of products were contaminated as a result of cross-contamination from a contaminated retail environment. By using a previously reported L. monocytogenes risk assessment model that uses product-specific growth kinetic parameters, cross-contamination of deli ham and turkey was estimated to increase the relative risk of listeriosis-associated deaths by 5.9- and 6.1-fold, respectively, for contamination from other products and by 4.9- and 5.8-fold, respectively, for contamination from the retail environment. Sensitivity and scenario analyses indicated that the frequency of cross-contamination at retail from any source (other food products or environment) was the most important factor affecting the relative risk of listeriosis-associated deaths. Overall, our data indicate that retail-level cross-contamination of ready-to-eat deli meats with L. monocytogenes has the potential to considerably increase the risk of human listeriosis cases and deaths, and thus precise estimates of cross-contamination frequency are critical for accurate risk assessments.

Quantitative risk assessment of listeriosis-associated deaths due to Listeria monocytogenes contamination of deli meats originating from manufacture and retail.

The objective of this study was to estimate the relative risk of listeriosis-associated deaths attributable to Listeria monocytogenes contamination in ham and turkey formulated without and with growth inhibitors (GIs). Two contamination scenarios were investigated: (i) prepackaged deli meats with contamination originating solely from manufacture at a frequency of 0.4% (based on reported data) and (ii) retail-sliced deli meats with contamination originating solely from retail at a frequency of 2.3% (based on reported data). Using a manufacture-to-consumption risk assessment with product-specific growth kinetic parameters (i.e., lag phase and exponential growth rate), reformulation with GIs was estimated to reduce human listeriosis deaths linked to ham and turkey by 2.8- and 9-fold, respectively, when contamination originated at manufacture and by 1.9- and 2.8-fold, respectively, for products contaminated at retail. Contamination originating at retail was estimated to account for 76 and 63% of listeriosis deaths caused by ham and turkey, respectively, when all products were formulated without GIs and for 83 and 84% of listeriosis deaths caused by ham and turkey, respectively, when all products were formulated with GIs. Sensitivity analyses indicated that storage temperature was the most important factor affecting the estimation of per annum relative risk. Scenario analyses suggested that reducing storage temperature in home refrigerators to consistently below 7 degrees C would greatly reduce the risk of human listeriosis deaths, whereas reducing storage time appeared to be less effective. Overall, our data indicate a critical need for further development and implementation of effective control strategies to reduce L. monocytogenes contamination at the retail level.

T-REX: software for the processing and analysis of T-RFLP data.

BACKGROUND: Despite increasing popularity and improvements in terminal restriction fragment length polymorphism (T-RFLP) and other microbial community fingerprinting techniques, there are still numerous obstacles that hamper the analysis of these datasets. Many steps are required to process raw data into a format ready for analysis and interpretation. These steps can be time-intensive, error-prone, and can introduce unwanted variability into the analysis. Accordingly, we developed T-REX, free, online software for the processing and analysis of T-RFLP data. RESULTS: Analysis of T-RFLP data generated from a multiple-factorial study was performed with T-REX. With this software, we were able to i) label raw data with attributes related to the experimental design of the samples, ii) determine a baseline threshold for identification of true peaks over noise, iii) align terminal restriction fragments (T-RFs) in all samples (i.e., bin T-RFs), iv) construct a two-way data matrix from labeled data and process the matrix in a variety of ways, v) produce several measures of data matrix complexity, including the distribution of variance between main and interaction effects and sample heterogeneity, and vi) analyze a data matrix with the additive main effects and multiplicative interaction (AMMI) model. CONCLUSION: T-REX provides a free, platform-independent tool to the research community that allows for an integrated, rapid, and more robust analysis of T-RFLP data.

Polarizable interaction potential for water from coupled cluster calculations. I. Analysis of dimer potential energy surface.

A six-dimensional interaction potential for the water dimer has been fitted to ab initio interaction energies computed at 2510 dimer configurations. These energies were obtained by combining the supermolecular second-order energies extrapolated to the complete basis set limit from up to quadruple-zeta quality basis sets with the contribution from the coupled-cluster method including single, double, and noniterative triple excitations computed in a triple-zeta quality basis set. All basis sets were augmented by diffuse functions and supplemented by midbond functions. The energies have been fitted using an analytic form with the induction component represented by a polarizable term, making the potential directly transferable to clusters and the bulk phase. Geometries and energies of stationary points on the potential surface agree well with the results of high-level ab initio geometry optimizations.

Polarizable interaction potential for water from coupled cluster calculations. II. Applications to dimer spectra, virial coefficients, and simulations of liquid water.

The six-dimensional CC-pol interaction potential for the water dimer was used to predict properties of the dimer and of liquid water, in the latter case after being supplemented by a nonadditive potential. All the results were obtained purely from first principles, i.e., without any fitting to experimental data. Calculations of the vibration-rotation-tunneling levels of (H(2)O)(2) and (D(2)O)(2), a very sensitive test of the potential surface, gave results in good agreement with experimental high-resolution spectra. Also the virial coefficients and properties of liquid water agree well with measured values. The present model performs better than published force fields for water in a simultaneous reproduction of experimental data for dimer spectra, virials, and properties of the liquid.

Potential energy surface for cyclotrimethylene trinitramine dimer from symmetry-adapted perturbation theory.

We present a potential energy surface (PES) for the cyclotrimethylene trinitramine (RDX) dimer obtained using symmetry-adapted perturbation theory based on the Kohn-Sham density functional theory (DFT) description of the monomers [SAPT(DFT)]. More than a thousand dimer configurations were computed using an augmented double-zeta-quality basis set supplemented by bond functions. The ab initio interaction energies were used to obtain a six-dimensional analytic fit of the interaction PES. The geometries and energies of the minima on the PES have been found from the fit. The decomposition of the PES into physical components provided by the SAPT(DFT) method has been analyzed. The PES was then used in molecular dynamics simulations of the RDX crystal. The predicted crystal density is in an excellent agreement with experiment.

Predictions of the properties of water from first principles.

A force field for water has been developed entirely from first principles, without any fitting to experimental data. It contains both pairwise and many-body interactions. This force field predicts the properties of the water dimer and of liquid water in excellent agreement with experiments, a previously elusive objective. Precise knowledge of the intermolecular interactions in water will facilitate a better understanding of this ubiquitous substance.

Interaction potential for water dimer from symmetry-adapted perturbation theory based on density functional description of monomers.

A new six-dimensional interaction potential for the water dimer has been obtained by fitting interaction energies computed at 2510 geometries using a variant of symmetry-adapted perturbation theory (SAPT) based on density functional theory (DFT) description of monomers, referred to as SAPT(DFT). The stationary points, second virial coefficient, vibration-rotation-tunneling spectrum, and structure of liquid water computed with the new potential are in very good agreement with experiment and advanced ab initio calculations, confirming the high level of accuracy provided by SAPT(DFT).

Potential energy surface for the benzene dimer and perturbational analysis of pi-pi interactions.

We present a complete 6-dimensional potential energy surface for the benzene dimer obtained using symmetry-adapted perturbation theory (SAPT) of intermolecular interactions based on Kohn-Sham's description of monomers. Ab initio calculations were performed for 491 dimer geometries in a triple-zeta-quality basis set supplemented by bond functions. An accurate analytic fit to the ab initio results has been developed and low-energy stationary points on the potential energy surface have been found. We have determined that there are three minima on the surface. Two of them, the tilted T-shape and the parallel-displaced, are nearly isoenergetic with interaction energies of -2.77 and -2.74 kcal/mol, respectively. The third minimum, a twisted edge-to-edge conformation, is significantly less attractive, with the interaction energy equal to -1.82 kcal/mol. Both the T-shape and sandwich geometries, sometimes assumed to be minima, are shown to be only saddle points. The potential energy surface is extremely flat between the two lowest minima, the barrier being only 0.10 kcal/mol above the global minimum. The second-virial coefficient obtained with the new potential agrees well with experimental results over a wide range of temperatures. The SAPT approach rigorously decomposes the interaction energy into physical components. The relative importance of these components has been analyzed.

Density-Fitting Method in Symmetry-Adapted Perturbation Theory Based on Kohn-Sham Description of Monomers.

We present a new implementation of symmetry-adapted perturbation theory of intermolecular interactions based on Kohn-Sham description of monomers. With density-fitting of molecular integrals, the scaling of the computational cost of the method is reduced from the sixth to the fifth power of the system size. Computational requirements of some operations scaling as the fifth power have also been significantly reduced. The new method allows an accurate treatment of molecules consisting of as many as a few dozen of atoms, using both nonhybrid and hybrid density functionals.