Discussion: Prioritize perennial grain development for sustainable food production and environmental benefits

Perennial grains have potential to contribute to ecological intensification of food production by enabling the direct harvest of human-edible crops without requiring annual cycles of disturbance and replanting. Studies of prototype perennial grains and other herbaceous perennials point to the ability of agroecosystems including these crops to protect water quality, enhance wildlife habitat, build soil quality, and sequester soil carbon. However, genetic improvement of perennial grain candidates has been hindered by limited investment due to uncertainty about whether the approach is viable. As efforts to develop perennial grain crops have expanded in past decades, critiques of the approach have arisen. With a recent report of perennial rice producing yields equivalent to those of annual rice over eight consecutive harvests, many theoretical concerns have been alleviated. Some valid questions remain over the timeline for new crop development, but we argue these may be mitigated by implementation of recent technological advances in crop breeding and genetics such as low-cost genotyping, genomic selection, and genome editing. With aggressive research investment in the development of new perennial grain crops, they can be developed and deployed to provide atmospheric greenhouse gas reductions.

Origin of current intermediate wheatgrass germplasm being developed for Kernza grain production

Intermediate wheatgrass (IWG, Thinopyrum intermedium [Host] Barkworth & D. R. Dewey) has been developed as a perennial grain crop to provide ecosystem services, environmental benefits, and human food. Grain and products derived from IWG varieties improved for food production have been marketed under the registered trademark, Kernza. In the 1980s, a joint breeding effort between the Rodale Institute (RI) and the Big Flats Plant Material Center used IWG plant introductions (PI) from the National Plant Germplasm System (NPGS) and recurrent phenotypic selection to improve populations of IWG with the goal of developing a perennial grain. Initial selections were provided to The Land Institute where they were subsequently improved for grain production, yet the identity of the founder material of improved, food-grade IWG has not been publicly documented. Recently recovered original documents have been used to reconstruct the early breeding program to identify the most likely 20 PIs that form the founders of modern food-grade IWG. Molecular data using genotyping-by-sequencing in current elite breeding material, remnant seed from the initial RI selections, and preserved sample material have provided supporting evidence for the historical records. The genetic origin for food-grade IWG is focused between the Black Sea and Caspian Sea in the Stavropol region of Russia, with smaller contributions likely from collections as distant as Kazakhstan in the east to Turkey in the west. This work connects the flow of germplasm and utility of NPGS PIs to present day IWG grain cultivars being developed in multiple breeding programs around the world.

Genetic architecture and QTL selection response for Kernza perennial grain domestication traits

Analysis of multi-year breeding program data revealed that the genetic architecture of an intermediate wheatgrass population was highly polygenic for both domestication and agronomic traits, supporting the use of genomic selection for new crop domestication. Perennial grains have the potential to provide food for humans and decrease the negative impacts of annual agriculture. Intermediate wheatgrass (IWG, Thinopyrum intermedium, Kernza®) is a promising perennial grain candidate that The Land Institute has been breeding since 2003. We evaluated four consecutive breeding cycles of IWG from 2016 to 2020 with each cycle containing approximately 1100 unique genets. Using genotyping-by-sequencing markers, quantitative trait loci (QTL) were mapped for 34 different traits using genome-wide association analysis. Combining data across cycles and years, we found 93 marker-trait associations for 16 different traits, with each association explaining 0.8-5.2% of the observed phenotypic variance. Across the four cycles, only three QTL showed an F_ST differentiation > 0.15 with two corresponding to a decrease in floret shattering. Additionally, one marker associated with brittle rachis was 216 bp from an ortholog of the btr2 gene. Power analysis and quantitative genetic theory were used to estimate the effective number of QTL, which ranged from a minimum of 33 up to 558 QTL for individual traits. This study suggests that key agronomic and domestication traits are under polygenic control and that molecular methods like genomic selection are needed to accelerate domestication and improvement of this new crop.

Applied phenomics and genomics for improving barley yellow dwarf resistance in winter wheat

Barley yellow dwarf is one of the major viral diseases of cereals. Phenotyping barley yellow dwarf in wheat is extremely challenging due to similarities to other biotic and abiotic stresses. Breeding for resistance is additionally challenging as the wheat primary germplasm pool lacks genetic resistance, with most of the few resistance genes named to date originating from a wild relative species. The objectives of this study were to (1) evaluate the use of high-throughput phenotyping to improve barley yellow dwarf assessment; (2) identify genomic regions associated with barley yellow dwarf resistance; and (3) evaluate the ability of genomic selection models to predict barley yellow dwarf resistance. Up to 107 wheat lines were phenotyped during each of 5 field seasons under both insecticide treated and untreated plots. Across all seasons, barley yellow dwarf severity was lower within the insecticide treatment along with increased plant height and grain yield compared with untreated entries. Only 9.2% of the lines were positive for the presence of the translocated segment carrying the resistance gene Bdv2. Despite the low frequency, this region was identified through association mapping. Furthermore, we mapped a potentially novel genomic region for barley yellow dwarf resistance on chromosome 5AS. Given the variable heritability of the trait (0.211–0.806), we obtained a predictive ability for barley yellow dwarf severity ranging between 0.06 and 0.26. Including the presence or absence of Bdv2 as a covariate in the genomic selection models had a large effect for predicting barley yellow dwarf but almost no effect for other observed traits. This study was the first attempt to characterize barley yellow dwarf using field-high-throughput phenotyping and apply genomic selection to predict disease severity. These methods have the potential to improve barley yellow dwarf characterization, additionally identifying new sources of resistance will be crucial for delivering barley yellow dwarf resistant germplasm.

Improving Wheat Yield Prediction Using Secondary Traits and High-Density Phenotyping Under Heat-Stressed Environments

A primary selection target for wheat (Triticum aestivum) improvement is grain yield. However, the selection for yield is limited by the extent of field trials, fluctuating environments, and the time needed to obtain multiyear assessments. Secondary traits such as spectral reflectance and canopy temperature (CT), which can be rapidly measured many times throughout the growing season, are frequently correlated with grain yield and could be used for indirect selection in large populations particularly in earlier generations in the breeding cycle prior to replicated yield testing. While proximal sensing data collection is increasingly implemented with high-throughput platforms that provide powerful and affordable information, efficient and effective use of these data is challenging. The objective of this study was to monitor wheat growth and predict grain yield in wheat breeding trials using high-density proximal sensing measurements under extreme terminal heat stress that is common in Bangladesh. Over five growing seasons, we analyzed normalized difference vegetation index (NDVI) and CT measurements collected in elite breeding lines from the International Maize and Wheat Improvement Center at the Regional Agricultural Research Station, Jamalpur, Bangladesh. We explored several variable reduction and regularization techniques followed by using the combined secondary traits to predict grain yield. Across years, grain yield heritability ranged from 0.30 to 0.72, with variable secondary trait heritability (0.0-0.6), while the correlation between grain yield and secondary traits ranged from -0.5 to 0.5. The prediction accuracy was calculated by a cross-fold validation approach as the correlation between observed and predicted grain yield using univariate and multivariate models. We found that the multivariate models resulted in higher prediction accuracies for grain yield than the univariate models. Stepwise regression performed equal to, or better than, other models in predicting grain yield. When incorporating all secondary traits into the models, we obtained high prediction accuracies (0.58-0.68) across the five growing seasons. Our results show that the optimized phenotypic prediction models can leverage secondary traits to deliver accurate predictions of wheat grain yield, allowing breeding programs to make more robust and rapid selections.

Nested association mapping reveals the genetic architecture of spike emergence and anthesis timing in intermediate wheatgrass

Intermediate wheatgrass (Thinopyrum intermedium) is an outcrossing, cool season grass species currently undergoing direct domestication as a perennial grain crop. Though many traits are selection targets, understanding the genetic architecture of those important for local adaptation may accelerate the domestication process. Nested association mapping (NAM) has proven useful in dissecting the genetic control of agronomic traits many crop species, but its utility in primarily outcrossing, perennial species has yet to be demonstrated. Here, we introduce an intermediate wheatgrass NAM population developed by crossing ten phenotypically divergent donor parents to an adapted common parent in a reciprocal manner, yielding 1,168 F₁ progeny from 10 families. Using genotyping by sequencing, we identified 8,003 SNP markers and developed a population-specific consensus genetic map with 3,144 markers across 21 linkage groups. Using both genomewide association mapping and linkage mapping combined across and within families, we characterized the genetic control of flowering time. In the analysis of two measures of maturity across four separate environments, we detected as many as 75 significant QTL, many of which correspond to the same regions in both analysis methods across 11 chromosomes. The results demonstrate a complex genetic control that is variable across years, locations, traits, and within families. The methods were effective at detecting previously identified QTL, as well as new QTL that align closely to the well-characterized flowering time orthologs from barley, including Ppd-H1 and Constans. Our results demonstrate the utility of the NAM population for understanding the genetic control of flowering time and its potential for application to other traits of interest.

Genomic prediction enables rapid selection of high-performing genets in an intermediate wheatgrass breeding program

In an era of constrained and depleted natural resources, perennial grains could provide sustainable food production along with beneficial ecosystem services like reduced erosion and increased atmospheric carbon capture. Intermediate wheatgrass (IWG) [Thinopyrum intermedium (Host) Barkworth & D. R. Dewey subsp. intermedium] has been undergoing continuous breeding for domestication to develop a perennial grain crop since the 1980s. As a perennial, IWG has required 2-5 yr per selection generation, but starting in 2017, genomic selection (GS) was initiated in the breeding program at The Land Institute, Salina, KS (TLI), enabling one complete cycle per year. For each cycle, ∼4,000 seedlings were profiled using genotyping-by-sequencing (GBS) and genomic estimated breeding values (GEBVs) were calculated. Selection based on GEBVs identified ∼100 individuals to advance as parents each generation, while validation populations of 1,000-1,200 genets for GS model training were also selected using the genomic relationship matrix to represent genetic diversity in each cycle. The selected parents were randomly intermated in a greenhouse crossing block to form the subsequent cycle, while the validation populations were transplanted to irrigated and nonirrigated field sites for phenotypic evaluations in the following years. For priority breeding traits of seed mass, free threshing, and nonshattering, correlations between predicted values and observed data were >.5. The realized selection differential ranged from 11-23% for selected traits, and the expected genetic gains for these traits, including spike yield, ranged from 6 to 14% per year. Genomic selection is a powerful tool to speed the domestication and development of IWG and other perennial crops with extended breeding timelines.

Development of whole-genome prediction models to increase the rate of genetic gain in intermediate wheatgrass (Thinopyrum intermedium) breeding

The development of perennial grain crops is driven by the vision of simultaneous food production and enhanced ecosystem services. Typically, perennial crops like intermediate wheatgrass (IWG)[Thinopyrum intermedium (Host) Barkworth & D.R Dewey] have low seed yield and other detrimental traits. Next-generation sequencing has made genomic selection (GS) a tractable and viable breeding method. To investigate how an IWG breeding program may use GS, we evaluated 3,658 genets over 2 yr for 46 traits to build a training population. Six statistical models were used to evaluate the non-replicated data, and a model using autoregressive order 1 (AR1) spatial correction for rows and columns combined with the genomic relationship matrix provided the highest estimates of heritability. Genomic selection models were built from 18,357 single nucleotide polymorphism markers via genotyping-by-sequencing, and a 20-fold cross-validation showed high predictive ability for all traits (r > .80). Predictive abilities improved with increased training population size and marker numbers, even with larger amounts of missing data per marker. On the basis of these results, we propose a GS breeding method that is capable of completing one cycle per year compared with a minimum of 2 yr per cycle with phenotypic selection. We estimate that this breeding approach can increase the rate of genetic gain up to 2.6× above phenotypic selection for spike yield in IWG, allowing GS to enable rapid domestication and improvement of this crop. These breeding methods should be transferable to other species with similar long breeding cycles or limited capacity for replicated observations.

Sequenced-based paternity analysis to improve breeding and identify self-incompatibility loci in intermediate wheatgrass (Thinopyrum intermedium)

KEY MESSAGE

Paternity assignment and genome-wide association analyses for fertility were applied to a Thinopyrum intermedium breeding program. A lack of progeny between combinations of parents was associated with loci near self-incompatibility genes. In outcrossing species such as intermediate wheatgrass (IWG, Thinopyrum intermedium), polycrossing is often used to generate novel recombinants through each cycle of selection, but it cannot track pollen-parent pedigrees and it is unknown how self-incompatibility (SI) genes may limit the number of unique crosses obtained. This study investigated the potential of using next-generation sequencing to assign paternity and identify putative SI loci in IWG. Using a reference population of 380 individuals made from controlled crosses of 64 parents, paternity was assigned with 92% agreement using Cervus software. Using this approach, 80% of 4158 progeny (n = 3342) from a polycross of 89 parents were assigned paternity. Of the 89 pollen parents, 82 (92%) were represented with 1633 unique full-sib families representing 42% of all potential crosses. The number of progeny per successful pollen parent ranged from 1 to 123, with number of inflorescences per pollen parent significantly correlated to the number of progeny (r = 0.54, p < 0.001). Shannon's diversity index, assessing the total number and representation of families, was 7.33 compared to a theoretical maximum of 8.98. To test our hypothesis on the impact of SI genes, a genome-wide association study of the number of progeny observed from the 89 parents identified genetic effects related to non-random mating, including marker loci located near putative SI genes. Paternity testing of polycross progeny can impact future breeding gains by being incorporated in breeding programs to optimize polycross methodology, maintain genetic diversity, and reveal genetic architecture of mating patterns.

Enhancing Crop Domestication Through Genomic Selection, a Case Study of Intermediate Wheatgrass

Perennial grains could simultaneously provide food for humans and a host of ecosystem services, including reduced erosion, minimized nitrate leaching, and increased carbon capture. Yet most of the world's food and feed is supplied by annual grains. Efforts to domesticate intermediate wheatgrass (Thinopyrumn intermedium, IWG) as a perennial grain crop have been ongoing since the 1980's. Currently, there are several breeding programs within North America and Europe working toward developing IWG into a viable crop. As new breeding efforts are established to provide a widely adapted crop, questions of how genomic and phenotypic data can be used among sites and breeding programs have emerged. Utilizing five cycles of breeding data that span 8 years and two breeding programs, University of Minnesota, St. Paul, MN, and The Land Institute, Salina, KS, we developed genomic selection (GS) models to predict IWG traits. Seven traits were evaluated with free-threshing seed, seed mass, and non-shattering being considered domestication traits while agronomic traits included spike yield, spikelets per inflorescence, plant height, and spike length. We used 6,199 genets - unique, heterozygous, individual plants - that had been profiled with genotyping-by-sequencing, resulting in 23,495 SNP markers to develop GS models. Within cycles, the predictive ability of GS was high, ranging from 0.11 to 0.97. Across-cycle predictions were generally much lower, ranging from -0.22 to 0.76. The prediction ability for domestication traits was higher than agronomic traits, with non-shattering and free threshing prediction abilities ranging from 0.27 to 0.75 whereas spike yield had prediction abilities ranging from -0.22 to 0.26. These results suggest that progress to reduce shattering and increase the percent free-threshing grain can be made irrespective of the location and breeding program. While site-specific programs may be required for agronomic traits, synergies can be achieved in rapidly improving key domestication traits for IWG. As other species are targeted for domestication, these results will aid in rapidly domesticating new crops.

Genome mapping of quantitative trait loci (QTL) controlling domestication traits of intermediate wheatgrass (Thinopyrum intermedium)

Allohexaploid (2n = 6x = 42) intermediate wheatgrass (Thinopyrum intermedium), abbreviated IWG, is an outcrossing perennial grass belonging to the tertiary gene pool of wheat. Perenniality would be valuable option for grain production, but attempts to introgress this complex trait from wheat-Thinopyrum hybrids have not been commercially successful. Efforts to breed IWG itself as a dual-purpose forage and grain crop have demonstrated useful progress and applications, but grain yields are significantly less than wheat. Therefore, genetic and physical maps have been developed to accelerate domestication of IWG. Herein, these maps were used to identify quantitative trait loci (QTLs) and candidate genes associated with IWG grain production traits in a family of 266 full-sib progenies derived from two heterozygous parents, M26 and M35. Transgressive segregation was observed for 17 traits related to seed size, shattering, threshing, inflorescence capacity, fertility, stem size, and flowering time. A total of 111 QTLs were detected in 36 different regions using 3826 genotype-by-sequence markers in 21 linkage groups. The most prominent QTL had a LOD score of 15 with synergistic effects of 29% and 22% over the family means for seed retention and percentage of naked seeds, respectively. Many QTLs aligned with one or more IWG gene models corresponding to 42 possible domestication orthogenes including the wheat Q and RHT genes. A cluster of seed-size and fertility QTLs showed possible alignment to a putative Z self-incompatibility gene, which could have detrimental grain-yield effects when genetic variability is low. These findings elucidate pathways and possible hurdles in the domestication of IWG.

Combining High-Throughput Phenotyping and Genomic Information to Increase Prediction and Selection Accuracy in Wheat Breeding

Genomics and phenomics have promised to revolutionize the field of plant breeding. The integration of these two fields has just begun and is being driven through big data by advances in next-generation sequencing and developments of field-based high-throughput phenotyping (HTP) platforms. Each year the International Maize and Wheat Improvement Center (CIMMYT) evaluates tens-of-thousands of advanced lines for grain yield across multiple environments. To evaluate how CIMMYT may utilize dynamic HTP data for genomic selection (GS), we evaluated 1170 of these advanced lines in two environments, drought (2014, 2015) and heat (2015). A portable phenotyping system called 'Phenocart' was used to measure normalized difference vegetation index and canopy temperature simultaneously while tagging each data point with precise GPS coordinates. For genomic profiling, genotyping-by-sequencing (GBS) was used for marker discovery and genotyping. Several GS models were evaluated utilizing the 2254 GBS markers along with over 1.1 million phenotypic observations. The physiological measurements collected by HTP, whether used as a response in multivariate models or as a covariate in univariate models, resulted in a range of 33% below to 7% above the standard univariate model. Continued advances in yield prediction models as well as increasing data generating capabilities for both genomic and phenomic data will make these selection strategies tractable for plant breeders to implement increasing the rate of genetic gain.

Epidemiology of morbidity and mortality in US and Canadian recreational scuba diving

OBJECTIVES

This study investigates morbidity and mortality suffered by divers in the USA and Canada.

STUDY DESIGN

Prospectively recruited probability-weighted sample for estimating the national burden of injury and a weighted retrospective survey for estimating exposure.

METHODS

The National Electronic Surveillance System and Canadian Hospitals Injury Reporting and Prevention Program (CHIRPP) were searched for scuba diving injuries. The Divers Alert Network diving fatality database was searched for deaths, and Sports and Fitness Industry Association estimates for diving were obtained from annual surveys.

RESULTS

In the USA, there were an estimated 1394 emergency department (ED) presentations annually for scuba-related injuries. The majority (80%) were treated and/or released. There were an estimated 306 million dives made by the US residents 2006-2015 and concurrently 563 recreational diving deaths, a fatality rate of 0.18 per 10⁵ dives and 1.8 per 10⁵ diver-years. There were 658 diving deaths in the US 2006-2015 and 13,943 ED presentations for scuba injuries, giving a ratio of 47 diving deaths in the USA for every 1000 ED presentations. There were 98 cases of scuba-related injuries identified in the CHIRPP data. The prevalence of scuba-related injuries for patients aged 3-17 years was 1.5 per 10⁵ cases, and the prevalence of scuba-related injuries to patients 18-62 years was 16.5 per 10⁵ cases.

DISCUSSION

In Canada and the USA, only one out of every 10,000 ED presentations is due to a scuba-related injury. That there are 47 deaths for every 1000 ED presentations for scuba injuries speaks to the relatively unforgiving environment in which scuba diving takes place. For 1.8 deaths per million recreational dives, mortality in scuba diving is nonetheless relatively low.

At-a-glance, Emergency department surveillance of thermal burns and scalds, electronic Canadian Hospitals Injury Reporting and Prevention Program, 2013

INTRODUCTION

Although fatality and hospitalization rates for burns in Canada have declined over time, less serious cases still commonly present to the emergency department (ED).

METHODS

The Canadian Hospitals Injury Reporting and Prevention Program (CHIRPP) is an injury and poisoning surveillance system administered by the Public Health Agency of Canada, operating in emergency departments of 17 hospitals.

RESULTS

Overall, cases reported in 2013 were scalds and contact burns from hot objects. The leading direct causes of scalds were hot beverages and hot water. The leading causes of contact burns were stoves/ovens and fireplaces/accessories. While the overall proportion of burns was highest among females, males comprised a higher proportion of burns from all mechanisms except scalds.

Status report - The Canadian Hospitals Injury Reporting and Prevention Program: a dynamic and innovative injury surveillance system

This status report on the Canadian Hospitals Injury Reporting and Prevention Program (CHIRPP), an emergency department-based injury and poisoning surveillance system, describes the result of migrating from a centralized data entry and coding process to a decentralized process, the web-based eCHIRPP system, in 2011. This secure system is improving the CHIRPP's overall flexibility and timeliness, which are key attributes of an effective surveillance system. The integrated eCHIRPP platform enables near real-time data entry and access, has user-friendly data management and analysis tools, and allows for easier communication and connectivity across the CHIRPP network through an online collaboration centre. Current pilot testing of automated data monitoring and trend analysis tools-designed to monitor and flag incoming data according to predefined criteria (for example, a new consumer product)-is revealing eCHIRPP's potential for providing early warnings of new hazards, issues and trends.

Molecular-scale remnants of the liquid-gas transition in supercritical polar fluids

An electronically coarse-grained model for water reveals a persistent vestige of the liquid-gas transition deep into the supercritical region. A crossover in the density dependence of the molecular dipole arises from the onset of nonpercolating hydrogen bonds. The crossover points coincide with the Widom line in the scaling region but extend farther, tracking the heat capacity maxima, offering evidence for liquidlike and gaslike state points in a "one-phase" fluid. The effect is present even in dipole-limit models, suggesting that it is common for all molecular liquids exhibiting dipole enhancement in the liquid phase.

Development of a PCR-free electrochemical point of care test for clinical detection of methicillin resistant Staphylococcus aureus (MRSA)

An MRSA assay requiring neither labeling nor amplification of target DNA has been developed. Sequence specific binding of fragments of bacterial genomic DNA is detected at femtomolar concentrations using electrochemical impedance spectroscopy (EIS). This has been achieved using systematic optimisation of probe chemistry (PNA self-assembled monolayer film on gold electrode), electrode film structure (the size and nature of the chemical spacer) and DNA fragmentation, as these are found to play an important role in assay performance. These sensitivity improvements allow the elimination of the PCR step and DNA labeling and facilitate the development of a simple and rapid point of care test for MRSA. Assay performance is then evaluated and specific direct detection of the MRSA diagnostic mecA gene from genomic DNA, extracted directly from bacteria without further treatment is demonstrated for bacteria spiked into saline (10(6) cells per mL) on gold macrodisc electrodes and into human wound fluid (10(4) cells per mL) on screen printed gold electrodes. The latter detection level is particularly relevant to clinical requirements and point of care testing where the general threshold for considering a wound to be infected is 10(5) cells per mL. By eliminating the PCR step typically employed in nucleic acid assays, using screen printed electrodes and achieving sequence specific discrimination under ambient conditions, the test is extremely simple to design and engineer. In combination with a time to result of a few minutes this means the assay is well placed for use in point of care testing.

Nonlinear optical methods for cellular imaging and localization

Of all the ways in which complex materials (including many biological systems) can be explored, imaging is perhaps the most powerful because delivering high information content directly. This is particular relevant in aspects of cellular localization where the physical proximity of molecules is crucial in biochemical processes. A great deal of effort in imaging has been spent on enabling chemically selective imaging so that only specific features are revealed. This is almost always achieved by adding fluorescent chemical labels to specific molecules. Under appropriate illumination conditions only the molecules (via their labels) will be visible. The technique is simple and elegant but does suffer from fundamental limitations: (1) the fluorescent labels may fade when illuminated (a phenomenon called photobleaching) thereby constantly decreasing signal contrast over the course of image acquisition. To combat photobleaching one must reduce observation times or apply unfavourably low excitation levels all of which reduce the information content of images; (2) the fluorescent species may be deactivated by various environmental factors (the general term is fluorescence quenching); (3) the presence of fluorescent labels may introduce unexpected complications or may interfere with processes of interest (4) Some molecules of interest cannot be labelled. In these circumstances we require a fundamentally different strategy. One of the most promising alternative is based on a technique called Coherent Anti-Stokes Raman scattering (CARS). CARS is a fundamentally more complex process than is fluorescence and the experimental procedures and optical systems required to deliver high quality CARS images are intricate. However, the rewards are correspondingly very high: CARS probes the chemically distinct vibrations of the constituent molecules in a complex system and is therefore also chemically selective as are fluorescence-based methods. Moreover,the potentially severe problems of fluorescence bleaching and quenching are circumvented and high-resolution three dimensional images can be obtained on completely unlabelled specimens. We review here aspects of CARS and Multiphoton fluorescence techniques to cellular localization and measurement.

Electronically coarse-grained model for water

We introduce an electronically coarse-grained description of water representing all long range, many-body electronic responses via an embedded quantum oscillator. Leading-order response coefficients and gas phase electrostatic moments are exactly reproduced. Molecular dynamics, using electronic path integral sampling, shows that this framework is sufficient for a realistic liquid to emerge naturally with transferability extending further to nonambient state points and to the free water surface. The model allows the strength of many-body dispersion and polarization to be adjusted independently and these are found to have significant effects on the condensed phase.

Report summary - Injury in Review, 2012 Edition: Spotlight on Road and Transport Safety

Injury in Review, 2012 Edition: Spotlight on Road and Transport Safety, the first national public health report of its kind, synthesizes road- and transport-related injury statistics from a variety of sources. It profiles injury patterns among Canadians aged up to 24 years, explains risks and protective factors, and makes recommendations for action. The findings inform the development of targeted injury prevention efforts.

Report summary - injury in review, 2012 edition: spotlight on road and transport safety

Injury in Review, 2012 Edition: Spotlight on Road and Transport Safety, the first national public health report of its kind, synthesizes road- and transport-related injury statistics from a variety of sources. It profiles injury patterns among Canadians aged up to 24 years, explains risks and protective factors, and makes recommendations for action. The findings inform the development of targeted injury prevention efforts.

Structural properties of liquid N-methylacetamide via ab initio, path integral, and classical molecular dynamics

In order to better understand the physical interactions that stabilize protein secondary structure, the neat liquid state of a peptidic fragment, N-methylacetamide (NMA), was studied using computer simulation. Three different descriptions of the molecular liquid were examined: an empirical force field treatment with classical nuclei, an empirical force field treatment with quantum mechanical nuclei, and an ab initio density functional theory (DFT) treatment. The DFT electronic structure was evaluated using the BLYP approximate functional and a plane wave basis set. The different physical effects probed by the three models, such as quantum dispersion, many-body polarization, and nontrivial charge distributions on the liquid properties, were compared. Much of the structural ordering in the liquid is characterized by hydrogen bonded chains of NMA molecules. Modest structural differences are present among the three models of liquid NMA. The average molecular dipole in the liquid under the ab initio treatment, however, is enhanced by 60% over the gas phase value.

Structure and hydrogen bonding in neat N-methylacetamide: classical molecular dynamics and Raman spectroscopy studies of a liquid of peptidic fragments

The results of classical molecular dynamics (MD) simulations and Raman spectroscopy studies of neat liquid N-methylacetamide (NMA), the simplest model system relevant to the peptides, are reported as a function of temperature and pressure. The MD simulations predict that near ambient conditions, the molecules form a hydrogen bond network consisting primarily of linear chains. Both the links between molecules within the hydrogen-bonded chains and the associations between chains are stabilized by weak methyl-donated "improper" hydrogen bonds. The three-dimensional structural motifs observed in the liquid show some similarity to protein beta-sheets. The temperature and pressure dependence of the hydrogen bond network, as probed by the mode frequency of the experimentally determined amide-I Raman band, blue shifts on heating and red shifts under compression, respectively, suggesting weakened and enhanced hydrogen bonding in response to temperature and pressure increases. Disruption of the hydrogen-bonding network is clearly observed in the simulation data as temperature is increased, whereas the improper hydrogen bonding is enhanced under compression to reduce the energetic cost of increasing the packing fraction. Because of the neglect of polarizability in the molecular model, the computed dielectric constant is underestimated compared to the experimental value, indicating that the simulation may underestimate dipolar coupling in the liquid.

Structure of Aqueous Proline via Parallel Tempering Molecular Dynamics and Neutron Diffraction

The structure of aqueous L-proline amino acid has been the subject of much debate centering on the validity of various proposed models, differing widely in the extent to which local and long-range correlations are present. Here, aqueous proline is investigated by atomistic, replica exchange molecular dynamics simulations, and the results are compared to neutron diffraction and small angle neutron scattering (SANS) data, which have been reported recently (McLain, S.; Soper, A.; Terry, A.; Watts, A. J. Phys. Chem. B 2007, 111, 4568). Comparisons between neutron experiments and simulation are made via the static structure factor S(Q) which is measured and computed from several systems with different H/D isotopic compositions at a concentration of 1:20 molar ratio. Several different empirical water models (TIP3P, TIP4P, and SPC/E) in conjunction with the CHARMM22 force field are investigated. Agreement between experiment and simulation is reasonably good across the entire Q range although there are significant model-dependent variations in some cases. In general, agreement is improved slightly upon application of approximate quantum corrections obtained from gas-phase path integral simulations. Dimers and short oligomeric chains formed by hydrogen bonds (frequently bifurcated) coexist with apolar (hydrophobic) contacts. These emerge as the dominant local motifs in the mixture. Evidence for long-range association is more equivocal: No long-range structures form spontaneously in the MD simulations, and no obvious low-Q signature is seen in the SANS data. Moreover, associations introduced artificially to replicate a long-standing proposed mesoscale structure for proline correlations as an initial condition are annealed out by parallel tempering MD simulations. However, some small residual aggregates do remain, implying a greater degree of long-range order than is apparent in the SANS data.

The stability and characteristics of a DNA Holliday junction switch

A Holliday junction (HJ) consists of four DNA double helices, with a branch point discontinuity at the intersection of the component strands. At low ionic strength, the HJ adopts an open conformation, with four widely spaced arms, primarily due to strong electrostatic repulsion between the phosphate groups on the backbones. At high ionic strength, screening of this repulsion induces a switch to a more compact (closed) junction conformation. Fluorescent labelling with dyes placed on the HJ arms allows this conformational switch to be detected optically using fluorescence resonance energy transfer (FRET), producing a sensitive fluorescent output of the switch state. This paper presents a systematic and quantitative survey of the switch characteristics of such a labelled HJ. A short HJ (arm length 8 bp) is shown to be prone to dissociation at low switching ion concentration, whereas an HJ of arm length 12 bp is shown to be stable over all switching ion concentrations studied. The switching characteristics of this HJ have been systematically and quantitatively studied for a variety of switching ions, by measuring the required ion concentration, the sharpness of the switching transition and the fluorescent output intensity of the open and closed states. This stable HJ is shown to have favourable switch characteristics for a number of inorganic switching ions, making it a promising candidate for use in nanoscale biomolecular switch devices.

Time-Resolved FRET and FLIM of Four-way DNA Junctions

Conformational transitions in a 4-way DNA junction when titrated with ionic solutions are studied using time-resolved fluorescence resonance energy transfer. Parameters characterising the transition in terms of critical ion concentration (c1/2) and the Hill coefficient for ion binding are obtained by fitting a simple two-state model using steady-state spectra. Data obtained from a fluorescence lifetime plate reader and analysed by fitting a single exponential to donor fluorescence lifetime decays are shown to be in good agreement with the parameters obtained from steady-state measurements. Fluorescence lifetimes, however, offer advantages, particularly in being independent of fluorophore concentration, output intensity, inhomogeneity in the excitation source and output wavelength. We demonstrate preliminary FRET-FLIM images of DNA junction solutions obtained using a picosecond gated CCD which are in agreement with results from a fluorescence lifetime plate reader. The results suggest that time-resolved FRET-FLIM is sensitive to subtle structural changes and may be useful in assays based on 4-way DNA junctions.

Segregation in aqueous methanol enhanced by cooling and compression

Molecular segregation in methanol-water mixtures is studied across a wide concentration range as a function of temperature and pressure. Cluster distributions obtained from both neutron diffraction and molecular dynamics simulations point to significantly enhanced segregation as the mixtures are cooled or compressed. This evolution toward greater molecular heterogenity in the mixture accounts for the observed changes in the water-water radial distribution function and there are indications also of a change in the topology of the water clusters. The observed behavior is consistent with an approach to an upper critical solution point. Such a point would appear to be "hidden" below the freezing line, thereby precluding observation of the two-fluid region.

Methanol-water solutions: A bi-percolating liquid mixture

An extensive series of neutron diffraction experiments and molecular dynamics simulations has shown that mixtures of methanol and water exhibit extended structures in solution despite the components being fully miscible in all proportions. Of particular interest is a concentration region (methanol mole fraction between 0.27 and 0.54) where both methanol and water appear to form separate, percolating networks. This is the concentration range where many transport properties and thermodynamic excess functions reach extremal values. The observed concentration dependence of several of these material properties of the solution may therefore have a structural origin.

Single Molecule Fluorescence Imaging and Its Application to the Study of DNA Condensation

Single molecule fluorescence imaging incorporated with optical tweezers and a laminar flow cell has been used to monitor the kinetic process of DNA condensation induced by spermidine. It was found that at least two steps were involved in the condensation process of the hydrodynamically-stretched linear DNA; a lag period followed by a rapid collapse of DNA. The lag time increased with the flow speed and the collapse time remained short within the range of the flow speed studied. The effect of salt concentration on the condensation process was examined, and the results suggest that the longer lag time observed in the higher salt buffer probably results from the displacement of bound cations and rearrangement of spermidine on the DNA. The flow-speed dependence of the lag time suggests that a nucleation event at the free end of the DNA, i.e. formation of a loop, may play a vital role in the kinetic process of condensation.

Molecular segregation observed in a concentrated alcohol-water solution

When a simple alcohol such as methanol or ethanol is mixed with water, the entropy of the system increases far less than expected for an ideal solution of randomly mixed molecules. This well-known effect has been attributed to hydrophobic headgroups creating ice-like or clathrate-like structures in the surrounding water, although experimental support for this hypothesis is scarce. In fact, an increasing amount of experimental and theoretical work suggests that the hydrophobic headgroups of alcohol molecules in aqueous solution cluster together. However, a consistent description of the details of this self-association is lacking. Here we use neutron diffraction with isotope substitution to probe the molecular-scale structure of a concentrated alcohol water mixture (7:3 molar ratio). Our data indicate that most of the water molecules exist as small hydrogen-bonded strings and clusters in a 'fluid' of close-packed methyl groups, with water clusters bridging neighbouring methanol hydroxyl groups through hydrogen bonding. This behaviour suggests that the anomalous thermodynamics of water alcohol systems arises from incomplete mixing at the molecular level and from retention of remnants of the three-dimensional hydrogen-bonded network structure of bulk water.

Colloid-liquid-crystal composites: An unusual soft solid

Mixing model colloidal particles with a thermotropic nematic liquid crystal results in a soft solid with significant storage modulus (G' approximately 10(3)-10(5) Pa). The soft solid comprises a network of particle aggregates, formed by the exclusion of particles from emergent nematic domains as the mixture is cooled below the isotropic-nematic transition. The unusually high storage modulus of the colloid-liquid-crystal composites may be due to the local frustration of nematic order within the particle aggregates. The birefringent soft solid is potentially important as a switchable electro-optical material that can be readily handled and processed.

Mitral regurgitation after atrioventricular node catheter ablation for atrial fibrillation and heart failure: acute hemodynamic features

BACKGROUND

Radiofrequency catheter ablation of the atrioventricular node and pacemaker insertion have been associated with occasional development of mitral regurgitation (MR). Ventricular pacing might result in MR if (1) left ventricular (LV) compliance is decreased and/or (2) mitral valve leaflet apposition is disturbed. We studied acute hemodynamic changes resulting from initiation of ventricular pacing in patients undergoing ablation.

METHODS AND RESULTS

Thirteen patients (10 men and 3 women) with a mean age of 73.4 +/- 8. 6 years, with chronic atrial fibrillation and congestive heart failure, had permanent ventricular pacemaker insertion with lead placement at the right ventricular (RV) apex. The following hemodynamic recordings were obtained before ablation (during atrial fibrillation) and then immediately after ablation (during RV pacing): heart rate, mean arterial pressure, LV end-diastolic pressure (LVEDP), mean pulmonary capillary wedge pressure, V-wave amplitude, and cardiac index. Presence of MR was assessed by V-wave amplitude and the results of LV angiography. In patients who had MR, recordings were also obtained during temporary ventricular pacing from the RV outflow tract (RVOT). As a group there were no significant changes in any hemodynamic measurement. Before ablation, mild MR by LV angiogram was present in 5 patients, but none had large V-wave amplitude. After ablation, mild MR was present by LV angiogram in 6 patients, and in 3 of these patients large V-wave amplitude developed (mean amplitude 42.7 +/- 2.2 mm Hg; assigned to group 1). This was associated with a rise in LVEDP in 1 patient (consistent with reduced LV compliance), but LVEDP was unchanged in the other 2 patients (suggesting abnormal mitral valve leaflet apposition). All patients in group 1 exhibited a fall in V-wave amplitude when the pacing site was moved to the RVOT.

CONCLUSIONS

Both reduced LV compliance and disturbed mitral valve leaflet apposition contribute to MR after ablation. MR is reduced by pacing from the RVOT.

Effect of antiphospholipid antibodies in women undergoing in-vitro fertilization: role of heparin and aspirin

To describe the prevalence of antiphospholipid antibodies in women undergoing in-vitro fertilization (IVF) and to determine if heparin and aspirin affect implantation rates, 191 women with a history of infertility undergoing IVF were prospectively tested for antiphospholipid antibodies. This was a two-centre, non-randomized comparison of women with positive antiphospholipid antibodies receiving heparin and aspirin versus standard treatment. An enzyme-linked immunosorbent assay, with referenced standards and known positive and negative sera on each plate, was utilized to measure antibodies to cardiolipin, phosphatidylinositol, phosphatidylglycerol, phosphatidylserine and phosphatidylethanolamine. Statistical analyses of results included analysis of variance and Fisher's two-tailed exact test. Antiphospholipid antibodies were detected in 18.8% of patients undergoing IVF compared with only 5.5% in the 200 normal controls, 26% in 200 women with recurrent pregnancy loss, and 32% in 200 women with systemic lupus erythematosus. In conclusion, antiphospholipid antibodies were found more frequently in women undergoing IVF than in the normal control population. Although implantation rates appeared higher in the group of women treated with heparin and aspirin, no statistically significant differences were detected in implantation, pregnancy and ongoing pregnancy rates between those who received standard therapy and those treated with heparin and aspirin.