Moment estimators of relatedness from low-depth whole-genome sequencing data.

BACKGROUND: Estimating relatedness is an important step for many genetic study designs. A variety of methods for estimating coefficients of pairwise relatedness from genotype data have been proposed. Both the kinship coefficient [Formula: see text] and the fraternity coefficient [Formula: see text] for all pairs of individuals are of interest. However, when dealing with low-depth sequencing or imputation data, individual level genotypes cannot be confidently called. To ignore such uncertainty is known to result in biased estimates. Accordingly, methods have recently been developed to estimate kinship from uncertain genotypes. RESULTS: We present new method-of-moment estimators of both the coefficients [Formula: see text] and [Formula: see text] calculated directly from genotype likelihoods. We have simulated low-depth genetic data for a sample of individuals with extensive relatedness by using the complex pedigree of the known genetic isolates of Cilento in South Italy. Through this simulation, we explore the behaviour of our estimators, demonstrate their properties, and show advantages over alternative methods. A demonstration of our method is given for a sample of 150 French individuals with down-sampled sequencing data. CONCLUSIONS: We find that our method can provide accurate relatedness estimates whilst holding advantages over existing methods in terms of robustness, independence from external software, and required computation time. The method presented in this paper is referred to as LowKi (Low-depth Kinship) and has been made available in an R package ( https://github.com/genostats/LowKi ).

Whole-Exome Sequencing in the Isolated Populations of Cilento from South Italy.

The present study describes the genetic architecture of the isolated populations of Cilento, through the analysis of exome sequence data of 245 representative individuals of these populations. By annotating the exome variants and cataloguing them according to their frequency and functional effects, we identified 347,684 variants, 67.4% of which are rare and low frequency variants, and 1% of them (corresponding to 319 variants per person) are classified as high functional impact variants; also, 39,946 (11.5% of the total) are novel variants, for which we determined a significant enrichment for deleterious effects. By comparing the allele frequencies in Cilento with those from the Tuscan population from the 1000 Genomes Project Phase 3, we highlighted an increase in allele frequency in Cilento especially for variants which map to genes involved in extracellular matrix formation and organization. Furthermore, among the variants showing increased frequency we identified several known rare disease-causing variants. By different population genetics analyses, we corroborated the status of the Cilento populations as genetic isolates. Finally, we showed that exome data of Cilento represents a useful local reference panel capable of improving the accuracy of genetic imputation, thus adding power to genetic studies of human traits in these populations.

Survival of Mycobacterium avium subsp. paratuberculosis in the intermediate and final digestion products of biogas plants.

AIMS: To evaluate the survival of Mycobacterium avium subsp. paratuberculosis (MAP) during anaerobic digestion (AD), we studied two different biogas plants loaded with manure and slurry from paratuberculosis-infected dairy herds. METHODS AND RESULTS: Both plants were operating under mesophilic conditions, the first with a single digester and the second with a double digester. Mycobacterium avium subsp. paratuberculosis detection was performed by sampling each stage of the process, specifically the prefermenter, fermenter, liquid digestate and solid digestate stages, for 11 months. In both plants, MAP was isolated from the prefermenter stage. Only the final products, the solid and liquid digestates, of the one-stage plant showed viable MAP, while no viable MAP was detected in the digestates of the two-stage plant. CONCLUSIONS: Mycobacterium avium subsp. paratuberculosis showed a significant decrease during subsequent steps of the AD process, particularly in the two-stage plant. We suggest that the second digester maintained the digestate under anaerobic conditions for a longer period of time, thus reducing MAP survival and MAP load under the culture detection limit. SIGNIFICANCE AND IMPACT OF THE STUDY: Our data are unable to exclude the presence of MAP in the final products of the biogas plants, particularly those products from the single digester; therefore, the use of digestates as fertilizers is a real concern related to the possible environmental contamination with MAP.

Body composition, leg length and blood pressure in a rural Italian population: a test of the capacity-load model.

BACKGROUND AND AIMS: Whereas adult weight or body mass index (BMI) are directly associated with blood pressure (BP), birth weight is inversely associated with BP. The scenario for height is more complex, as both tall and short stature have been associated with higher BP. We used a theoretical model treating sitting height (SH) and tissue masses (fat mass, lean mass) as components of metabolic load, and leg length (LL) as a marker of homeostatic metabolic capacity. We predicted that decreased capacity and increased load would be independently associated with increased BP.. METHODS AND RESULTS: Anthropometry, body composition (bio-electrical impedance analysis) and BP were measured in 601 adults (228 male) aged 20-91 years from three hill villages in southern Italy. Multiple regression analysis was used to investigate associations of body composition and anthropometry with BP. Adjusting for age, systolic BP (SBP) was associated with lean mass in males, and with adiposity in females, whereas diastolic BP (DBP) was associated with fat mass in both sexes. Associations of LL and SH with BP were in opposite directions. LL was inversely associated with SBP and DBP in males, with a similar trend evident in females. SH was directly associated with SBP and DBP in females, and with DBP in males. CONCLUSIONS: Consistent with our theoretical model, metabolic load is associated with increased BP, though differently between the sexes, whereas metabolic capacity is independently associated with lower BP. Our findings suggest that early growth improves hemodynamic tolerance of high metabolic load in adulthood..

Meta-analysis of genome-wide association studies for personality.

Personality can be thought of as a set of characteristics that influence people's thoughts, feelings and behavior across a variety of settings. Variation in personality is predictive of many outcomes in life, including mental health. Here we report on a meta-analysis of genome-wide association (GWA) data for personality in 10 discovery samples (17,375 adults) and five in silico replication samples (3294 adults). All participants were of European ancestry. Personality scores for Neuroticism, Extraversion, Openness to Experience, Agreeableness and Conscientiousness were based on the NEO Five-Factor Inventory. Genotype data of ≈ 2.4M single-nucleotide polymorphisms (SNPs; directly typed and imputed using HapMap data) were available. In the discovery samples, classical association analyses were performed under an additive model followed by meta-analysis using the weighted inverse variance method. Results showed genome-wide significance for Openness to Experience near the RASA1 gene on 5q14.3 (rs1477268 and rs2032794, P=2.8 × 10(-8) and 3.1 × 10(-8)) and for Conscientiousness in the brain-expressed KATNAL2 gene on 18q21.1 (rs2576037, P=4.9 × 10(-8)). We further conducted a gene-based test that confirmed the association of KATNAL2 to Conscientiousness. In silico replication did not, however, show significant associations of the top SNPs with Openness and Conscientiousness, although the direction of effect of the KATNAL2 SNP on Conscientiousness was consistent in all replication samples. Larger scale GWA studies and alternative approaches are required for confirmation of KATNAL2 as a novel gene affecting Conscientiousness.

Meta-analysis of genome-wide association studies identifies common variants in CTNNA2 associated with excitement-seeking.

The tendency to seek stimulating activities and intense sensations define excitement-seeking, a personality trait akin to some aspects of sensation-seeking. This trait is a central feature of extraversion and is a component of the multifaceted impulsivity construct. Those who score high on measures of excitement-seeking are more likely to smoke, use other drugs, gamble, drive recklessly, have unsafe/unprotected sex and engage in other risky behaviors of clinical and social relevance. To identify common genetic variants associated with the Excitement-Seeking scale of the Revised NEO Personality Inventory, we performed genome-wide association studies in six samples of European ancestry (N=7860), and combined the results in a meta-analysis. We identified a genome-wide significant association between the Excitement-Seeking scale and rs7600563 (P=2 × 10(-8)). This single-nucleotide polymorphism maps within the catenin cadherin-associated protein, alpha 2 (CTNNA2) gene, which encodes for a brain-expressed α-catenin critical for synaptic contact. The effect of rs7600563 was in the same direction in all six samples, but did not replicate in additional samples (N=5105). The results provide insight into the genetics of excitement-seeking and risk-taking, and are relevant to hyperactivity, substance use, antisocial and bipolar disorders.

Angiogenesis and biomarkers of cardiovascular risk in adults with metabolic syndrome.

OBJECTIVES: Metabolic syndrome (MetSyn) is associated with an increased risk of atherosclerosis and fatal cardiovascular events. Angiogenesis is thought to contribute to this risk as it might be involved in the progression of atherosclerotic plaques. We investigated the levels of circulating biomarkers of angiogenesis and cardiovascular risk in adults with MetSyn and assessed their association with established metabolic risk factors. DESIGN: The Genetic Park project is a highly inclusive cross-sectional survey (about 80% of residents) conducted in three isolated populations in Southern Italy. A total of 1000 men and women (age range: 18-98 years) were included in the analysis. Anthropometric and blood pressure measurements were recorded. Metabolic and cardiovascular biomarkers included glucose, triglycerides, total cholesterol, HDL, vascular endothelial growth factor, placental growth factor (PlGF), soluble fms-like tyrosine kinase-1, high-sensitivity C-reactive protein, high-sensitivity troponin T (hs-TnT) and N-terminal pro-brain natriuretic peptide (NT-proBNP). RESULT: Subjects with MetSyn had higher levels of PlGF and NT-proBNP after adjustment for age, smoking and body mass index. Circulating levels of PlGF, hs-TnT and NT-proBNP were directly related to the number of criteria of MetSyn, and this association interacted with gender. There was a strong correlation between ageing and cardiovascular risk. CONCLUSIONS: The increase in circulating levels of biomarkers of angiogenesis and cardiac function in subjects with MetSyn mirrors the pathophysiological changes occurring in the cardiovascular system. Over time, these changes might accelerate the formation and progression of atherosclerotic plaques and contribute significantly to cardiovascular morbidity and mortality risk.

Unbalanced expression of HLA-A and -B antigens: a specific feature of cutaneous melanoma and other non-hemopoietic malignancies reverted by IFN-gamma.

Conflicting evidences suggested that levels of HLA-A and -B antigens expressed on normal and neoplastic cells of given individuals are genetically predetermined, or, on the other hand, regulated by molecular mechanisms generating the down-regulated expression of HLA-B antigens frequently observed on melanoma cells. In our study, we quantitated, both at the protein and mRNA level, the amounts of HLA-A and -B antigens constitutively expressed on 23 primary cultures of metastatic melanomas and on autologous peripheral blood mononuclear cells (PBMC). Flow cytometric analyses identified a significantly (p < 0.01) lower expression of HLA-B antigens on melanoma cell cultures but not on autologous PBMC. Consistently, lower amounts of HLA-B antigens mRNA were detected by RNase protection assay exclusively in neoplastic cells. This unbalanced expression of HLA-A and -B antigens was readily reverted by interferon (IFN)-gamma but not by the DNA hypomethylating agent 5-aza-2'-deoxycytidine in 4 melanoma cell cultures investigated. Significantly (p < 0.05) lower levels of HLA-B antigens were also detected on cells from solid malignancies of different histotypes but not on neoplastic cells from hemopoietic neoplasms; levels of HLA-B antigens were rapidly up-regulated by IFN-gamma exclusively on non-hemopoietic transformed cells. Together, these data strongly argue against a genetic predetermination of the amounts of HLA-A and -B antigens expressed on normal and neoplastic cells of distinct melanoma patients and suggest that constitutively low levels of HLA-B antigens are a specific feature of non-hemopoietic transformed cells that is controlled by common regulatory mechanism(s) and that is possibly shared by non-hemopoietic normal cells.

Downstream sequence adjacent to AUG affects translation of chloramphenicol acetyl transferase in eukaryotic cells.

The CAT gene is widely used as a reporter in eukaryotic systems because of the efficient translation of its mRNA. We report here that a sequence occurring in the CAT mRNA at +15 nucleotides from CAT AUG is essential for translation. This sequence includes a stem-loop structure, which, however, exhibits a calculated stability significantly lower than that required for a hairpin to act as an enhancer of translation in vitro. Replacement of this region with the corresponding sequence from mRNAs that are normally translated in eukaryotic systems drastically reduced translation of CAT in COS cells, although the consensus sequence around the AUG, known to be required for high-level translation initiation, was conserved. These observations may be relevant for the exploitation of the CAT reporter system for analysis of the mechanisms of translation initiation by means of fusion constructs.

Biosynthesis and immunobiochemical characterization of gp17/GCDFP-15. A glycoprotein from seminal vesicles and from breast tumors, in HeLa cells and in Pichia pastoris yeast.

The gp17 factor is a secretory product of human seminal vesicle cells which binds to CD4 and acts as a potent inhibitor of T lymphocyte apoptosis induced by CD4 crosslinking and subsequent T-cell receptor (TCR) activation. The protein is identical to gross cystic disease fluid protein-15 (GCDFP-15), a breast tumor secretory marker PIP (prolactin inducible protein), a prolactin-controlled and androgen-controlled protein; secretory actin binding protein (SABP), a seminal plasma actin binding protein and extra-parotid glycoprotein (EP-GP), a secretory protein from the salivary gland. The structure of this protein has not yet been elucidated and no biological function has been clearly attributed to date. Expression of recombinant gp17/GCDFP-15 cDNA in bacteria and insect cells leads to the production of a misfolded insoluble protein. In this study, we describe the production of gp17/GCDFP-15 in two different eukaryotic systems, namely HeLa cells and the Pichia pastoris yeast. Using constructs in which gp17/GCDFP-15 was tagged with enhanced green fluorescent protein (EGFP) in various combinations, we observed expression only when the fusion protein was directed to the secretory compartment by the correct signal peptide. The resulting fluorescent protein was inefficiently secreted, thus suggesting that gp17/GCDFP-15 is not appropriately post-translationally processed and/or transported in HeLa cells. The use of the P. pastoris secretory pathway allowed instead the accumulation in the culture medium of a GCDFP-15/gp17 species which retained the ability to bind to CD4 and also most of the biochemical and immunological properties of the native protein. The production of an active recombinant molecule opens the way to correlate the structural properties of this peculiar factor to its ability to bind several proteins, including CD4, and to block CD4-mediated T cell programmed death.

Differential expression of insulin-dependent diabetes mellitus-associated HLA-DQA1 alleles in vivo.

The strong association of HLA-DQ genes with insulin-dependent diabetes mellitus (IDDM) susceptibility is persuasive evidence of their central role in the etiology of this autoimmune disease. Among other possibilities, it has been proposed that an unbalanced expression of IDDM-associated DQA, and/or DQB alleles may lead to alterations in the composition of alpha beta heterodimers and preferential expression of a particular heterodimer on the antigen-presenting cell surface, leading to self-recognition. In this report, we demonstrate the differential expression of DQA1 alleles in vivo, in particular of the two diabetogenic alleles DQA1*0301 and DQA1*0501. Family studies suggest that unequal HLA-DQA1 allele expression in heterozygous individuals is not associated in cis with the HLA-DQA1 gene, but may be affected by trans-acting determinant(s). We also discuss the segregation of this phenotype in IDDM-affected members. Furthermore, we examined historical samples of PBL from an IDDM-affected individual and an HLA-identical unaffected sibling acting in a kidney transplant program as donor and recipient, respectively. This analysis allowed us to establish that unbalanced expression of DQA1*0301 and DQA1*0501 can be induced by microenvironmental conditions. Inducible differential expression of HLA-DQA1 alleles may account for the discordance in the outcome of autoimmune disease in monozygotic twins and HLA-identical siblings.

Polymorphism in the 5' terminal region of the mRNA of HLA-DQA1 gene: identification of four groups of transcripts and their association with polymorphism in the alpha 1 domain.

Relative to other loci in the MHC, the HLA-DQ locus exhibits an exceptional degree of polymorphism of both A1 and B1 genes, particularly in the region coding for alpha and beta chains. Diversification of the association between different alpha and beta molecules either in cis or in trans contributes to the structural diversity of the repertoire of cell-surface class II protein's in the population. In addition, structural allelic polymorphisms in the 5' regulatory region of both DQB1 and DQA1 shows several linkage groups with respect to the allelic coding sequence of the respective genes. We describe here the allelic polymorphism in the DQA1 mRNA structure located at the 5' untranslated terminal region. This portion of the mRNA molecule represents, in many genes, a cis-acting regulatory sequence playing a role in the posttranscriptional mechanisms by which gene expression can be modulated. Based on detailed transcriptional analysis, we have been able to define at least four groups of transcripts in DQA1. The mRNA variability was associated with the polymorphism of the second exon of the DQA1 gene, coding for the alpha 1 domain and not with the DNA polymorphism in the 5' regulatory region.

Control of nucleo-cytoplasmic HLA-DRA mRNA partitioning by interaction of a retention signal with compartmentalized proteins.

Transport of mRNA from the nucleus to the cytoplasm is still a poorly understood process in which RNA signal sequences and cognate RNA-binding proteins may be involved. We have analysed the transport of the mRNA encoded by HLA-DRA, a member of the immunologically important MHC class II multigene family. We report that, in transient transfection experiments, HLA-DRA mRNA molecules encompassing a signal situated in the 3' untranslated region predominantly accumulate in the nucleus. We also show that the RNA sequence involved interacts with compartmentalized proteins of either nuclear or cytoplasmic origin. Deletion of the mRNA region encompassing this retention site results both in the abrogation of protein binding and in the release of HLA-DRA mRNA into the cytoplasm. In addition, we have found that the distribution of these HLA-DRA mRNA binding proteins is different in different cell types; in particular, their pattern of expression in Ntera-2, a human teratocarcinoma cell line, is distinct from that observed in Raji, a human B-lymphoma cell line, and is modulated by growth in retinoic acid. We conclude that recognition of a mRNA retention signal by proteins located in different compartments on either side of the nuclear membrane may regulate the nucleo-cytoplasmic partitioning of HLA-DRA transcripts and, perhaps, of MHC class II mRNA in general.