Genome-wide SNP discovery in walnut with an AGSNP pipeline updated for SNP discovery in allogamous organisms

ABSTRACT: BACKGROUND: A genome-wide set of single nucleotide polymorphisms (SNPs) is a valuable resource in genetic research and breeding and is usually developed by re-sequencing a genome. If a genome sequence is not available, an alternative strategy must be used. We previously reported the development of a pipeline (AGSNP) for genome-wide SNP discovery in coding sequences and other single-copy DNA without a complete genome sequence in self-pollinating (autogamous) plants. Here we updated this pipeline for SNP discovery in outcrossing (allogamous) species and demonstrated its efficacy in SNP discovery in walnut (Juglans regia L.). RESULTS: The first step in the original implementation of the AGSNP pipeline was the construction of a reference sequence and the identification of single-copy sequences in it. To identify single-copy sequences, multiple genome equivalents of short SOLiD reads of another individual were mapped to shallow genome coverage of long Sanger or Roche 454 reads making up the reference sequence. The relative depth of SOLiD reads was used to filter out repeated sequences from single-copy sequences in the reference sequence. The second step was a search for SNPs between SOLiD reads and the reference sequence. Polymorphism within the mapped SOLiD reads would have precluded SNP discovery; hence both individuals had to be homozygous. The AGSNP pipeline was updated here for using SOLiD or other type of short reads of a heterozygous individual for these two principal steps. A total of 32.6X walnut genome equivalents of SOLiD reads of vegetatively propagated walnut scion cultivar 'Chandler' were mapped to 48,661 'Chandler' bacterial artificial chromosome (BAC) end sequences (BESs) produced by Sanger sequencing during the construction of a walnut physical map. A total of 22,799 putative SNPs were initially identified. A total of 6,000 Infinium II type SNPs evenly distributed along the walnut physical map were selected for the construction of an Infinium BeadChip, which was used to genotype a walnut mapping population having 'Chandler' as one of the parents. Genotyping results were used to adjust the filtering parameters of the updated AGSNP pipeline. With the adjusted filtering criteria, 69.6% of SNPs discovered with the updated pipeline were real and could be mapped on the walnut genetic map. A total of 13,439 SNPs were discovered by BES re-sequencing. BESs harboring SNPs were in 677 FPC contigs covering 98% of the physical map of the walnut genome. CONCLUSION: The updated AGSNP pipeline is a versatile SNP discovery tool for a high-throughput, genome-wide SNP discovery in both autogamous and allogamous species. With this pipeline, a large set of SNPs were identified in a single walnut cultivar.     

Net primary production in three bioenergy crop systems following land conversion

Aims Identifying the amount of production and the partitioning to above- and belowground biomass is generally the first step toward selecting bioenergy systems. There are very few existing studies on the dynamics of production following land conversion. The objectives of this study were to (i) determine the differences in aboveground net primary production (ANPP), belowground net primary production (BNPP), shoot-to-root ratio (S:R) and leaf area index in three bioenergy crop systems and (ii) evaluate the production of these three systems in two different land use conversions.Methods This investigation included biometric analysis of NPP on three agricultural sites converted from conservation reserve program (CRP) management to bioenergy crop production (corn, switchgrass and prairie mix) and three sites converted from traditional agriculture production to bioenergy crop production.Important findings The site converted from conventional agriculture produced smaller ANPP in corn (19.03±1.90 standard error [SE] Mg ha^-1 year^-1) than the site converted from CRP to corn (24.54±1.43 SE Mg ha^-1 year^-1). The two land conversions were similar in terms of ANPP for switchgrass (4.88±0.43 SE for CRP and 2.04±0.23 SE Mg ha^-1 year^-1 for agriculture) and ANPP for prairie mix (4.70±0.50 SE for CRP and 3.38±0.33 SE Mg ha^-1 year^-1 for agriculture). The BNPP at the end of the growing season in all the bioenergy crop systems was not significantly different (P = 0.75, N = 8).     

Acute vascular responses to isometric handgrip exercise and effects of training in persons medicated for hypertension

Previous work from our laboratory demonstrated that isometric handgrip (IHG) training improved local, endothelium-dependent vasodilation in medicated hypertensives [McGowan CL (PhD Thesis), 2006; McGowan et al. Physiologist 47: 285, 2004]. We investigated whether changes in the capacity of smooth muscle to dilate (regardless of endothelial factors) influenced this training-induced change, and we examined the acute vascular responses to a single bout of IHG. Seventeen subjects performed four 2-min unilateral IHG contractions at 30% of maximal voluntary effort, three times a week for 8 wk. Pre- and posttraining, brachial artery flow-mediated dilation (FMD, an index of endothelium-dependent vasodilation) and nitroglycerin-mediated maximal vasodilation (an index of endothelium-independent vasodilation) were measured in the exercised arm by using ultrasound before and immediately after acute IHG exercise. IHG training resulted in improved resting brachial FMD (P < 0.01) and no change in nitroglycerin-mediated maximal vasodilation. Pre- and posttraining, brachial artery FMD decreased following an acute bout of IHG exercise (normalized to peak shear rate, pre-, before IHG exercise: 0.01 +/- 0.002, after IHG exercise: 0.008 +/- 0.002%/s(-1); post-, before IHG exercise: 0.020 +/- 0.003, after IHG exercise: 0.010 +/- 0.003%/s(-1); P < 0.01). Posttraining, resting brachial artery FMD improved yet nitroglycerin-mediated maximal vasodilation was unchanged in persons medicated for hypertension. This suggests that the training-induced improvements in the resting brachial artery FMD were not due to underlying changes in the forearm vasculature. Acute IHG exercise attenuated brachial artery FMD, and although this impairment may be interpreted as hazardous to medicated hypertensives with already dysfunctional endothelium, the effects appear transient as repeated exposure to the IHG stimulus improved resting endothelium-dependent vasodilation.     

PERK is responsible for the increased phosphorylation of eIF2alpha and the severe inhibition of protein synthesis after transient global brain ischemia

Reperfusion after global brain ischemia results initially in a widespread suppression of protein synthesis in neurons that is due to inhibition of translation initiation as a result of the phosphorylation of the alpha-subunit of eukaryotic initiation factor 2 (eIF2). To address the role of the eIF2alpha kinase RNA-dependent protein kinase-like endoplasmic reticulum kinase (PERK) in the reperfused brain, transgenic mice with a targeted disruption of the Perk gene were subjected to 20 min of forebrain ischemia followed by 10 min of reperfusion. In wild-type mice, phosphorylated eIF2alpha was detected in the non-ischemic brain and its levels were elevated threefold after 10 min of reperfusion. Conversely, there was no phosphorylated eIF2alpha detected in the non-ischemic transgenic mice and there was no sizeable rise in phosphorylated eIF2alpha levels in the forebrain after ischemia and reperfusion. Moreover, there was a substantial rescue of protein translation in the reperfused transgenic mice. Neither group showed any change in total eIF2alpha, phosphorylated eukaryotic elongation factor 2 or total eukaryotic elongation factor 2 levels. These data demonstrate that PERK is responsible for the large increase in phosphorylated eIF2alpha and the suppression of translation early in reperfusion after transient global brain ischemia.     

Insular Organization of Gene Space in Grass Genomes

Wheat and maize genes were hypothesized to be clustered into islands but the hypothesis was not statistically tested. The hypothesis is statistically tested here in four grass species differing in genome size, Brachypodium distachyon, Oryza sativa, Sorghum bicolor, and Aegilops tauschii. Density functions obtained under a model where gene locations follow a homogeneous Poisson process and thus are not clustered are compared with a model-free situation quantified through a non-parametric density estimate. A simple homogeneous Poisson model for gene locations is not rejected for the small O. sativa and B. distachyon genomes, indicating that genes are distributed largely uniformly in those species, but is rejected for the larger S. bicolor and Ae. tauschii genomes, providing evidence for clustering of genes into islands. It is proposed to call the gene islands “gene insulae” to distinguish them from other types of gene clustering that have been proposed. An average S. bicolor and Ae. tauschii insula is estimated to contain 3.7 and 3.9 genes with an average intergenic distance within an insula of 2.1 and 16.5 kb, respectively. Inter-insular distances are greater than 8 and 81 kb and average 15.1 and 205 kb, in S. bicolor and Ae. tauschii, respectively. A greater gene density observed in the distal regions of the Ae. tauschii chromosomes is shown to be primarily caused by shortening of inter-insular distances. The comparison of the four grass genomes suggests that gene locations are largely a function of a homogeneous Poisson process in small genomes. Nonrandom insertions of LTR retroelements during genome expansion creates gene insulae, which become less dense and further apart with the increase in genome size. High concordance in relative lengths of orthologous intergenic distances among the investigated genomes including the maize genome suggests functional constraints on gene distribution in the grass genomes.     

Induction of humoral and CD8+ T cell responses are required for protection against lethal Ebola virus infection

Ebola virus (EBOV)-like particles (eVLP), composed of the EBOV glycoprotein and matrix viral protein (VP)40 with a lipid membrane, are a highly efficacious method of immunization against EBOV infection. The exact requirements for immunity against EBOV infection are poorly defined at this time. The goal of this work was to determine the requirements for EBOV immunity following eVLP vaccination. Vaccination of BALB/c or C57BL/6 mice with eVLPs in conjunction with QS-21 adjuvant resulted in mixed IgG subclass responses, a Th1-like memory cytokine response, and protection from lethal EBOV challenge. Further, this vaccination schedule led to the generation of both CD4(+) and CD8(+) IFN-gamma(+) T cells recognizing specific peptides within glycoprotein and VP40. The transfer of both serum and splenocytes, but not serum or splenocytes alone, from eVLP-vaccinated mice conferred protection against lethal EBOV infection in these studies. B cells were required for eVLP-mediated immunity to EBOV because B cell-deficient mice vaccinated with eVLPs were not protected from lethal EBOV challenge. We also found that CD8(+), but not CD4(+), T cells are absolutely required for eVLP-mediated protection against EBOV infection. Further, eVLP-induced protective mechanisms were perforin-independent, but IFN-gamma-dependent. Taken together, both EBOV-specific humoral and cytotoxic CD8(+) T cell responses are critical to mediate protection against filoviruses following eVLP vaccination.     

Discovery and mapping of wheat Ph1 suppressors

Pairing between wheat (Triticum turgidum and T. aestivum) homeologous chromosomes is prevented by the expression of the Ph1 locus on the long arm of chromosome 5B. The genome of Aegilops speltoides suppresses Ph1 expression in wheat x Ae. speltoides hybrids. Suppressors with major effects were mapped as Mendelian loci on the long arms of Ae. speltoides chromosomes 3S and 7S. The chromosome 3S locus was designated Su1-Ph1 and the chromosome 7S locus was designated Su2-Ph1. A QTL with a minor effect was mapped on the short arm of chromosome 5S and was designated QPh.ucd-5S. The expression of Su1-Ph1 and Su2-Ph1 increased homeologous chromosome pairing in T. aestivum x Ae. speltoides hybrids by 8.4 and 5.8 chiasmata/cell, respectively. Su1-Ph1 was completely epistatic to Su2-Ph1, and the two genes acting together increased homeologous chromosome pairing in T. aestivum x Ae. speltoides hybrids to the same level as Su1-Ph1 acting alone. QPh.ucd-5S expression increased homeologous chromosome pairing by 1.6 chiasmata/cell in T. aestivum x Ae. speltoides hybrids and was additive to the expression of Su2-Ph1. It is hypothesized that the products of Su1-Ph1 and Su2-Ph1 affect pairing between homeologous chromosomes by regulating the expression of Ph1 but the product of QPh.ucd-5S may primarily regulate recombination between homologous chromosomes.     

Influence of ventrolateral surface of medulla on reflex tracheal constriction

To assess the role of structures located superficially near the ventrolateral surface of the medulla on the reflex constriction of tracheal smooth muscle that occurs when airway and pulmonary receptors are stimulated mechanically or chemically, experiments were conducted in alpha-chloralose-anesthetized, paralyzed, and artificially ventilated cats. Pressure changes within a bypassed segment of the trachea were used as an index of alterations smooth muscle tone. The effects of focal cooling of the intermediate areas or topically applied lidocaine on the ventral surface of the medulla on the response of the trachea to mechanical and chemical stimulation of airway receptors were examined. Atropine abolished tracheal constriction induced by mechanical stimulation of the carina or aerosolized histamine, showing that the responses were mediated over vagal pathways. Moderate cooling of the intermediate area (20 degrees C) or local application of lidocaine significantly decreased the tracheal constrictive response to mechanical activation of airway receptors. Furthermore, when the trachea was constricted by histamine, cooling of the intermediate area significantly diminished the increased tracheal tone, whereas rewarming restored tracheal tone to the previous level. These findings suggest that under the conditions of the experiments the ventral surface of the medulla plays an important role in constriction of the trachea by inputs from intrapulmonary receptors and in the modulation of parasympathetic outflow to airway smooth muscle.