Simultaneously accounting for population structure, genotype by environment interaction, and spatial variation in marker-trait associations in sugarcane

Few association mapping studies have simultaneously accounted for population structure, genotype by environment interaction (GEI), and spatial variation. In this sugarcane association mapping study we tested models accounting for these factors and identified the impact that each model component had on the list of markers declared as being significantly associated with traits. About 480 genotypes were evaluated for cane yield and sugar content at three sites and scored with DArT markers. A mixed model was applied in analysis of the data to simultaneously account for the impacts of population structure, GEI, and spatial variation within a trial. Two forms of the DArT marker data were used in the analysis: the standard discrete data (0, 1) and a continuous DArT score, which is related to the marker dosage. A large number of markers were significantly associated with cane yield and sugar content. However, failure to account for population structure, GEI, and (or) spatial variation produced both type I and type II errors, which on the one hand substantially inflated the number of significant markers identified (especially true for failing to account for GEI) and on the other hand resulted in failure to detect markers that could be associated with cane yield or sugar content (especially when failing to account for population structure). We concluded that association mapping based on trials from one site or analysis that failed to account for GEI would produce many trial-specific associated markers that would have low value in breeding programs.     

Clustered O-Glycans of IgA1: DEFINING MACRO- AND MICROHETEROGENEITY BY USE OF ELECTRON CAPTURE/TRANSFER DISSOCIATION*

IgA nephropathy (IgAN) is the most common primary glomerulonephritis in the world. Aberrantly glycosylated IgA1, with galactose (Gal)-deficient hinge region (HR) O-glycans, plays a pivotal role in the pathogenesis of the disease. It is not known whether the glycosylation defect occurs randomly or preferentially at specific sites. We have described the utility of activated ion-electron capture dissociation (AI-ECD) mass spectrometric analysis of IgA1 O-glycosylation. However, locating and characterizing the entire range of O-glycan attachment sites are analytically challenging due to the clustered serine and threonine residues in the HR of IgA1 heavy chain. To address this problem, we analyzed all glycoforms of the HR glycopeptides of a Gal-deficient IgA1 myeloma protein, mimicking the aberrant IgA1 in patients with IgAN, by use of a combination of IgA-specific proteases + trypsin and AI-ECD Fourier transform ion cyclotron resonance (FT-ICR) tandem mass spectrometry (MS/MS). The IgA-specific proteases provided a variety of IgA1 HR fragments that allowed unambiguous localization of all O-glycosylation sites in the six most abundant glycoforms, including the sites deficient in Gal. Additionally, this protocol was adapted for on-line liquid chromatography (LC)-AI-ECD MS/MS and LC-electron transfer dissociation MS/MS analysis. Our results thus represent a new clinically relevant approach that requires ECD/electron transfer dissociation-type fragmentation to define the molecular events leading to pathogenesis of a chronic kidney disease. Furthermore, this work offers generally applicable principles for the analysis of clustered sites of O-glycosylation.Glycosylation is one of the most common post-translational modifications of proteins. It is estimated that over half of mammalian proteins are glycosylated. Patients with several autoimmune disorders, chronic inflammatory diseases, and some infectious diseases exhibit abnormal glycosylation of serum immunoglobulins and other glycoproteins (1–5). The biological functions of these modifications in health and disease have become a significant area of interest in biomedical research (6). A subset of these glycoproteins has clustered sites of O-glycosylation with serine- and threonine-rich stretches within the amino acid sequence. Mucins, such as membrane-associated MUC1, are perhaps the best known family of proteins that are heavily O-glycosylated. Their altered expression and aberrant glycosylation have made them potential targets as biomarkers for early detection of cancer (7). Immunoglobulin A1 (IgA1)1 contains both O- and N-glycans (Fig. 1). Aberrant O-glycosylation of IgA1 is involved in the pathogenesis of IgA nephropathy (IgAN) and the closely related Henoch-Schönlein purpura nephritis (1, 8). Interestingly, the aberrantly glycosylated molecules, IgA1 in IgAN and MUC1 in cancer, are recognized by the immune system as neoepitopes as evidenced by formation of specific antibodies (9–11). Mucin-like bacterial surface proteins exhibit similar properties: the molecules have clustered bacterial O-glycans that mediate cellular adhesion, and blocking antibodies target these glycan-containing epitopes (12).Open in a separate windowFig. 1.IgA1 structural elements. IgA1 has N-linked glycans (filled circles) and O-linked glycans (open circles). The O-glycosylated sites are in the HR between the first and second constant region domains of the heavy chains. The HR is a Pro-rich segment with nine possible sites of O-glycan attachment. Underlined serine and threonine residues are usually glycosylated (31). Arrows show cleavage sites of trypsin and IgA-specific proteases.An O-glycosylated protein from a single source contains a population of variably O-glycosylated isoforms that show a distinct distribution of microheterogeneity of the O-glycan chains in terms of number, sites of attachment, and composition. Characterizing these clustered sites and understanding how the distributions change under different biological conditions or disease states are an analytical challenge. Enzymatic or chemical release of O-glycans is not selective. The heterogeneity, composition, and quantitative aspects of different O-glycan chains can be assessed and quantified by gas chromatographic and/or mass spectrometric techniques. However, the site-specific information and context of location and composition of adjacent chains are lost. Carbohydrate-specific lectin analysis of O-glycoproteins can provide information on glycan composition and comparative differences between samples, such as those from healthy controls and patients with various disease states. We have successfully demonstrated this in the analysis of IgA1 O-glycans from patients with IgAN versus healthy controls and disease controls (13–15). This included proximal assessment of sites with galactose (Gal)-deficient O-glycans after digests with IgA-specific proteases (8). Several studies have demonstrated the value of mass spectrometry (MS) in identifying Gal-deficient IgA1 in patients with IgAN (16–21), including our work that demonstrated the first direct localization of native sites of O-glycan chains in the hinge region (HR) of IgA1 by use of electron capture dissociation (ECD) (20, 22). ECD and the more recently developed electron transfer dissociation (ETD) have been used to identify sites of O-glycosylation on a variety of proteins (23–26). This includes the analysis of sites of O-glycosylation by on-line LC-ECD/ETD MS/MS methods (23, 26, 27).IgAN is the most common primary glomerulonephritis worldwide (28) with about 20–40% of patients developing end stage renal failure. It is characterized by mesangial deposits of IgA1-containing immune complexes (28). The distinctive O-glycan chains of IgA1 molecules play a pivotal role in the pathogenesis of IgAN (1, 10, 14–16, 29, 30). IgA1 contains an HR between the first and second heavy chain constant region domains with a high content of Ser, Thr, and Pro. This segment usually has three to five O-glycan chains per HR (31) (see Fig. 1). Aberrantly glycosylated IgA1, deficient in Gal in some of the O-glycans in the HR, in serum is rare in healthy individuals but is present at elevated levels in IgAN patients (13, 15). This distinctive IgA1 is in circulating immune complexes (8, 10, 15) and in the glomerular deposits of IgAN patients (16, 29). The absence of Gal apparently leads to the exposure of neoepitopes, including terminal and sialylated N-acetylgalactosamine (GalNAc) residues (9, 10). These epitopes are recognized by naturally occurring anti-glycan IgG or IgA1 antibodies and, consequently, circulating immune complexes are formed (9, 10, 15) that can deposit in the glomerular mesangia. To identify the pathogenic forms of IgA1, a thorough analysis of O-glycan microheterogeneity, including identification of the attachment sites, will be required.In this work, we demonstrate the complete analysis of O-glycoform microheterogeneity and site localization of the glycoforms in a naturally Gal-deficient IgA1 (Ale) myeloma protein that mimics the nephritogenic IgA1 in patients with IgAN (8, 9). Reversed phase (RP) LC FT-ICR MS successfully identified 10 distinct IgA1 HR fragments representing >99% of total IgA1. AI-ECD of the six most abundant IgA1 HR glycoforms (>95% of total IgA1) was accomplished with three distinct IgA-specific protease + trypsin digestions, identifying sites of Gal deficiency across four distinct IgA1 O-glycoforms. Based on the success of the ECD fragmentation of these IgA1 HR fragments, we adapted the analysis for on-line LC-MS/MS methods for both ECD and ETD. The variety of IgA1 HR proteolytic fragments provides a practical set of guidelines for the ECD/ETD analysis of clustered sites of O-glycosylation on this and other proteins. These results also provide insight into the order of attachment of the O-glycans in the IgA1 HR.     

Independent wheat B and G genome origins in outcrossing Aegilops progenitor haplotypes

The origin of modern wheats involved alloploidization among related genomes. To determine if Aegilops speltoides was the donor of the B and G genomes in AABB and AAGG tetraploids, we used a 3-tiered approach. Using 70 amplified fragment length polymorphism (AFLP) loci, we sampled molecular diversity among 480 wheat lines from their natural habitats encompassing all S genome Aegilops, the putative progenitors of wheat B and G genomes. Fifty-nine Aegilops representatives for S genome diversity were compared at 375 AFLP loci with diploid, tetraploid, and 11 nulli-tetrasomic Triticum aestivum wheat lines. B genome-specific markers allowed pinning the origin of the B genome to S chromosomes of A. speltoides, while excluding other lineages. The outbreeding nature of A. speltoides influences its molecular diversity and bears upon inferences of B and G genome origins. Haplotypes at nuclear and chloroplast loci ACC1, G6PDH, GPT, PGK1, Q, VRN1, and ndhF for approximately 70 Aegilops and Triticum lines (0.73 Mb sequenced) reveal both B and G genomes of polyploid wheats as unique samples of A. speltoides haplotype diversity. These have been sequestered by the AABB Triticum dicoccoides and AAGG Triticum araraticum lineages during their independent origins.     

Shedding light on host niches: label‐free in situ detection of Mycobacterium gordonae via carotenoids in macrophages by Raman microspectroscopy

Macrophages are the primary habitat of pathogenic mycobacteria during infections. Current research about the host–pathogen interaction on the cellular level is still going on. The present study proves the potential of Raman microspectroscopy as a label‐free and non‐invasive method to investigate intracellular mycobacteria in situ. Therefore, macrophages were infected with Mycobacterium gordonae, a mycobacterium known to cause inflammation linked to intracellular survival in macrophages. Here, we show that Raman maps provided spatial and spectral information about the position of bacteria within determined cell margins of macrophages in two‐dimensional scans and in three‐dimensional image stacks. Simultaneously, the relative intracellular concentration and distributions of cellular constituents such as DNA, proteins and lipids provided phenotypic information about the infected macrophages. Locations of bacteria outside or close to the outer membrane of the macrophages were notably different in their spectral pattern compared with intracellular once. Furthermore, accumulations of bacteria inside of macrophages exhibit distinct spectral/molecular information because of the chemical composition of the intracellular microenvironment. The data show that the connection of microscopically and chemically gained information provided by Raman microspectroscopy offers a new analytical way to detect and to characterize the mycobacterial infection of macrophages.     

Modelling the genetic architecture of flowering time control in barley through nested association mapping

Background

Barley, globally the fourth most important cereal, provides food and beverages for humans and feed for animal husbandry. Maximizing grain yield under varying climate conditions largely depends on the optimal timing of flowering. Therefore, regulation of flowering time is of extraordinary importance to meet future food and feed demands. We developed the first barley nested association mapping (NAM) population, HEB-25, by crossing 25 wild barleys with one elite barley cultivar, and used it to dissect the genetic architecture of flowering time.

Results

Upon cultivation of 1,420 lines in multi-field trials and applying a genome-wide association study, eight major quantitative trait loci (QTL) were identified as main determinants to control flowering time in barley. These QTL accounted for 64% of the cross-validated proportion of explained genotypic variance (p_G). The strongest single QTL effect corresponded to the known photoperiod response gene Ppd-H1. After sequencing the causative part of Ppd-H1, we differentiated twelve haplotypes in HEB-25, whereof the strongest exotic haplotype accelerated flowering time by 11 days compared to the elite barley haplotype. Applying a whole genome prediction model including main effects and epistatic interactions allowed predicting flowering time with an unmatched accuracy of 77% of cross-validated p_G.

Conclusions

The elaborated causal models represent a fundamental step to explain flowering time in barley. In addition, our study confirms that the exotic biodiversity present in HEB-25 is a valuable toolbox to dissect the genetic architecture of important agronomic traits and to replenish the elite barley breeding pool with favorable, trait-improving exotic alleles.

     

Genetic map of triticale compiling DArT, SSR, and AFLP markers

A set of 90 doubled haploid (DH) lines derived from F(1) plants that originated from a cross between × Triticosecale Wittm. 'Saka3006' and ×Triticosecale Wittm. 'Modus', via wide crossing with maize, were used to create a genetic linkage map of triticale. The map has 21 linkage groups assigned to the A, B, and R genomes including 155 simple sequence repeat (SSR), 1385 diversity array technology (DArT), and 28 amplified fragment length polymorphism (AFLP) markers covering 2397 cM with a mean distance between two markers of 4.1 cM. Comparative analysis with wheat consensus maps revealed that triticale chromosomes of the A and B genomes were represented by 15 chromosomes, including combinations of 2AS.2AL#, 2AL#2BL, 6AS.6AL#, and 2BS.6AL# instead of 2A, 2B, and 6A. In respect to published maps of rye, substantial rearrangements were found also for chromosomes 1R, 2R, and 3R of the rye genome. Chromosomes 1R and 2R were truncated and the latter was linked with 3R. A nonhomogeneous distribution of markers across the triticale genome was observed with evident bias (48%) towards the rye genome. This genetic map may serve as a reference linkage map of triticale for efficient studies of structural rearrangements, gene mapping, and marker-assisted selection.     

Novel SSR markers from BAC-end sequences, DArT arrays and a comprehensive genetic map with 1,291 marker loci for chickpea (Cicer arietinum L.)

Chickpea (Cicer arietinum L.) is the third most important cool season food legume, cultivated in arid and semi-arid regions of the world. The goal of this study was to develop novel molecular markers such as microsatellite or simple sequence repeat (SSR) markers from bacterial artificial chromosome (BAC)-end sequences (BESs) and diversity arrays technology (DArT) markers, and to construct a high-density genetic map based on recombinant inbred line (RIL) population ICC 4958 (C. arietinum)×PI 489777 (C. reticulatum). A BAC-library comprising 55,680 clones was constructed and 46,270 BESs were generated. Mining of these BESs provided 6,845 SSRs, and primer pairs were designed for 1,344 SSRs. In parallel, DArT arrays with ca. 15,000 clones were developed, and 5,397 clones were found polymorphic among 94 genotypes tested. Screening of newly developed BES-SSR markers and DArT arrays on the parental genotypes of the RIL mapping population showed polymorphism with 253 BES-SSR markers and 675 DArT markers. Segregation data obtained for these polymorphic markers and 494 markers data compiled from published reports or collaborators were used for constructing the genetic map. As a result, a comprehensive genetic map comprising 1,291 markers on eight linkage groups (LGs) spanning a total of 845.56 cM distance was developed (http://cmap.icrisat.ac.in/cmap/sm/cp/thudi/). The number of markers per linkage group ranged from 68 (LG 8) to 218 (LG 3) with an average inter-marker distance of 0.65 cM. While the developed resource of molecular markers will be useful for genetic diversity, genetic mapping and molecular breeding applications, the comprehensive genetic map with integrated BES-SSR markers will facilitate its anchoring to the physical map (under construction) to accelerate map-based cloning of genes in chickpea and comparative genome evolution studies in legumes.     

DArT for high-throughput genotyping of Cassava (Manihot esculenta) and its wild relatives

Understanding the distribution of genetic diversity within and among individuals, populations, species and gene pools is crucial for the efficient management of germplasm collections. Molecular markers are playing an increasing role in germplasm characterization, yet their broad application is limited by the availability of markers, the costs and the low throughput of existing technologies. This is particularly true for crops of resource-poor farmers such as cassava, Manihot esculenta. Here we report on the development of Diversity Arrays Technology (DArT) for cassava. DArT uses microarrays to detect DNA polymorphism at several hundred genomic loci in a single assay without relying on DNA sequence information. We tested three complexity reduction methods and selected the two that generated genomic representations with the largest frequency of polymorphic clones (PstI/TaqI: 14.6%, PstI/BstNI: 17.2%) to produce large genotyping arrays. Nearly 1,000 candidate polymorphic clones were detected on the two arrays. The performance of the PstI/TaqI array was validated by typing a group of 38 accessions, 24 of them in duplicate. The average call rate was 98.1%, and the scoring reproducibility was 99.8%. DArT markers displayed fairly high polymorphism information content (PIC) values and revealed genetic relationships among the samples consistent with the information available on these samples. Our study suggests that DArT offers advantages over current technologies in terms of cost and speed of marker discovery and analysis. It can therefore be used to genotype large germplasm collections.Electronic Supplementary Material Supplementary material is available for this article at     

16S rRNA gene-based oligonucleotide microarray for environmental monitoring of the betaproteobacterial order "Rhodocyclales"

For simultaneous identification of members of the betaproteobacterial order "Rhodocyclales" in environmental samples, a 16S rRNA gene-targeted oligonucleotide microarray (RHC-PhyloChip) consisting of 79 probes was developed. Probe design was based on phylogenetic analysis of available 16S rRNA sequences from all cultured and as yet uncultured members of the "Rhodocyclales." The multiple nested probe set was evaluated for microarray hybridization with 16S rRNA gene PCR amplicons from 29 reference organisms. Subsequently, the RHC-PhyloChip was successfully used for cultivation-independent "Rhodocyclales" diversity analysis in activated sludge from an industrial wastewater treatment plant. The implementation of a newly designed "Rhodocyclales"-selective PCR amplification system prior to microarray hybridization greatly enhanced the sensitivity of the RHC-PhyloChip and thus enabled the detection of "Rhodocyclales" populations with relative abundances of less than 1% of all bacteria (as determined by fluorescence in situ hybridization) in the activated sludge. The presence of as yet uncultured Zoogloea-, Ferribacterium/Dechloromonas-, and Sterolibacterium-related bacteria in the industrial activated sludge, as indicated by the RHC-PhyloChip analysis, was confirmed by retrieval of their 16S rRNA gene sequences and subsequent phylogenetic analysis, demonstrating the suitability of the RHC-PhyloChip as a novel monitoring tool for environmental microbiology.