Unexpected Relationships and Inbreeding in HapMap Phase III Populations

Correct annotation of the genetic relationships between samples is essential for population genomic studies, which could be biased by errors or omissions. To this end, we used identity-by-state (IBS) and identity-by-descent (IBD) methods to assess genetic relatedness of individuals within HapMap phase III data. We analyzed data from 1,397 individuals across 11 ethnic populations. Our results support previous studies (Pemberton et al., 2010; Kyriazopoulou-Panagiotopoulou et al., 2011) assessing unknown relatedness present within this population. Additionally, we present evidence for 1,657 novel pairwise relationships across 9 populations. Surprisingly, significant Cotterman''s coefficients of relatedness K1 (IBD1) values were detected between pairs of known parents. Furthermore, significant K2 (IBD2) values were detected in 32 previously annotated parent-child relationships. Consistent with a hypothesis of inbreeding, regions of homozygosity (ROH) were identified in the offspring of related parents, of which a subset overlapped those reported in previous studies (Gibson et al. 2010; Johnson et al. 2011). In total, we inferred 28 inbred individuals with ROH that overlapped areas of relatedness between the parents and/or IBD2 sharing at a different genomic locus between a child and a parent. Finally, 8 previously annotated parent-child relationships had unexpected K0 (IBD0) values (resulting from a chromosomal abnormality or genotype error), and 10 previously annotated second-degree relationships along with 38 other novel pairwise relationships had unexpected IBD2 (indicating two separate paths of recent ancestry). These newly described types of relatedness may impact the outcome of previous studies and should inform the design of future studies relying on the HapMap Phase III resource.     

Identification of mitogen-activated protein/extracellular signal-responsive kinase kinase 2 as a novel partner of the scaffolding protein human homolog of disc-large

Human disc-large homolog (hDlg), also known as synapse-associated protein 97, is a scaffold protein, a member of the membrane-associated guanylate kinase family, implicated in neuronal synapses and epithelial-epithelial cell junctions whose expression and function remains poorly characterized in most tissues, particularly in the vasculature. In human vascular tissues, hDlg is highly expressed in smooth muscle cells (VSMCs). Using the yeast two-hybrid system to screen a human aorta cDNA library, we identified mitogen-activated protein/extracellular signal-responsive kinase (ERK) kinase (MEK)2, a member of the ERK cascade, as an hDlg binding partner. Site-directed mutagenesis showed a major involvement of the PSD-95, disc-large, ZO-1 domain-2 of hDlg and the C-terminal sequence RTAV of MEK2 in this interaction. Coimmunoprecipitation assays in both human VSMCs and human embryonic kidney 293 cells, demonstrated that endogenous hDlg physically interacts with MEK2 but not with MEK1. Confocal microscopy suggested a colocalization of the two proteins at the inner layer of the plasma membrane of confluent human embryonic kidney 293 cells, and in a perinuclear area in human VSMCs. Additionally, hDlg also associates with the endoplasmic reticulum and microtubules in these latter cells. Taken together, these findings allow us to hypothesize that hDlg acts as a MEK2-specific scaffold protein for the ERK signaling pathway, and may improve our understanding of how scaffold proteins, such as hDlg, differentially tune MEK1/MEK2 signaling and cell responses.     

Engineering glycoproteins for use as pharmaceuticals

The fact that glycosylation is not a significant process in prokaryotes means that many of the proteins produced by genetically engineered bacteria are not identical to their eukaryotic counterparts. Although glycosylation affects the physical, chemical and biological nature of proteins, its pharmacological value in potential protein pharmaceuticals is not easy to predict. However, the development of mammalian cell culture methods for expressing recombinant DNA-derived glycoproteins will permit further studies in the field.     

Extent and diversity of human alternative splicing established by complementary database annotation and microarray analysis

Alternative splicing generates functional diversity in higher organisms through alternative first and last exons, skipped and included exons, intron retentions and alternative donor, and acceptor sites. In large-scale microarray studies in humans and the mouse, emphasis so far has been placed on exon-skip events, leaving the prevalence and importance of other splice types largely unexplored. Using a new human splice variant database and a genome-wide microarray to probes thousands of splice events of each type, we measured differential expression of splice types across six pair of diverse cell lines and validated the database annotation process. Results suggest that splicing in humans is more complex than simple exon-skip events, which account for a minority of splicing differences. The relative frequency of differential expression of the splice types correlates with what is found by our annotation efforts. In conclusion, alternative splicing in human cells is considerably more complex than the canonical example of the exon skip. The complementary approaches of genome-wide annotation of alternative splicing in human and design of genome-wide splicing microarrays to measure differential splicing in biological samples provide a powerful high-throughput tool to study the role of alternative splicing in human biology.     

Identification of candidate genes for dyslexia susceptibility on chromosome 18

Background

Six independent studies have identified linkage to chromosome 18 for developmental dyslexia or general reading ability. Until now, no candidate genes have been identified to explain this linkage. Here, we set out to identify the gene(s) conferring susceptibility by a two stage strategy of linkage and association analysis.

Methodology/Principal Findings

Linkage analysis: 264 UK families and 155 US families each containing at least one child diagnosed with dyslexia were genotyped with a dense set of microsatellite markers on chromosome 18. Association analysis: Using a discovery sample of 187 UK families, nearly 3000 SNPs were genotyped across the chromosome 18 dyslexia susceptibility candidate region. Following association analysis, the top ranking SNPs were then genotyped in the remaining samples. The linkage analysis revealed a broad signal that spans approximately 40 Mb from 18p11.2 to 18q12.2. Following the association analysis and subsequent replication attempts, we observed consistent association with the same SNPs in three genes; melanocortin 5 receptor (MC5R), dymeclin (DYM) and neural precursor cell expressed, developmentally down-regulated 4-like (NEDD4L).

Conclusions

Along with already published biological evidence, MC5R, DYM and NEDD4L make attractive candidates for dyslexia susceptibility genes. However, further replication and functional studies are still required.     

Site-specific chlorophyll reference conditions for lakes in Northern and Western Europe

The Water Framework Directive (WFD) requires EU Member States to assess the “ecological status” of surface waters. As a component of ecological status, many European countries are developing a classification scheme for chlorophyll concentrations as a measure of phytoplankton biomass. The chlorophyll classification must be based on the degree of divergence of a water body from an appropriate baseline or ‘reference condition’. This article describes the development of a series of regression models for predicting reference chlorophyll concentrations on a site-specific basis. For model development, a large dataset of European lakes considered to be in reference condition, 466 lakes in total, was assembled. Data were included from 12 European countries, but lakes from Northern and Western Europe dominated and made up 92% of all reference lakes. Data have been collated on chlorophyll concentration, altitude, mean depth, alkalinity, humic type, surface area and geographical region. Regression models were developed for estimating site-specific reference chlorophyll concentrations from significant predictor ‘typology’ variables. Reference chlorophyll concentrations were found to vary along a number of environmental gradients. Concentrations increased with colour and alkalinity and decreased with lake depth and altitude. Forward selection was used to identify independent explanatory variables in regression models for predicting site-specific reference chlorophyll concentrations. Depth was selected as an explanatory variable in all models. Alkalinity was included in models for low colour and humic lakes and altitude was included in models for low colour and very humic lakes. Uncertainty in the models was quite high and arises from errors in the data used to develop the models (including natural temporal and spatial variability in data) and also from additional explanatory variables not considered in the models, particularly nutrient concentrations, retention time and grazing. Despite these uncertainties, site-specific reference conditions are still recommended in preference to type-specific reference conditions, as they use the individual characteristics of a site known to influence phytoplankton biomass, rather than adopt standards set to generally represent a large population of lakes of a particular type. For this reason, site-specific reference conditions should result in reduced error in ecological status classifications, particularly for lakes close to typology boundaries.     

Increased Expression of a Phloem Membrane Protein Encoded by NHL26 Alters Phloem Export and Sugar Partitioning in Arabidopsis

The complex process of phloem sugar transport involves symplasmic and apoplasmic events. We characterized Arabidopsis thaliana lines ectopically expressing a phloem-specific gene encoding NDR1/HIN1-like26 (NHL26), a putative membrane protein. NHL26 overexpressor plants grew more slowly than wild-type plants, accumulated high levels of carbohydrates in mature leaves, and had a higher shoot biomass, contrasting with slower root growth and a lower seed yield. Similar effects were observed when NHL26 was overexpressed in companion cells, under the control of a companion cell–specific promoter. The soluble sugar content of the phloem sap and sink organs was lower than that in the wild type, providing evidence of a sugar export defect. This was confirmed in a phloem-export assay with the symplastic tracer carboxyfluorescein diacetate. Leaf sugar accumulation was accompanied by higher organic acid, amino acid, and protein contents, whereas analysis of the metabolite profile of phloem sap exudate revealed no change in amino acid or organic acid content, indicating a specific effect on sugar export. NHL26 was found to be located in the phloem plasmodesmata and the endoplasmic reticulum. These findings reveal that NHL26 accumulation affects either the permeability of plasmodesmata or sugar signaling in companion cells, with a specific effect on sugar export.     

Glycosylated cell-penetrating peptides and their conjugates to a proapoptotic peptide: preparation by click chemistry and cell viability studies

Cell-penetrating peptides (CPPs), which are usually short basic peptides, are able to cross cell membranes and convey bioactive cargoes inside cells. CPPs have been widely used to deliver inside cells peptides, proteins, and oligonucleotides; however, their entry mechanisms still remain controversial. A major problem concerning CPPs remains their lack of selectivity to target a specific type of cell and/or an intracellular component. We have previously shown that myristoylation of one of these CPPs affected the intracellular distribution of the cargo. We report here on the synthesis of glycosylated analogs of the cell-penetrating peptide (R6/W3): Ac-RRWWRRWRR-NH₂. One, two, or three galactose(s), with or without a spacer, were introduced into the sequence of this nonapeptide via a triazole link, the Huisgen reaction being achieved on a solid support. Four of these glycosylated CPPs were coupled via a disulfide bridge to the proapoptotic KLAK peptide, (KLAKLAKKLAKLAK), which alone does not enter into cells. The effect on cell viability and the uptake efficiency of different glycosylated conjugates were studied on CHO cells and were compared to those of the nonglycosylated conjugates: (R6/W3)S-S-KLAK and penetratinS-S-KLAK. We show that glycosylation significantly increases the cell viability of CHO cells compared to the nonglycosylated conjugates and concomitantly decreases the internalization of the KLAK cargo. These results suggest that glycosylation of CPP may be a key point in targeting specific cells.     

Insights into interactions between the alpha-helical region of the salmon calcitonin antagonists and the human calcitonin receptor using photoaffinity labeling

Fish-like calcitonins (CTs), such as salmon CT (sCT), are widely used clinically in the treatment of bone-related disorders; however, the molecular basis for CT binding to its receptor, a class II G protein-coupled receptor, is not well defined. In this study we have used photoaffinity labeling to identify proximity sites between CT and its receptor. Two analogues of the antagonist sCT(8-32) containing a single photolabile p-benzoyl-l-phenylalanine (Bpa) residue in position 8 or 19 were used. Both analogues retained high affinity for the CT receptor and potently inhibited agonist-induced cAMP production. The [Bpa(19)]sCT(8-32) analogue cross-linked to the receptor at or near the equivalent cross-linking site of the full-length peptide, within the fragment Cys(134)-Lys(141) (within the amino terminus of the receptor, adjacent to transmembrane 1) (Pham, V., Wade, J. D., Purdue, B. W., and Sexton, P. M. (2004) J. Biol. Chem. 279, 6720-6729). In contrast, proteolytic mapping and mutational analysis identified Met(49) as the cross-linking site for [Bpa(8)]sCT(8-32). This site differed from the previously identified cross-linking site of the agonist [Bpa(8)]human CT (Dong, M., Pinon, D. I., Cox, R. F., and Miller, L. J. (2004) J. Biol. Chem. 279, 31177-31182) and may provide evidence for conformational differences between interaction with active and inactive state receptors. Molecular modeling suggests that the difference in cross-linking between the two Bpa(8) analogues can be accounted for by a relatively small change in peptide orientation. The model was also consistent with cooperative interaction between the receptor amino terminus and the receptor core.