Search

Abstract Background Type 1 diabetes mellitus (T1DM) is a prototypic endocrine autoimmune disease resulting from an immune-mediated destruction of pancreatic insulin-secreting $$\beta$$ β cells. A comprehensive immune cell phenotype evaluation in T1DM has not been performed thus far at the single-cell level. Methods In this cross-sectional analysis, we generated a single-cell transcriptomic dataset of peripheral blood mononuclear cells (PBMCs) from 46 manifest T1DM (stage 3) cases and 31 matched controls. Results We surprisingly detected profound alterations in circulatory immune cells (1784 dysregulated genes in 13 immune cell types), far exceeding the count in the comparator systemic autoimmune disease SLE. Genes upregulated in T1DM were involved in WNT signaling, interferon signaling and migration of T/NK cells, antigen presentation by B cells, and monocyte activation. A significant fraction of these differentially expressed genes were also altered in T1DM pancreatic islets. We used the single-cell data to construct a T1DM metagene z-score (TMZ score) that distinguished cases and controls and classified patients into molecular subtypes. This score correlated with known prognostic immune markers of T1DM, as well as with drug response in clinical trials. Conclusions Our study reveals a surprisingly strong systemic dimension at the level of immune cell network in T1DM, defines disease-relevant molecular subtypes, and has the potential to guide non-invasive test development and patient stratification.

Abstract X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.

Abstract RNAseq data can be used to infer genetic variants, yet its use for estimating genetic population structure remains underexplored. Here, we construct a freely available computational tool (RGStraP) to estimate RNAseq-based genetic principal components (RG-PCs) and assess whether RG-PCs can be used to control for population structure in gene expression analyses. Using whole blood samples from understudied Nepalese populations and the Geuvadis study, we show that RG-PCs had comparable results to paired array-based genotypes, with high genotype concordance and high correlations of genetic principal components, capturing subpopulations within the dataset. In differential gene expression analysis, we found that inclusion of RG-PCs as covariates reduced test statistic inflation. Our paper demonstrates that genetic population structure can be directly inferred and controlled for using RNAseq data, thus facilitating improved retrospective and future analyses of transcriptomic data.

ABSTRACT A previous investigation has elucidated the landscape of Mtb genomic diversity and transmission dynamics in Ho Chi Minh City, Vietnam. Here, we expand the scope of this survey by adding a substantial number of additional genomes (total sample size: 2,542) and phenotypic drug susceptibility data for the majority of isolates. We aim to explore the prevalence and evolutionary dynamics of drug resistance and our ability to predict drug resistance from sequencing data. Among isolates tested phenotypically against first-line drugs, we observed high rates of streptomycin [STR, 37.7% (N = 573/1,520)] and isoniazid resistance [INH, 25.7% (N = 459/1,786)] and lower rates of resistance to rifampicin [RIF, 4.9% (N = 87/1,786)] and ethambutol [EMB, 4.2% (N = 75/1,785)]. Relative to global benchmarks, resistance to STR and INH was predicted accurately when applying the TB-Profiler algorithm to whole genome sequencing data (sensitivities of 0.81 and 0.87, respectively), while resistance to RIF and EMB was predicted relatively poorly (sensitivities of 0.70 and 0.44, respectively). Exploring the evolution of drug resistance revealed the main phylogenetic lineages to display differing dynamics and tendencies to evolve resistance via mutations in certain genes. The Beijing sublineage L2.2.1 was found to acquire de novo resistance mutations more frequently than isolates from other lineages and to suffer no apparent fitness cost acting to impede the transmission of resistance. Mutations conferring resistance to INH and STR arose earlier, on average, than those conferring resistance to RIF and are now more widespread across the phylogeny. The high prevalence of “background” INH resistance, combined with high rates of RIF mono-resistance (20.7%, N = 18/87), suggests that rapid assays for INH resistance will be valuable in this setting. These tests will allow the detection of INH mono-resistance and will allow multi-drug-resistant isolates to be distinguished from isolates with RIF mono-resistance. IMPORTANCE Drug-resistant tuberculosis (TB) infection is a growing and potent concern, and combating it will be necessary to achieve the WHO’s goal of a 95% reduction in TB deaths by 2035. While prior studies have explored the evolution and spread of drug resistance, we still lack a clear understanding of the fitness costs (if any) imposed by resistance-conferring mutations and the role that Mtb genetic lineage plays in determining the likelihood of resistance evolution. This study offers insight into these questions by assessing the dynamics of resistance evolution in a high-burden Southeast Asian setting with a diverse lineage composition. It demonstrates that there are clear lineage-specific differences in the dynamics of resistance acquisition and transmission and shows that different lineages evolve resistance via characteristic mutational pathways.

Introduction: Educational attainment, widely used in epidemiologic studies as a surrogate for socioeconomic status, is a predictor of cardiovascular health outcomes.Methods: A two-stage genome-wide meta-analysis of low-density lipoprotein cholesterol (LDL), high-density lipoprotein cholesterol (HDL), and triglyceride (TG) levels was performed while accounting for gene-educational attainment interactions in up to 226,315 individuals from five population groups. We considered two educational attainment variables: “Some College” (yes/no, for any education beyond high school) and “Graduated College” (yes/no, for completing a 4-year college degree). Genome-wide significant (p < 5 × 10−8) and suggestive (p < 1 × 10−6) variants were identified in Stage 1 (in up to 108,784 individuals) through genome-wide analysis, and those variants were followed up in Stage 2 studies (in up to 117,531 individuals).Results: In combined analysis of Stages 1 and 2, we identified 18 novel lipid loci (nine for LDL, seven for HDL, and two for TG) by two degree-of-freedom (2 DF) joint tests of main and interaction effects. Four loci showed significant interaction with educational attainment. Two loci were significant only in cross-population analyses. Several loci include genes with known or suggested roles in adipose (FOXP1, MBOAT4, SKP2, STIM1, STX4), brain (BRI3, FILIP1, FOXP1, LINC00290, LMTK2, MBOAT4, MYO6, SENP6, SRGAP3, STIM1, TMEM167A, TMEM30A), and liver (BRI3, FOXP1) biology, highlighting the potential importance of brain-adipose-liver communication in the regulation of lipid metabolism. An investigation of the potential druggability of genes in identified loci resulted in five gene targets shown to interact with drugs approved by the Food and Drug Administration, including genes with roles in adipose and brain tissue.Discussion: Genome-wide interaction analysis of educational attainment identified novel lipid loci not previously detected by analyses limited to main genetic effects.

Abstract Background The popular statistics-based Genome-wide association studies (GWAS) have provided deep insights into the field of complex disorder genetics. However, its clinical applicability to predict disease/trait outcomes remains unclear as statistical models are not designed to make predictions. This study employs statistics-free machine-learning (ML)-optimized polygenic risk score (PRS) to complement existing GWAS and bring the prediction of disease/trait outcomes closer to clinical application. Rheumatoid Arthritis (RA) was selected as a model disease to demonstrate the robustness of ML in disease prediction as RA is a prevalent chronic inflammatory joint disease with high mortality rates, affecting adults at the economic prime. Early identification of at-risk individuals may facilitate measures to mitigate the effects of the disease. Methods This study employs a robust ML feature selection algorithm to identify single nucleotide polymorphisms (SNPs) that can predict RA from a set of training data comprising RA patients and population control samples. Thereafter, selected SNPs were evaluated for their predictive performances across 3 independent, unseen test datasets. The selected SNPs were subsequently used to generate PRS which was also evaluated for its predictive capacity as a sole feature. Results Through robust ML feature selection, 9 SNPs were found to be the minimum number of features for excellent predictive performance (AUC > 0.9) in 3 independent, unseen test datasets. PRS based on these 9 SNPs was significantly associated with (P 0.9) of RA in the 3 unseen datasets. A RA ML-PRS calculator of these 9 SNPs was developed ( https://xistance.shinyapps.io/prs-ra/ ) to facilitate individualized clinical applicability. The majority of the predictive SNPs are protective, reside in non-coding regions, and are either predicted to be potentially functional SNPs (pfSNPs) or in high linkage disequilibrium (r2 > 0.8) with un-interrogated pfSNPs. Conclusions These findings highlight the promise of this ML strategy to identify useful genetic features that can robustly predict disease and amenable to translation for clinical application.

Abstract Background Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery. Results To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3–5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism. Conclusions Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

Abstract Background Skeletal dysplasia is a heterogeneous group of disorders. Spondyloepiphyseal dysplasias comprise one subgroup. Deficiency of carbohydrate sulfotransferase 3 has been reported in a small number of patients with recessively inherited spondyloepiphyseal dysplasia with joint dislocation, short stature and scoliosis. We report here molecular and clinical findings of affected individuals in three consanguineous Pakistani families. Affected individuals in all three families had a uniform phenotype including severe short stature, multiple dislocated joints, progressive scoliosis and facial dysmorphism. Methods Clinical evaluation was done for three unrelated families. Radiological survey of bones was completed for patients from two of the families. Whole exome sequencing index patients from each family was performed followed by Sanger sequencing for validation of segregation of identified variants in respective families. In-silico analysis for determining pathogenicity of identified variants and conservation was done. Results Whole-exome sequencing revealed biallelic variants c.590 T > C;p.(Leu197Pro), c.603C > A;p.(Tyr201Ter) and c.661C > T;p.(Arg221Cys) in CHST3 (NM_004273.5) in the three families with eight, five and two affected individuals, respectively. Contrary to previous reports, affected individuals in none of the families exhibited a hearing loss. Conclusion We describe genotypic and phenotypic findings of three unrelated families with spondyloepiphyseal dysplasia. Our study confirms phenotypic variability and adds to the genotypic spectrum of spondyloepiphyseal dysplasia.

Summary: Genomic researchers increasingly utilize commercial cloud service providers (CSPs) to manage data and analytics needs. CSPs allow researchers to grow Information Technology (IT) infrastructure on demand to overcome bottlenecks when combining large datasets. However, without adequate security controls, the risk of unauthorized access may be higher for data stored on the cloud. Additionally, regulators are mandating data access patterns and specific security protocols for the storage and use of genomic data. While CSP provides tools for security and regulatory compliance, building the necessary controls required for cloud solutions is not trivial. Research Assets Provisioning and Tracking Online Repository (RAPTOR) by the Genome Institute of Singapore is a cloud-native genomics data repository and analytics platform that implements a “five-safes” framework to provide security and governance controls to data contributors and users, leveraging CSP for sharing and analysis of genomic datasets without the risk of security breaches or running afoul of regulations.

Background: To evaluate the utility of polygenic risk scores (PRSs) in identifying high-risk individuals, different publicly available PRSs for breast (n=85), prostate (n=37), colorectal (n=22), and lung cancers (n=11) were examined in a prospective study of 21,694 Chinese adults. Methods: We constructed PRS using weights curated in the online PGS Catalog. PRS performance was evaluated by distribution, discrimination, predictive ability, and calibration. Hazard ratios (HR) and corresponding confidence intervals (CI) of the common cancers after 20 years of follow-up were estimated using Cox proportional hazard models for different levels of PRS. Results: A total of 495 breast, 308 prostate, 332 female-colorectal, 409 male-colorectal, 181 female-lung, and 381 male-lung incident cancers were identified. The area under receiver operating characteristic curve for the best-performing site-specific PRS were 0.61 (PGS000873, breast), 0.70 (PGS00662, prostate), 0.65 (PGS000055, female-colorectal), 0.60 (PGS000734, male-colorectal), 0.56 (PGS000721, female-lung), and 0.58 (PGS000070, male-lung), respectively. Compared to the middle quintile, individuals in the highest cancer-specific PRS quintile were 64% more likely to develop cancers of the breast, prostate, and colorectal. For lung cancer, the lowest cancer-specific PRS quintile was associated with 28–34% decreased risk compared to the middle quintile. In contrast, the HR observed for quintiles 4 (female-lung: 0.95 [0.61–1.47]; male-lung: 1.14 [0.82–1.57]) and 5 (female-lung: 0.95 [0.61–1.47]) were not significantly different from that for the middle quintile. Conclusions: Site-specific PRSs can stratify the risk of developing breast, prostate, and colorectal cancers in this East Asian population. Appropriate correction factors may be required to improve calibration. Funding: This work is supported by the National Research Foundation Singapore (NRF-NRFF2017-02), PRECISION Health Research, Singapore (PRECISE) and the Agency for Science, Technology and Research (A*STAR). WP Koh was supported by National Medical Research Council, Singapore (NMRC/CSA/0055/2013). CC Khor was supported by National Research Foundation Singapore (NRF-NRFI2018-01). Rajkumar Dorajoo received a grant from the Agency for Science, Technology and Research Career Development Award (A*STAR CDA - 202D8090), and from Ministry of Health Healthy Longevity Catalyst Award (HLCA20Jan-0022).

Abstract Background Family history, and genetic and non-genetic risk factors can stratify women according to their individual risk of developing breast cancer. The extent of overlap between these risk predictors is not clear. Methods In this case-only analysis involving 7600 Asian breast cancer patients diagnosed between age 30 and 75 years, we examined identification of high-risk patients based on positive family history, the Gail model 5-year absolute risk [5yAR] above 1.3%, breast cancer predisposition genes (protein-truncating variants [PTV] in ATM, BRCA1, BRCA2, CHEK2, PALB2, BARD1, RAD51C, RAD51D, or TP53), and polygenic risk score (PRS) 5yAR above 1.3%. Results Correlation between 5yAR (at age of diagnosis) predicted by PRS and the Gail model was low (r=0.27). Fifty-three percent of breast cancer patients (n=4041) were considered high risk by one or more classification criteria. Positive family history, PTV carriership, PRS, or the Gail model identified 1247 (16%), 385 (5%), 2774 (36%), and 1592 (21%) patients who were considered at high risk, respectively. In a subset of 3227 women aged below 50 years, the four models studied identified 470 (15%), 213 (7%), 769 (24%), and 325 (10%) unique patients who were considered at high risk, respectively. For younger women, PRS and PTVs together identified 745 (59% of 1276) high-risk individuals who were not identified by the Gail model or family history. Conclusions Family history and genetic and non-genetic risk stratification tools have the potential to complement one another to identify women at high risk.

Abstract Background and Aim Chronic hepatitis C virus (HCV) infection, long‐term alcohol use, cigarette smoking, and obesity are the major risk factors for hepatocellular carcinoma (HCC) in the United States, but the disease risk varies substantially among individuals with these factors, suggesting host susceptibility to and gene–environment interactions in HCC. To address genetic susceptibility to HCC, we conducted a genome‐wide association study (GWAS). Methods Two case‐control studies on HCC were conducted in the United States. DNA samples were genotyped using the Illumian microarray chip with over 710 000 single nucleotide polymorphisms (SNPs). We compared these SNPs between 705 HCC cases and 1455 population controls for their associations with HCC and verified our findings in additional studies. Results In this GWAS, we found that two SNPs were associated with HCC at P

RhoGTPase regulators play a key role in the development of the nervous system, and their dysfunction can result in brain malformation and associated disorders. Several guanine nucleotide exchange factors (GEF) have been linked to neurodevelopmental disorders. In line with this, ARHGEF17 has been recently linked as a risk gene to intracranial aneurysms. Here we report siblings of a consanguineous Pakistani family with biallelic variants in the ARHGEF17 gene associated with a neurodevelopmental disorder with intellectual disability, speech delay and motor dysfunction but not aneurysms. Cranial MRI performed in one patient revealed generalized brain atrophy with an enlarged ventricular system, thin corpus callosum and microcephaly. Whole exome sequencing followed by Sanger sequencing in two of the affected individuals revealed a homozygous missense variant (g.11:73021307, c.1624C>T (NM_014786.4), p.R542W) in the ARHGEF17 gene. This variant is in a highly conserved DCLK1 phosphorylation consensus site (I/L/V/F/M]RRXX[pS/pT][I/L/M/V/F) of the protein. Our report expands the phenotypic spectrum of ARHGEF17 variants from increased intracranial aneurysm risk to neurodevelopmental disease and thereby add ARHGEF17 to the list of GEF genes involved in neurodevelopmental disorders.

Xuling Chang et al. study whether low frequency variants are associated with leukocyte telomere length and whether these variants predispose to incident cancers and mortalities among East Asians. They find that three East Asian specific variants at POT1, TERF1 and STN1 loci are associated with leukocyte telomere length and report a mechanistic pathway linking TERF1, leukocyte telomere length and incident colon cancer.

Primary open-angle glaucoma (POAG) is highly heritable, yet not well understood from a genetic perspective. Here, the authors perform a meta-analysis of genome-wide association studies in 34,179 POAG cases, identifying 44 previously unreported risk loci and mapping effects across multiple ethnicities.

Here, the authors report a multi-ethnic genome wide association meta-analysis of 12 studies from the International Cataract Genetics Consortium. They find six new loci associated with age-related nuclear cataract, in addition to replicating the association at CRYAA, and suggest a strong genetic link between age-related nuclear and congenital cataracts.

Abstract Background Shorter telomere length (TL) has been associated with poor health behaviors, increased risks of chronic diseases and early mortality. Excessive shortening of telomere is a marker of accelerated aging and can be influenced by oxidative stress and nutritional deficiency. Plasma n6:n3 polyunsaturated fatty acid (PUFA) ratio may impact cell aging. Increased dietary intake of marine n-3 PUFA is associated with reduced telomere attrition. However, the effect of plasma PUFA on leukocyte telomere length (LTL) and its interaction with genetic variants are not well established. Methods A nested coronary artery disease (CAD) case-control study comprising 711 cases and 638 controls was conducted within the Singapore Chinese Health Study (SCHS). Samples genotyped with the Illumina ZhongHua-8 array. Plasma n-3 and n-6 PUFA were quantified using mass spectrometry (MS). LTL was measured with quantitative PCR method. Linear regression was used to test the association between PUFA and LTL. The interaction between plasma PUFAs and genetic variants was assessed by introducing an additional term (PUFA×genetic variant) in the regression model. Analysis was carried out in cases and controls separately and subsequently meta-analyzed using the inverse-variance weighted method. We further assessed the association of PUFA and LTL with CAD risk by Cox Proportional-Hazards model and whether the effect of PUFA on CAD was mediated through LTL by using structural equation modeling. Results Higher n6:n3 ratio was significantly associated with shorter LTL (p = 0.018) and increased CAD risk (p = 0.005). These associations were mainly driven by elevated plasma total n-3 PUFAs, especially eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) (p

Some immune-mediated diseases may originate in early childhood. The authors mapped eQTLs and response eQTLs to various stimuli in neonatal myeloid cells and T cells, and revealed their potential role in immune-mediated diseases using colocalisation and Mendelian randomisation.

Polygenic risk scores predict the likelihood that an individual will develop a certain cancer, however these are often specific for a given population. Here, the authors show that a risk score developed to assess the risk of breast cancer in European women can also predict risk in Asian populations.

Inbreeding depression has been observed in many different species, but in humans a systematic analysis has been difficult so far. Here, analysing more than 1.3 million individuals, the authors show that a genomic inbreeding coefficient (FROH) is associated with disadvantageous outcomes in 32 out of 100 traits tested.

Shortening of leukocyte telomere length (LTL) is associated with age and increased risk for various chronic diseases. Here, the authors report genome-wide association studies for LTL in Singaporean Chinese populations and find that longer LTL associates with less severe outcomes of respiratory disease phenotypes.

Purpose: To describe a case of plateau iris syndrome (PIS) and angle-closure glaucoma (ACG) in a patient with nail-patella syndrome (NPS). Observation: A 33 year-old woman of Slovakian ancestry from Norway with a history of NPS presented with angle-closure secondary to plateau iris. At the time of her NPS diagnosis, she had no ocular pathology. Genetic testing revealed a rare de novo mutation in LMX1B [c.668G>C (p.Arg223Pro)]. Two years later, she experienced acute bilateral ocular pain and blurred vision in the setting of one year of reported visual loss. Subsequent ophthalmic examinations revealed closed angles and plateau iris OU with ACG OD and angle-closure OS. Perimetry showed superonasal visual field defects OD and no defects OS. Ocular coherence tomography (OCT) revealed thinning of the inferior pole of the optic nerve OD. Medical management proved ineffective. A laser peripheral iridotomy (LPI) OD was performed, without resolution of the angle-closure, and a diagnosis of plateau iris syndrome (PIS) was made. She was then treated with an argon laser peripheral iridoplasty (ALPI) and clear lens extraction with a posterior chamber intraocular lens (PCIOL) and goniosynechialysis OD, but her IOP remained elevated OU. She was referred to New York Eye and Ear Infirmary of Mount Sinai, where an LPI OS was performed, but angle-closure persisted, consistent with PIS. An ALPI OS with touch-up was performed, and her IOP normalized with dark-room gonioscopy revealing open angles OU. Conclusions And importance: NPS has been associated with ocular hypertension (OHTN) and open-angle glaucoma (OAG); however, to our knowledge, no association between NPS and angle-closure has previously been reported. The case described here, of a patient with a rare de novo mutation and ocular findings of PIS with associated ACG, represents a novel genetic and clinical presentation of NPS.

Abstract Introduction Osteogenesis imperfecta (OI) is a clinically and genetically heterogeneous disease with skeletal fragility and variable extra-skeletal manifestations. To date several point mutations in 18 different genes causing different types of OI have been identified. Mutations in WNT1 compromise activity of the osteoblasts leading to disturbed bone mass accrual, fragility fractures and progressive skeletal abnormalities. The present study was conducted to determine the underlying genetic cause of an autosomal recessive skeletal dysplasia in a large consanguineous family from Chinute, Pakistan. Materials and methods Blood was collected from 24 individuals of affected family along with clinical data. Homozygosity mapping was performed to confirm consanguinity. SNPs were identified, followed by whole exome and Sanger sequencing. In silico characterization of WNT1 mutation was performed using multiple platforms. Results Nine affected family members exhibited severe bone deformities, recurrent fractures, short stature and low bone mineral density. SNP array data revealed homozygous segments > 1 Mb in length accounting for 2.1–12.7% of the genome in affected individuals and their siblings and a single 6,344,821 bp homozygous region in all affected individuals on chromosome 12q12-q13. This region includes two potential OI candidate genes WNT1 and VDR. We did whole-exome sequencing for both genes in two patients and identified a novel damaging missense mutation in exon 4 of WNT1: c.1168G > T (NM_005430) resulting in p.G324C. Sanger sequencing confirmed segregation of mutation with the disease in family. Conclusion We report a novel mutation responsible for OI and our investigation expands the spectrum of disease-causing WNT1 mutations and the resulting OI phenotypes.

Reduced central corneal thickness (CCT) is observed in common eye diseases as well as in rare Mendelian disorders. Here, in a cross-ancestry GWAS, the authors identify 19 novel genetic loci associated with CCT, a subset of which is involved in rare corneal or connective tissue disorders.

Emmanuelle Souzeau, who contributed to analysis of data, was inadvertently omitted from the author list in the originally published version of this Article. This has now been corrected in both the PDF and HTML versions of the Article.

Abstract Background Osteopetrosis is a rare inherited bone disorder mainly described as an increased bone density caused by defective osteoclastic bone resorption. To date, genetic variants of eleven genes have been reported so far to be associated with different types of osteopetrosis. However, malignant infantile osteopetrosis, a lethal form of the disease, is mostly (50%) caused by mutation(s) in TCIRG1 gene. In this study, we investigated a consanguineous Pakistani family clinically and genetically to elucidate underlying molecular basis of the infantile osteopetrosis. Methods DNA samples from five family members were subjected to SNP-array based whole genome homozygosity mapping. Data was analyzed and potentially pathogenic mutation was identified by Sanger sequencing of two affected as well as three phenotypically healthy individuals in the family. The significance of identified pathogenic variation and its impact on protein structure and function was studied using various bioinformatics tools. Results DNA samples from five family members were subjected to genome-wide SNP array genotyping and homozygosity mapping which identified ~4 Mb region on chr11 harboring the TCIRG1 gene. Sanger sequencing unveiled a novel homozygous deletion c. 624delC in exon 6 of the TCIRG1 gene encodes a3 subunit of V-ATPase complex. The identified deletion resulted in a frame shift producing a truncated protein of 208 aa. In silico analysis of premature termination of the a3 subunit of V-ATPase complex revealed deleterious effects on the protein structure, predicting impaired or complete loss of V-ATPase function causing infantile osteopetrosis. Conclusions Since a3 subunit of V-ATPase complex plays a crucial role in bone resorption process, structurally abnormal a3 subunit might have adversely affected bone resorption process, leading to infantile osteopetrosis in Pakistani family. Therefore, the present study not only expands the genotypic spectrum of osteopetrosis but also improve understandings of the role of V-ATPase a3 subunit in bone resorption process. Moreover, our findings should help in genetic counseling and provide further insight into the disease pathogenesis and potential targeted therapy.

Abstract Background Preliminary evidence suggests that changes in plasma brain-derived neurotrophic factor (BDNF) levels may contribute to the occurrence of chemotherapy-associated cognitive impairment (CACI), and a previous study suggested that carriers of the BDNF Met homozygous genotype are protected from CACI. Methods This multicenter, prospective cohort study involved chemotherapy-receiving early-stage breast cancer (ESBC) patients. Self-perceived cognitive function was longitudinally assessed using the validated FACT-Cog (ver. 3) across three time points: Prior to chemotherapy (T1), during chemotherapy (T2), and at the end of chemotherapy (T3). Plasma BDNF levels were quantified using enzyme-linked immunosorbent assay. Genotyping was performed using Sanger Sequencing. Results A total of 51 chemotherapy-receiving ESBC patients (mean age: 52.6 ± 9.5 years) were recruited, and 11 patients (21.6%) reported subjective cognitive impairment post-chemotherapy. Overall, there was a reduction in median plasma BDNF levels over time (T1: 5423.0 pg/ml; T2: 5313.6 pg/ml; T3: 4050.3 pg/ml; p

Disturbance in lipid metabolism has been suggested as a major pathogenic factor for age-related macular degeneration (AMD). Conventional lipid measures have been inconsistently associated with AMD. Other factors that can alter lipid metabolism include lipoprotein phenotype and genetic mutations. We performed a case-control study to examine the association between lipoprotein profile and neovascular AMD (nAMD) and whether the cholesterylester transfer protein (CETP) D442G mutation modulates these associations. Patients with nAMD had significantly higher concentrations of HDL and IDL compared with controls. The increase in HDL particles in nAMD patients was driven by an excess of medium-sized particles. Concurrently, patients with nAMD also had lower Apo A-1, lower VLDL and chylomicron lipoprotein. Many of these associations showed a dose-dependent association between controls, early AMD cases, and nAMD cases. Adjustment for the presence of the D442G mutation at the CETP locus did not significantly alter the increased AMD risk associated with HDL particle concentration. AMD is associated with variation in many lipoprotein subclasses, including increased HDL and IDL particles and decreased Apo A-1, VLDL, and chylomicron particles. These data suggest widespread systemic disturbance in lipid metabolism in the pathogenesis of AMD, including possible alterations in lipoprotein carrier capacity.

Abstract Background A previous genome-wide association study identified 2 susceptibility loci for severe dengue at MICB rs3132468 and PLCE1 rs3740360 and further work showed these mutations to be also associated with less severe clinical presentations. The aim of this study was to determine if these specific loci were associated with laboratory features of dengue that correlate with clinical severity with the aim of elucidating the functional basis of these genetic variants. Methods This was a case-only analysis of laboratory-confirmed dengue patients obtained from 2 prospective cohort studies and 1 randomised clinical trial in Vietnam (Trial registration: ISRCTN ISRCTN03147572. Registered 24th July 2012). 2742 dengue cases were successfully genotyped at MICB rs3132468 and PLCE1 rs3740360. Laboratory variables were compared between genotypes and stratified by DENV serotype. Results The analysis showed no association between MICB and PLCE1 genotype and early viraemia level, platelet nadir, white cell count nadir, or maximum haematocrit in both overall analysis and in analysis stratified by serotype. Discussion The lack of an association between genotype and viremia level may reflect the sampling procedures within the included studies. The study findings mean that the functional basis of these mutations remains unclear. Trial registration ISRCTN ISRCTN03147572 . Registered 24th July 2012.

LOXL1 is a genetic risk factor for pseudoexfoliation syndrome of the eye but a causal variant has not been identified. Here, Pasutto et al., find intronic LOXL1 risk variants influence transcription factor binding and alternative splicing of LOXL1 in affected tissues reducing levels of LOXL1mRNA.

Hepatocellular carcinoma has one of the poorest survival rates amongst cancers. Here, the authors highlight the intra-tumour heterogeneity of this disease, finding that spatially closer tumour sectors are genetically more similar and that intra-hepatic metastasis is accompanied by rapid genetic diversification.

A consanguineous family from Pakistan was ascertained to have a novel deafness-dystonia syndrome with motor regression, ichthyosis-like features and signs of sensory neuropathy. By applying a combined strategy of linkage analysis and whole-exome sequencing in the presented family, a homozygous nonsense mutation, c.4G>T (p.Glu2*), in FITM2 was identified. FITM2 and its paralog FITM1 constitute an evolutionary conserved protein family involved in partitioning of triglycerides into cellular lipid droplets. Despite the role of FITM2 in neutral lipid storage and metabolism, no indications for lipodystrophy were observed in the affected individuals. In order to obtain independent evidence for the involvement of FITM2 in the human pathology, downregulation of the single Fitm ortholog, CG10671, in Drosophila melanogaster was pursued using RNA interference. Characteristics of the syndrome, including progressive locomotor impairment, hearing loss and disturbed sensory functions, were recapitulated in Drosophila, which supports the causative nature of the FITM2 mutation. Mutation-based genetic counseling can now be provided to the family and insight is obtained into the potential impact of genetic variation in FITM2.

Summary: Vaccines have revolutionized public health, preventing millions of deaths each year, particularly in childhood. Yet, there is considerable variability in the magnitude and persistence of vaccine-induced immunity. Maintenance of specific antibody is essential for continuity of vaccine-induced serological protection. We conducted a genome-wide association study into the persistence of immunity to three childhood vaccines: capsular group C meningococcal (MenC), Haemophilus influenzae type b, and tetanus toxoid (TT) vaccines. We detail associations between variants in a locus containing a family of signal-regulatory proteins and the persistence MenC immunity. We postulate a regulatory role for the lead SNP, with supporting epigenetic and expression quantitative trait loci data. Furthermore, we define associations between SNPs in the human leukocyte antigen (HLA) locus and the persistence of TT-specific immunity. Moreover, we describe four classical HLA alleles, HLA DRB1∗0301, HLA DQB1∗0201, HLA DQB1∗0602, and HLA DRB1∗1501, associated with TT-specific immunity, independent of the lead SNP association. : Genetic variants impact immune responses to antigenic stimuli. O’Connor et al. show that variants within the human leukocyte antigen (HLA) and signal-regulatory proteins loci are associated with persistence of immunity following immunization, suggesting these genes are involved in regulating immune responses to routine childhood immunization. Keywords: vaccine-induced immunity, GWAS, persistence of immunity, meningococcal disease, HLA, SIRPG

We report on three new patients with spondyloocular syndrome (SOS) in a consanguineous Pakistani family. All three patients present progressive generalized osteoporosis, short stature, recurrent fractures, hearing loss and visual impairments. WES revealed a novel homozygous frameshift variant in exon 11 of XYLT2 (NG 012175.1, NP_071450.2) resulting in loss of evolutionary conserved amino acid sequences (840 – 865/865) at C-terminus p.R840fs∗115. Sanger Sequencing confirmed the presence of the novel homozygous mutation in all three patients while the parents were heterozygous carriers of the mutation, in accordance with an autosomal recessive inheritance pattern. Only nine variants worldwide have previously been reported in XYLT2 in patients with SOS phenotype. These three patients with novel homozygous variant extend the genotypic and phenotypic spectrum of SOS.

Barrett’s oesophagus is a precancerous intestinal metaplasia that can progress to oesophageal adenocarcinoma. In this study, the authors isolate and characterize human Barrett’s stem cells and identify a specific genomic pedigree that supports the potential role of these cells as precursors of oesophageal adenocarcinoma.

Ethnic diversity has been long considered as one of the factors explaining why the severe forms of dengue are more prevalent in Southeast Asia than anywhere else. Here we take advantage of the admixed profile of Southeast Asians to perform coupled association-admixture analyses in Thai cohorts. For dengue shock syndrome (DSS), the significant haplotypes are located in genes coding for phospholipase C members (PLCB4 added to previously reported PLCE1), related to inflammation of blood vessels. For dengue fever (DF), we found evidence of significant association with CHST10, AHRR, PPP2R5E and GRIP1 genes, which participate in the xenobiotic metabolism signaling pathway. We conducted functional analyses for PPP2R5E, revealing by immunofluorescence imaging that the coded protein co-localizes with both DENV1 and DENV2 NS5 proteins. Interestingly, only DENV2-NS5 migrated to the nucleus, and a deletion of the predicted top-linking motif in NS5 abolished the nuclear transfer. These observations support the existence of differences between serotypes in their cellular dynamics, which may contribute to differential infection outcome risk. The contribution of the identified genes to the genetic risk render Southeast and Northeast Asian populations more susceptible to both phenotypes, while African populations are best protected against DSS and intermediately protected against DF, and Europeans the best protected against DF but the most susceptible against DSS.

Heavy alcohol consumption is an established risk factor for hypertension; the mechanism by which alcohol consumption impact blood pressure (BP) regulation remains unknown. We hypothesized that a genome-wide association study accounting for gene-alcohol consumption interaction for BP might identify additional BP loci and contribute to the understanding of alcohol-related BP regulation. We conducted a large two-stage investigation incorporating joint testing of main genetic effects and single nucleotide variant (SNV)-alcohol consumption interactions. In Stage 1, genome-wide discovery meta-analyses in ≈131K individuals across several ancestry groups yielded 3,514 SNVs (245 loci) with suggestive evidence of association (P < 1.0 x 10-5). In Stage 2, these SNVs were tested for independent external replication in ≈440K individuals across multiple ancestries. We identified and replicated (at Bonferroni correction threshold) five novel BP loci (380 SNVs in 21 genes) and 49 previously reported BP loci (2,159 SNVs in 109 genes) in European ancestry, and in multi-ancestry meta-analyses (P < 5.0 x 10-8). For African ancestry samples, we detected 18 potentially novel BP loci (P < 5.0 x 10-8) in Stage 1 that warrant further replication. Additionally, correlated meta-analysis identified eight novel BP loci (11 genes). Several genes in these loci (e.g., PINX1, GATA4, BLK, FTO and GABBR2) have been previously reported to be associated with alcohol consumption. These findings provide insights into the role of alcohol consumption in the genetic architecture of hypertension.

A low level of HDL cholesterol (HDL-C) is a common clinical scenario and an important marker for increased cardiovascular risk. Many patients with very low or very high HDL-C have a rare mutation in one of several genes, but identification of the molecular abnormality in patients with extreme HDL-C is rarely performed in clinical practice. We investigated the accuracy and diagnostic yield of a targeted next-generation sequencing (NGS) assay for extreme levels of HDL-C. We developed a targeted NGS panel to capture the exons, intron/exon boundaries, and untranslated regions of 26 genes with highly penetrant effects on plasma lipid levels. We sequenced 141 patients with extreme HDL-C levels and prioritized variants in accordance with medical genetics guidelines. We identified 35 pathogenic and probably pathogenic variants in HDL genes, including 21 novel variants, and performed functional validation on a subset of these. Overall, a molecular diagnosis was established in 35.9% of patients with low HDL-C and 5.2% with high HDL-C, and all prioritized variants identified by NGS were confirmed by Sanger sequencing. Our results suggest that a molecular diagnosis can be identified in a substantial proportion of patients with low HDL-C using targeted NGS.

In the original publication of this article [1], there are two errors in the article which the cDNA position of the pathogenic variant WNT1 p.Gly324Cys should be c.970G>T instead of c.1168G>T.

Multi-drug and extensively drug-resistant tuberculosis (MDR and XDR-TB) are problems that threaten public health worldwide. Only some genetic markers associated with drug-resistant TB are known. Whole-genome sequencing (WGS) is a promising tool for distinguishing between re-infection and persistent infection in isolates taken at different times from a single patient, but has not yet been applied in MDR and XDR-TB. We aim to detect genetic markers associated with drug resistance and distinguish between reinfection and persistent infection from MDR and XDR-TB patients based on WGS analysis. Samples of Mycobacterium tuberculosis (n = 7), serially isolated from 2 MDR cases and 1 XDR-TB case, were retrieved from Siriraj Hospital, Bangkok. The WGS analysis used an Illumina Miseq sequencer. In cases of persistent infection, MDR-TB isolates differed at an average of 2 SNPs across the span of 2-9 months whereas in the case of reinfection, isolates differed at 61 SNPs across 2 years. Known genetic markers associated with resistance were detected from strains susceptible to streptomycin (2/7 isolates), p-aminosalicylic acid (3/7 isolates) and fluoroquinolone drugs. Among fluoroquinolone drugs, ofloxacin had the highest phenotype-genotype concordance (6/7 isolates), whereas gatifloxcain had the lowest (3/7 isolates). A putative candidate SNP in Rv2477c associated with kanamycin and amikacin resistance was suggested for further validation. WGS provided comprehensive results regarding molecular epidemiology, distinguishing between persistent infection and reinfection in M/XDR-TB and potentially can be used for detection of novel mutations associated with drug resistance.

To determine the differential effects of genetic polymorphism in CFH and ARMS2 on risk of age-related macular degeneration (AMD) between phakic vs. pseudophakic/aphakic eyes.9,529 eyes of 4,918 participants from the Singapore Malay Eye Study and Singapore Indian Eye Study were analyzed. Participants had detailed eye examinations, including slit-lamp examinations and dilated fundus photography. AMD grading was performed according to the Wisconsin age-related maculopathy grading system. Lens status was defined. Single nucleotide polymorphisms (SNPs) rs10801555 (Y402H) within CFH and rs3750847 in ARMS2 were assessed. The main outcome measure was early AMD or any AMD.No significant associations between the CFH Y402H genotypes and early AMD were found in phakic individuals. In contrast, among pseudophakic/aphakic individuals, the CFH Y402H risk genotypes were significantly associated with higher odds of early AMD, with an OR of 1.57 (95% CI: 1.07-2.29) for GA genotype and 2.40 (95% CI: 1.25-4.61) for AA genotype, compared to those with GG genotype. There was significant interaction between pseudophakic/aphakic status and CFH Y402H variant on risk of early AMD (p = 0.037), adjusting for age, gender, and the first 5 genetic principal components. No significant interaction was found between lens status and ARMS2 rs3750847.CFH genetic polymorphism and pseudophakic/aphakic status may have a potential synergistic effect on early AMD, suggesting roles for the complement system and related pathways in the pathogenesis of AMD in eyes after cataract surgery.

Anterior chamber depth (ACD) is a key anatomical risk factor for primary angle closure glaucoma (PACG). We conducted a genome-wide association study (GWAS) on ACD to discover novel genes for PACG on a total of 5,308 population-based individuals of Asian descent. Genome-wide significant association was observed at a sequence variant within ABCC5 (rs1401999; per-allele effect size = -0.045 mm, P = 8.17 × 10(-9)). This locus was associated with an increase in risk of PACG in a separate case-control study of 4,276 PACG cases and 18,801 controls (per-allele OR = 1.13 [95% CI: 1.06-1.22], P = 0.00046). The association was strengthened when a sub-group of controls with open angles were included in the analysis (per-allele OR = 1.30, P = 7.45 × 10(-9); 3,458 cases vs. 3,831 controls). Our findings suggest that the increase in PACG risk could in part be mediated by genetic sequence variants influencing anterior chamber dimensions.

The influenza A virus is an RNA virus that is responsible for seasonal epidemics worldwide with up to five million cases of severe illness and 500,000 deaths annually according to the World Health Organization estimates. The factors associated with severe diseases are not well defined, but more severe disease is more often seen among persons aged >65 years, infants, pregnant women, and individuals of any age with underlying health conditions.Using gene expression microarrays, the transcriptomic profiles of influenza-infected patients with severe (N = 11), moderate (N = 40) and mild (N = 83) symptoms were compared with the febrile patients of unknown etiology (N = 73). We found that influenza-infected patients, regardless of their clinical outcomes, had a stronger induction of antiviral and cytokine responses and a stronger attenuation of NK and T cell responses in comparison with those with unknown etiology. More importantly, we found that both interferon and ubiquitination signaling were strongly attenuated in patients with the most severe outcomes in comparison with those with moderate and mild outcomes, suggesting the protective roles of these pathways in disease pathogenesis.The attenuation of interferon and ubiquitination pathways may associate with the clinical outcomes of influenza patients.

Measurement error of a phenotypic trait reduces the power to detect genetic associations. We examined the impact of sample size, allele frequency and effect size in presence of measurement error for quantitative traits. The statistical power to detect genetic association with phenotype mean and variability was investigated analytically. The non-centrality parameter for a non-central F distribution was derived and verified using computer simulations. We obtained equivalent formulas for the cost of phenotype measurement error. Effects of differences in measurements were examined in a genome-wide association study (GWAS) of two grading scales for cataract and a replication study of genetic variants influencing blood pressure. The mean absolute difference between the analytic power and simulation power for comparison of phenotypic means and variances was less than 0.005, and the absolute difference did not exceed 0.02. To maintain the same power, a one standard deviation (SD) in measurement error of a standard normal distributed trait required a one-fold increase in sample size for comparison of means, and a three-fold increase in sample size for comparison of variances. GWAS results revealed almost no overlap in the significant SNPs (p

The diagnosis of childhood neurological disorders remains challenging given the overlapping clinical presentation across subgroups and heterogeneous presentation within subgroups. To determine the underlying genetic cause of a severe neurological disorder in a large consanguineous Pakistani family presenting with severe scoliosis, anarthria and progressive neuromuscular degeneration, we performed genome-wide homozygosity mapping accompanied by whole-exome sequencing in two affected first cousins and their unaffected parents to find the causative mutation. We identified a novel homozygous splice-site mutation (c.3512+1G>A) in the ALS2 gene (NM_020919.3) encoding alsin that segregated with the disease in this family. Homozygous loss-of-function mutations in ALS2 are known to cause juvenile-onset amyotrophic lateral sclerosis (ALS), one of the many neurological conditions having overlapping symptoms with many neurological phenotypes. RT-PCR validation revealed that the mutation resulted in exon-skipping as well as the use of an alternative donor splice, both of which are predicted to cause loss-of-function of the resulting proteins. By examining 216 known neurological disease genes in our exome sequencing data, we also identified 9 other rare nonsynonymous mutations in these genes, some of which lie in highly conserved regions. Sequencing of a single proband might have led to mis-identification of some of these as the causative variant. Our findings established a firm diagnosis of juvenile ALS in this family, thus demonstrating the use of whole exome sequencing combined with linkage analysis in families as a powerful tool for establishing a quick and precise genetic diagnosis of complex neurological phenotypes.

The BLK and CD40 loci have been associated with Kawasaki disease (KD) in two genome-wide association studies (GWAS) conducted in a Taiwanese population of Han Chinese ancestry (Taiwanese) and in Japanese cohorts. Here we build on these findings with replication studies of the BLK and CD40 loci in populations of Korean and European descent. The BLK region was significantly associated with KD susceptibility in both populations. Within the BLK gene the rs2736340-located linkage disequilibrium (LD ) comprising the promoter and first intron was strongly associated with KD, with the combined results of Asian studies including Taiwanese, Japanese, and Korean populations (2,539 KD patients and 7,021 controls) providing very compelling evidence of association (rs2736340, OR = 1.498, 1.354-1.657; P = 4.74×10(-31)). We determined the percentage of B cells present in the peripheral blood mononuclear cell (PBMC) population and the expression of BLK in the peripheral blood leukocytes (leukocytes) of KD patients during the acute and convalescent stages. The percentage of B cells in the PBMC population and the expression of BLK in leukocytes were induced in patients in the acute stage of KD. In B cell lines derived from KD patients, and in purified B cells from KD patients obtained during the acute stage, those with the risk allele of rs2736340 expressed significantly lower levels of BLK. These results suggest that peripheral B cells play a pathogenic role during the acute stage of KD. Decreased BLK expression in peripheral blood B cells may alter B cell function and predispose individuals to KD. These associative data suggest a role for B cells during acute KD. Understanding the functional implications may facilitate the development of B cell-mediated therapy for KD.

A recent genome-wide association study (GWAS) identified susceptibility loci for dengue shock syndrome (DSS) at MICB rs3132468 and PLCE1 rs3740360. The aim of this study was to define the extent to which MICB (rs3132468) and PLCE1 (rs3740360) were associated with less severe clinical phenotypes of pediatric and adult dengue.3961 laboratory-confirmed dengue cases and 5968 controls were genotyped at MICB rs3132468 and PLCE1 rs3740360. Per-allele odds ratios (OR) with 95% confidence intervals (CI) were calculated for each patient cohort. Pooled analyses were performed for adults and paediatrics respectively using a fixed effects model.Pooled analysis of the paediatric and adult cohorts indicated a significant association between MICB rs3132468 and dengue cases without shock (OR = 1.15; 95%CI: 1.07 - 1.24; P = 0.0012). Similarly, pooled analysis of pediatric and adult cohorts indicated a significant association between dengue cases without shock and PLCE1 rs3740360 (OR = 0.92; 95%CI: 0.85 - 0.99; P = 0.018). We also note significant association between both SNPs (OR = 1.48; P = 0.0075 for MICB rs3132468 and OR = 0.75, P = 0.041 for PLCE1 rs3740360) and dengue in infants.This study confirms that the MICB rs3132468 and PLCE1 rs3740360 risk genotypes are not only associated with DSS, but are also associated with less severe clinical phenotypes of dengue, as well as with dengue in infants. These findings have implications for our understanding of dengue pathogenesis.

Tuberculosis (TB), caused by Mycobacterium tuberculosis (M.tb), remains the leading cause of mortality from a single infectious agent. Each year around 9 million individuals newly develop active TB disease, and over 2 billion individuals are latently infected with M.tb worldwide, thus being at risk of developing TB reactivation disease later in life. The underlying mechanisms and pathways of protection against TB in humans, as well as the dynamics of the host response to M.tb infection, are incompletely understood. We carried out whole-genome expression profiling on a cohort of TB patients longitudinally sampled along 3 time-points: during active infection, during treatment, and after completion of curative treatment. We identified molecular signatures involving the upregulation of type-1 interferon (α/β) mediated signaling and chronic inflammation during active TB disease in an Indonesian population, in line with results from two recent studies in ethnically and epidemiologically different populations in Europe and South Africa. Expression profiles were captured in neutrophil-depleted blood samples, indicating a major contribution of lymphocytes and myeloid cells. Expression of type-1 interferon (α/β) genes mediated was also upregulated in the lungs of M.tb infected mice and in infected human macrophages. In patients, the regulated gene expression-signature normalized during treatment, including the type-1 interferon mediated signaling and a concurrent opposite regulation of interferon-gamma. Further analysis revealed IL15RA, UBE2L6 and GBP4 as molecules involved in the type-I interferon response in all three experimental models. Our data is highly suggestive that the innate immune type-I interferon signaling cascade could be used as a quantitative tool for monitoring active TB disease, and provide evidence that components of the patient's blood gene expression signature bear similarities to the pulmonary and macrophage response to mycobacterial infection.