Genome-wide association mapping to candidate polymorphism resolution in the unsequenced barley genome

Although commonplace in human disease genetics, genome-wide association (GWA) studies have only relatively recently been applied to plants. Using 32 phenotypes in the inbreeding crop barley, we report GWA mapping of 15 morphological traits across ∼500 cultivars genotyped with 1,536 SNPs. In contrast to the majority of human GWA studies, we observe high levels of linkage disequilibrium within and between chromosomes. Despite this, GWA analysis readily detected common alleles of high penetrance. To investigate the potential of combining GWA mapping with comparative analysis to resolve traits to candidate polymorphism level in unsequenced genomes, we fine-mapped a selected phenotype (anthocyanin pigmentation) within a 140-kb interval containing three genes. Of these, resequencing the putative anthocyanin pathway gene HvbHLH1 identified a deletion resulting in a premature stop codon upstream of the basic helix-loop-helix domain, which was diagnostic for lack of anthocyanin in our association and biparental mapping populations. The methodology described here is transferable to species with limited genomic resources, providing a paradigm for reducing the threshold of map-based cloning in unsequenced crops.

Association mapping of partitioning loci in barley

Background: Association mapping, initially developed in human disease genetics, is now being applied to plant species. The model species Arabidopsis provided some of the first examples of association mapping in plants, identifying previously cloned flowering time genes, despite high population sub-structure. More recently, association genetics has been applied to barley, where breeding activity has resulted in a high degree of population sub-structure. A major genotypic division within barley is that between winter- and spring-sown varieties, which differ in their requirement for vernalization to promote subsequent flowering. To date, all attempts to validate association genetics in barley by identifying major flowering time loci that control vernalization requirement (VRN-H1 and VRN-H2) have failed. Here, we validate the use of association genetics in barley by identifying VRN-H1 and VRN-H2, despite their prominent role in determining population sub-structure. Results: By taking barley as a typical inbreeding crop, and seasonal growth habit as a major partitioning phenotype, we develop an association mapping approach which successfully identifies VRN-H1 and VRN-H2, the underlying loci largely responsible for this agronomic division. We find a combination of Structured Association followed by Genomic Control to correct for population structure and inflation of the test statistic, resolved significant associations only with VRN-H1 and the VRN-H2 candidate genes, as well as two genes closely linked to VRN-H1 (HvCSFs1 and HvPHYC). Conclusion: We show that, after employing appropriate statistical methods to correct for population sub-structure, the genome-wide partitioning effect of allelic status at VRN-H1 and VRN-H2 does not result in the high levels of spurious association expected to occur in highly structured samples. Furthermore, we demonstrate that both VRN-H1 and the candidate VRN-H2 genes can be identified using association mapping. Discrimination between intragenic VRN-H1 markers was achieved, indicating that candidate causative polymorphisms may be discerned and prioritised within a larger set of positive associations. This proof of concept study demonstrates the feasibility of association mapping in barley, even within highly structured populations. A major advantage of this method is that it does not require large numbers of genome-wide markers, and is therefore suitable for fine mapping and candidate gene evaluation, especially in species for which large numbers of genetic markers are either unavailable or too costly.

Construction of a Festuca pratensis BAC library for map-based cloning in Festulolium substitution lines

Introgression in Festulolium is a potentially powerful tool to isolate genes for a large number of traits which differ between Festuca pratensis Huds. and Lolium perenne L. Not only are hybrids between the two species fertile, but the two genomes can be distinguished by genomic in situ hybridisation and a high frequency of recombination occurs between homoeologous chromosomes and chromosome segments. By a programme of introgression and a series of backcrosses, L. perenne lines have been produced which contain small F. pratensis substitutions. This material is a rich source of polymorphic markers targeted towards any trait carried on the F. pratensis substitution not observed in the L. perenne background. We describe here the construction of an F. pratensis BAC library, which establishes the basis of a map-based cloning strategy in L. perenne. The library contains 49,152 clones, with an average insert size of 112 kbp, providing coverage of 2.5 haploid genome equivalents. We have screened the library for eight amplified fragment length polymorphism (AFLP) derived markers known to be linked to an F. pratensis gene introgressed into L. perenne and conferring a staygreen phenotype as a consequence of a mutation in primary chlorophyll catabolism. While for four of the markers it was possible to identify bacterial artificial chromosome (BAC) clones, the other four AFLPs were too repetitive to enable reliable identification of locus-specific BACs. Moreover, when the four BACs were partially sequenced, no obvious coding regions could be identified. This contrasted to BACs identified using cDNA sequences, when multiple genes were identified on the same BAC.

Weaning diet induces sustained metabolic phenotype shift in the pig and influences host response to Bifidobacterium lactis NCC2818

Background: The process of weaning causes a major shift in intestinal microbiota and is a critical period for developing appropriate immune responses in young mammals. Objective: To use a new systems approach to provide an overview of host metabolism and the developing immune system in response to nutritional intervention around the weaning period. Design: Piglets (n¼14) were weaned onto either an eggbased or soya-based diet at 3 weeks until 7 weeks, when all piglets were switched onto a fish-based diet. Half the animals on each weaning diet received Bifidobacterium lactis NCC2818 supplementation from weaning onwards. Immunoglobulin production from immunologically relevant intestinal sites was quantified and the urinary 1H NMR metabolic profile was obtained from each animal at post mortem (11 weeks). Results: Different weaning diets induced divergent and sustained shifts in the metabolic phenotype, which resulted in the alteration of urinary gut microbial co-metabolites, even after 4 weeks of dietary standardisation. B lactis NCC2818 supplementation affected the systemic metabolism of the different weaning diet groups over and above the effects of diet. Additionally, production of gut mucosa-associated IgA and IgM was found to depend upon the weaning diet and on B lactis NCC2818 supplementation. Conclusion: The correlation of urinary 1H NMR metabolic profile with mucosal immunoglobulin production was demonstrated, thus confirming the value of this multiplatform approach in uncovering non-invasive biomarkers of immunity. This has clear potential for translation into human healthcare with the development of urine testing as a means of assessing mucosal immune status. This might lead to early diagnosis of intestinal dysbiosis and with subsequent intervention, arrest disease development. This system enhances our overall understanding of pathologies under supra-organismal control.

Microbial−mammalian cometabolites dominate the age-associated urinary metabolic phenotype in Taiwanese and American populations

Understanding the metabolic processes associated with aging is key to developing effective management and treatment strategies for age-related diseases. We investigated the metabolic profiles associated with age in a Taiwanese and an American population. 1H NMR spectral profiles were generated for urine specimens collected from the Taiwanese Social Environment and Biomarkers of Aging Study (SEBAS; n = 857; age 54–91 years) and the Mid-Life in the USA study (MIDUS II; n = 1148; age 35–86 years). Multivariate and univariate linear projection methods revealed some common age-related characteristics in urinary metabolite profiles in the American and Taiwanese populations, as well as some distinctive features. In both cases, two metabolites—4-cresyl sulfate (4CS) and phenylacetylglutamine (PAG)—were positively associated with age. In addition, creatine and β-hydroxy-β-methylbutyrate (HMB) were negatively correlated with age in both populations (p < 4 × 10–6). These age-associated gradients in creatine and HMB reflect decreasing muscle mass with age. The systematic increase in PAG and 4CS was confirmed using ultraperformance liquid chromatography–mass spectrometry (UPLC–MS). Both are products of concerted microbial–mammalian host cometabolism and indicate an age-related association with the balance of host–microbiome metabolism.

Optimising the analysis of transcript data using high density oligonucleotide arrays and genomic DNA-based probe selection

Background: Affymetrix GeneChip arrays are widely used for transcriptomic studies in a diverse range of species. Each gene is represented on a GeneChip array by a probe- set, consisting of up to 16 probe-pairs. Signal intensities across probe- pairs within a probe-set vary in part due to different physical hybridisation characteristics of individual probes with their target labelled transcripts. We have previously developed a technique to study the transcriptomes of heterologous species based on hybridising genomic DNA (gDNA) to a GeneChip array designed for a different species, and subsequently using only those probes with good homology. Results: Here we have investigated the effects of hybridising homologous species gDNA to study the transcriptomes of species for which the arrays have been designed. Genomic DNA from Arabidopsis thaliana and rice (Oryza sativa) were hybridised to the Affymetrix Arabidopsis ATH1 and Rice Genome GeneChip arrays respectively. Probe selection based on gDNA hybridisation intensity increased the number of genes identified as significantly differentially expressed in two published studies of Arabidopsis development, and optimised the analysis of technical replicates obtained from pooled samples of RNA from rice. Conclusion: This mixed physical and bioinformatics approach can be used to optimise estimates of gene expression when using GeneChip arrays.

Using genomic DNA-based probe-selection to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species

High-density oligonucleotide (oligo) arrays are a powerful tool for transcript profiling. Arrays based on GeneChip® technology are amongst the most widely used, although GeneChip® arrays are currently available for only a small number of plant and animal species. Thus, we have developed a method to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species and tested the method by analysing the transcriptome of Brassica oleracea L., a species for which no GeneChip® array is available, using a GeneChip® array designed for Arabidopsis thaliana (L.) Heynh. Genomic DNA from B. oleracea was labelled and hybridised to the ATH1-121501 GeneChip® array. Arabidopsis thaliana probe-pairs that hybridised to the B. oleracea genomic DNA on the basis of the perfect-match (PM) probe signal were then selected for subsequent B. oleracea transcriptome analysis using a .cel file parser script to generate probe mask files. The transcriptional response of B. oleracea to a mineral nutrient (phosphorus; P) stress was quantified using probe mask files generated for a wide range of gDNA hybridisation intensity thresholds. An example probe mask file generated with a gDNA hybridisation intensity threshold of 400 removed > 68 % of the available PM probes from the analysis but retained >96 % of available A. thaliana probe-sets. Ninety-nine of these genes were then identified as significantly regulated under P stress in B. oleracea, including the homologues of P stress responsive genes in A. thaliana. Increasing the gDNA hybridisation intensity thresholds up to 500 for probe-selection increased the sensitivity of the GeneChip® array to detect regulation of gene expression in B. oleracea under P stress by up to 13-fold. Our open-source software to create probe mask files is freely available http://affymetrix.arabidopsis.info/xspecies/ webcite and may be used to facilitate transcriptomic analyses of a wide range of plant and animal species in the absence of custom arrays.

Metabolic phenotype modulation by caloric restriction in a lifelong dog study

Modeling aging and age-related pathologies presents a substantial analytical challenge given the complexity of gene−environment influences and interactions operating on an individual. A top-down systems approach is used to model the effects of lifelong caloric restriction, which is known to extend life span in several animal models. The metabolic phenotypes of caloric-restricted (CR; n = 24) and pair-housed control-fed (CF; n = 24) Labrador Retriever dogs were investigated by use of orthogonal projection to latent structures discriminant analysis (OPLS-DA) to model both generic and age-specific responses to caloric restriction from the 1H NMR blood serum profiles of young and older dogs. Three aging metabolic phenotypes were resolved: (i) an aging metabolic phenotype independent of diet, characterized by high levels of glutamine, creatinine, methylamine, dimethylamine, trimethylamine N-oxide, and glycerophosphocholine and decreasing levels of glycine, aspartate, creatine and citrate indicative of metabolic changes associated largely with muscle mass; (ii) an aging metabolic phenotype specific to CR dogs that consisted of relatively lower levels of glucose, acetate, choline, and tyrosine and relatively higher serum levels of phosphocholine with increased age in the CR population; (iii) an aging metabolic phenotype specific to CF dogs including lower levels of liproprotein fatty acyl groups and allantoin and relatively higher levels of formate with increased age in the CF population. There was no diet metabotype that consistently differentiated the CF and CR dogs irrespective of age. Glucose consistently discriminated between feeding regimes in dogs (≥312 weeks), being relatively lower in the CR group. However, it was observed that creatine and amino acids (valine, leucine, isoleucine, lysine, and phenylalanine) were lower in the CR dogs (<312 weeks), suggestive of differences in energy source utilization. 1H NMR spectroscopic analysis of longitudinal serum profiles enabled an unbiased evaluation of the metabolic markers modulated by a lifetime of caloric restriction and showed differences in the metabolic phenotype of aging due to caloric restriction, which contributes to longevity studies in caloric-restricted animals. Furthermore, OPLS-DA provided a framework such that significant metabolites relating to life extension could be differentiated and integrated with aging processes.

Prediction of Staphylococcus aureus antimicrobial resistance by whole-genome sequencing

Whole-genome sequencing (WGS) could potentially provide a single platform for extracting all the information required to predict an organism’s phenotype. However, its ability to provide accurate predictions has not yet been demonstrated in large independent studies of specific organisms. In this study, we aimed to develop a genotypic prediction method for antimicrobial susceptibilities. The whole genomes of 501 unrelated Staphylococcus aureus isolates were sequenced, and the assembled genomes were interrogated using BLASTn for a panel of known resistance determinants (chromosomal mutations and genes carried on plasmids). Results were compared with phenotypic susceptibility testing for 12 commonly used antimicrobial agents (penicillin, methicillin, erythromycin, clindamycin, tetracycline, ciprofloxacin, vancomycin, trimethoprim, gentamicin, fusidic acid, rifampin, and mupirocin) performed by the routine clinical laboratory. We investigated discrepancies by repeat susceptibility testing and manual inspection of the sequences and used this information to optimize the resistance determinant panel and BLASTn algorithm. We then tested performance of the optimized tool in an independent validation set of 491 unrelated isolates, with phenotypic results obtained in duplicate by automated broth dilution (BD Phoenix) and disc diffusion. In the validation set, the overall sensitivity and specificity of the genomic prediction method were 0.97 (95% confidence interval [95% CI], 0.95 to 0.98) and 0.99 (95% CI, 0.99 to 1), respectively, compared to standard susceptibility testing methods. The very major error rate was 0.5%, and the major error rate was 0.7%. WGS was as sensitive and specific as routine antimicrobial susceptibility testing methods. WGS is a promising alternative to culture methods for resistance prediction in S. aureus and ultimately other major bacterial pathogens.

Neonatal environment exerts a sustained influence on the development of the intestinal microbiota and metabolic phenotype

The postnatal environment, including factors such as weaning and acquisition of the gut microbiota, has been causally linked to the development of later immunological diseases such as allergy and autoimmunity, and has also been associated with a predisposition to metabolic disorders. We show that the very early-life environment influences the development of both the gut microbiota and host metabolic phenotype in a porcine model of human infants. Farmpiglets were nursed by their mothers for 1 day, before removal to highly controlled, individual isolators where they received formula milk until weaning at 21 days. The experiment was repeated, to create two batches, which differed only in minor environmental fluctuations during the first day. At day 1 after birth, metabolic profiling of serum by 1H nuclear magnetic resonance spectroscopy demonstrated significant, systemic, inter-batch variation which persisted until weaning. However, the urinary metabolic profiles demonstrated that significant inter-batch effects on 3-hydroxyisovalerate, trimethylamine-N-oxide and mannitol persisted beyond weaning to at least 35 days. Batch effects were linked to significant differences in the composition of colonic microbiota at 35 days, determined by 16 S pyrosequencing. Different weaning diets modulated both the microbiota and metabolic phenotype independently of the persistent batch effects. We demonstrate that the environment during the first day of life influences development of the microbiota and metabolic phenotype and thus should be taken into account when interrogating experimental outcomes. In addition, we suggest that intervention at this early time could provide ‘metabolic rescue’ for at-risk infants who have undergone aberrant patterns of initial intestinal colonisation.

Ibuprofen inhibits migration and proliferation of human coronary artery smooth muscle cells by inducing a differentiated phenotype: role of peroxisome proliferator-activated receptor y

Objectives: The search for agents that are capable of preventing restenosis and reduce the risk of late thrombosis is of utmost importance. In this study we aim to evaluate the in vitro effects of ibuprofen on proliferation and migration of human coronary artery smooth muscle cells (HCASMCs) and on human coronary artery endothelial cells (HCAECs) migration. Methods: Cell proliferation was evaluated by direct cell counting using trypan blue exclusion. Cell migration was assessed by wound healing “scratch” assay and by time lapse video-microscopy. Protein expression was assessed by immunoblotting, and morphological changes were studied by immunocytochemistry. The involvement of the PPARγ pathway was studied with the selective agonist troglitazone, and the use of highly selective antagonists of PPARγ such as PGF2α and GW9662. Results: We demonstrate that ibuprofen inhibits proliferation and migration of HCASMCs and induces a switch in HCASMCs towards a differentiated and contractile phenotype, and that these effects are mediated through the PPARγ pathway. Importantly we also show that the effects of ibuprofen are cell type specific as it does not affect migration and proliferation of endothelial cells. Conclusions: Taken together, our results suggest that ibuprofen could be an effective drug for the development of novel drug eluting stents, which could lead reduced rates of restenosis and potentially other complications of DES stent implantation.

Loss of Dictyostelium HSPC300 causes a scar-like phenotype and loss of SCAR protein.

HSPC300 is essential for most SCAR complex functions. The phenotype of HSPC300 knockouts is most similar to mutants in scar, not the other members of the SCAR complex, suggesting that HSPC300 acts most directly on SCAR itself.

Genomic signals of migration and continuity in Britain before the Anglo-Saxons

The purported migrations that have formed the peoples of Britain have been the focus of generations of scholarly controversy. However, this has not benefited from direct analyses of ancient genomes. Here we report nine ancient genomes (~1 x) of individuals from northern Britain: seven from a Roman era York cemetery, bookended by earlier Iron-Age and later Anglo-Saxon burials. Six of the Roman genomes show affinity with modern British Celtic populations, particularly Welsh, but significantly diverge from populations from Yorkshire and other eastern English samples. They also show similarity with the earlier Iron-Age genome, suggesting population continuity, but differ from the later Anglo-Saxon genome. This pattern concords with profound impact of migrations in the Anglo-Saxon period. Strikingly, one Roman skeleton shows a clear signal of exogenous origin, with affinities pointing towards the Middle East, confirming the cosmopolitan character of the Empire, even at its northernmost fringes.

Acute benefits of the microbial-derived isoflavone metabolite equol on arterial stiffness in men prospectively recruited according to equol producer phenotype: a double-blind randomized controlled trial

There is much speculation with regard to the potential cardioprotective benefits of equol, a microbial-derived metabolite of the isoflavone daidzein, which is produced in the large intestine after soy intake in 30% of Western populations. Although cross-sectional and retrospective data support favorable associations between the equol producer (EP) phenotype and cardiometabolic health, few studies have prospectively recruited EPs to confirm this association. The aim was to determine whether the acute vascular benefits of isoflavones differ according to EP phenotype and subsequently investigate the effect of providing commercially produced S-(–)equol to non-EPs. We prospectively recruited male EPs and non-EPs (n = 14/ group) at moderate cardiovascular risk into a double-blind, placebocontrolled crossover study to examine the acute effects of soy isoflavones (80-mg aglycone equivalents) on arterial stiffness [carotid-femoral pulse-wave velocity (cfPWV)], blood pressure, endothelial function (measured by using the EndoPAT 2000; Itamar Medical), and nitric oxide at baseline (0 h) and 6 and 24 h after intake. In a separate assessment, non-EPs consumed 40 mg S-(–)equol with identical vascular measurements performed 2 h after intake. After soy intake, cfPWV significantly improved in EPs at 24 h (cfPWV change from 0 h: isoflavone, 20.2 6 0.2 m/s; placebo, 0.6 6 0.2 m/s; P , 0.01), which was significantly associated with plasma equol concentrations (R = 20.36, P = 0.01). No vascular effects were observed in EPs at 6 h or in non-EPs at any time point. Similarly, no benefit of commercially produced S-(–)equol was observed in non-EPs despite mean plasma equol concentrations reaching 3.2 mmol/L. Acute soy intake improved cfPWV in EPs, equating to an 11–12% reduced risk of cardiovascular disease if sustained. However, a single dose of commercially produced equol had no cardiovascular benefits in non-EPs. These data suggest that the EP phenotype is critical in unlocking the vascular benefits of equol in men, and long-term trials should focus on confirming the implications of EP phenotype on cardiovascular health. This trial was registered at clinicaltrials.gov as NCT01530893. Am J Clin Nutr doi: 10.3945/ajcn.115.125690.

Genomic evidence for the Pleistocene and recent population history of Native Americans

How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we found that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (ka) and after no more than an 8000-year isolation period in Beringia. After their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 ka, one that is now dispersed across North and South America and the other restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative “Paleoamerican” relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model.