TMEM132: an ancient architecture of cohesin and immunoglobulin domains define a new family of neural adhesion molecules

Genetic analysis of a Kaijiang duck conservation population through genome-wide scan

1. The Kaijiang duck is a native Chinese breed known for its excellent egg laying performance, killing-out percentage (88.57%), and disease resistance. The assessment of population genetic structure is the basis for understanding the genetics of indigenous breeds and for their protection and management.2. In this study, whole-genome sequencing was performed on 60 Kaijiang ducks to identify genetic variations and investigate the population structure. Homozygosity (ROH) analysis was conducted to assess inbreeding levels in the population.3. The study revealed a moderate level of inbreeding, indicated by an average inbreeding coefficient of 0.1043. This may impact the overall genetic diversity.4. Genomic Regions of Interest identified included 168 genomic regions exhibiting high levels of autozygosity. These regions were associated with processes including muscle growth, pigmentation, neuromodulation, and growth and reproduction.5. The significance of these pathways indicated their potential role in shaping the desirable traits of the Kaijiang duck. These findings provide insights into the genetic basis of the Kaijiang duck's desirable traits and can inform future breeding and conservation efforts.

Multidimensional analysis of TMEM132A in pan-cancer: unveiling its potential as a biomarker for treatment response prediction

Background: TMEM132A is a transmembrane protein that regulates gastric cancer cell malignancy and overall survival in bladder cancer patients. However, while some studies have investigated the involvement of TMEM132A in specific cancers, further systematic studies are required to elucidate its specific mechanisms of action in different cancer types. Methods: We investigated the pan-cancer role of TMEM132A using several databases. We analyzed TMEM132A expression and its correlation with clinical survival, immune checkpoints, tumor stemness score, prognostic value, immunomodulators, genomic profiles, immunological characteristics, immunotherapy and functional enrichment. Results: First, it was observed that TMEM132A expression levels were higher in the majority of tumors compared to non-tumor tissues. In addition, high TMEM132A expression may have a higher prognostic value in some cancers. Furthermore, TMEM132A was significantly associated with immune checkpoints, immunomodulators, prognosis, immunomodulatory genes, tumor stemness score, cell function status and immune infiltration in most tumors. Further analysis of TMEM132A-related gene enrichment, mutation sites and types, RNA modification and genomic heterogeneity showed that the major mutations of TMEM132A were missense mutations and that TMEM132A plays a very important role in UCEC, LUAD and LIHC. Finally, these results suggest that high TMEM132A expression may be associated with a better response to specific immunotherapies. Conclusion: This comprehensive study uncovers an important function for TMEM132A in different types of cancer. It also has the potential to identify TMEM132A as a potential biomarker for predicting treatment response. This may help us to better understand how TMEM132A plays a role in cancer and provide valuable insights for developing personalised treatments.

SMARCA4 Loss and Mutated β-Catenin Induce Proliferative Lesions in the Murine Embryonic Cerebellum

Almost all medulloblastomas (MB) of the Wingless/Int-1 (WNT) type are characterized by hotspot mutations in CTNNB1, and mouse models have convincingly demonstrated the tumor-initiating role of these mutations. Additional alterations in SMARCA4 are detected in ∼20% of WNT MB, but their functional role is mostly unknown. We, therefore, amended previously described brain lipid binding protein (Blbp)-cre::Ctnnb1(ex3)fl/wt mice by the introduction of floxed Smarca4 alleles. Unexpectedly, mutated and thereby stabilized β-catenin on its own induced severe developmental phenotypes in male and female Blbp-cre::Ctnnb1(ex3)fl/wt mice in our hands, including a thinned cerebral cortex, hydrocephalus, missing cerebellar layering, and cell accumulations in the brainstem and cerebellum. An additional loss of SMARCA4 even resulted in prenatal death for most mice. Respective Blbp-cre::Ctnnb1(ex3)fl/wt::Smarca4fl/rec mutants (male and female) developed large proliferative lesions in the cerebellum evolving from E13.5 to E16.5. Histological and molecular analysis of these lesions by DNA methylation profiling and single-cell RNA sequencing suggested an origin in early undifferentiated SOX2-positive cerebellar progenitors. Furthermore, upregulated WNT signaling, altered actin/cytoskeleton organization, and reduced neuronal differentiation were evident in mutant cells. In vitro, cells harboring alterations in both Ctnnb1 and Smarca4 were negatively selected and did not show tumorigenic potential after transplantation in adult female recipient mice. However, in cerebellar explant cultures, mutant cells displayed significantly increased proliferation, suggesting an important role of the embryonic microenvironment in the development of lesions. Altogether, these results represent an important first step toward the unraveling of tumorigenic mechanisms induced by aberrant WNT signaling and SMARCA4 deficiency.

Heritability estimates and genome-wide association study of methane emission traits in Nellore cattle

The objectives of the present study were to estimate the heritability for daily methane emission (CH4) and residual daily methane emission (CH4res) in Nellore cattle, as well as to perform genome-wide association studies (GWAS) to identify genomic regions and candidate genes influencing the genetic variation of CH4 and CH4res. Methane emission phenotypes of 743 Nellore animals belonging to 3 breeding programs were evaluated. CH4 was measured using the sulfur hexafluoride (SF6) tracer technique (which involves an SF6 permeation tube introduced into the rumen, and an appropriate apparatus on each animal), and CH4res was obtained as the difference between observed CH4 and CH4 adjusted for dry matter intake. A total of 6,252 genotyped individuals were used for genomic analyses. Data were analyzed with a univariate animal model by the single-step GBLUP method using the average information restricted maximum likelihood (AIREML) algorithm. The effects of single nucleotide polymorphisms (SNPs) were obtained using a single-step GWAS approach. Candidate genes were identified based on genomic windows associated with quantitative trait loci (QTLs) related to the 2 traits. Annotation of QTLs and identification of candidate genes were based on the initial and final coordinates of each genomic window considering the bovine genome ARS-UCD1.2 assembly. Heritability estimates were of moderate to high magnitude, being 0.42 ± 0.09 for CH4 and 0.21 ± 0.09 for CH4res, indicating that these traits will respond rapidly to genetic selection. GWAS revealed 11 and 15 SNPs that were significantly associated (P < 10-6) with genetic variation of CH4 and CH4res, respectively. QTLs associated with feed efficiency, residual feed intake, body weight, and height overlapped with significant markers for the traits evaluated. Ten candidate genes were present in the regions of significant SNPs; 3 were associated with CH4 and 7 with CH4res. The identified genes are related to different functions such as modulation of the rumen microbiota, fatty acid production, and lipid metabolism. CH4 and CH4res presented sufficient genetic variation and may respond rapidly to selection. Therefore, these traits can be included in animal breeding programs aimed at reducing enteric methane emissions across generations.

TMEM132D and VIPR2 Polymorphisms as Genetic Risk Loci for Retinal Detachment: A Genome-Wide Association Study and Polygenic Risk Score Analysis

Purpose

Retinal detachment (RD) is a sight-threatening ocular disease caused by separation of the neurosensory retina from the underlying retinal pigment epithelium layer. Its genetic basis is unclear because of a limited amount of data. In this study, we aimed to identify genetic risk loci associated with RD in participants without diabetes mellitus and to construct a polygenic risk score (PRS) to predict the risk of RD.

Methods

A genome-wide association study was conducted using data from the Taiwan Biobank to identify RD risk loci. A total of 1533 RD cases and 106,270 controls were recruited, all of whom were Han Chinese. Replication studies were performed using data from the UK Biobank and Biobank Japan. To construct the PRS, a traditional clumping and thresholding method was performed and validated by fivefold cross-validation.

Results

Two novel loci with significant associations were identified. These two genes were TMEM132D (lead single nucleotide polymorphism [SNP]: rs264498, adjusted-P = 7.18 × 10-9) and VIPR2 (lead SNP: rs3812305, adjusted-P = 8.38 × 10-9). The developed PRS was effective in discriminating individuals at high risk of RD with a dose-response relationship. The quartile with the highest risk had an odds ratio of 1244.748 compared to the lowest risk group (95% confidence interval, 175.174-8844.892).

Conclusions

TMEM132D and VIPR2 polymorphisms are genetic candidates linked to RD in Han Chinese populations. Our proposed PRS was effective at discriminating high-risk from low-risk individuals.

Genetic Architecture of Ischaemic Strokes after COVID-19 Shows Similarities with Large Vessel Strokes

We aimed to analyse whether patients with ischaemic stroke (IS) occurring within eight days after the onset of COVID-19 (IS-COV) are associated with a specific aetiology of IS. We used SUPERGNOVA to identify genome regions that correlate between the IS-COV cohort (73 IS-COV cases vs. 701 population controls) and different aetiological subtypes. Polygenic risk scores (PRSs) for each subtype were generated and tested in the IS-COV cohort using PRSice-2 and PLINK to find genetic associations. Both analyses used the IS-COV cohort and GWAS from MEGASTROKE (67,162 stroke patients vs. 454,450 population controls), GIGASTROKE (110,182 vs. 1,503,898), and the NINDS Stroke Genetics Network (16,851 vs. 32,473). Three genomic regions were associated (p-value < 0.05) with large artery atherosclerosis (LAA) and cardioembolic stroke (CES). We found four loci targeting the genes PITX2 (rs10033464, IS-COV beta = 0.04, p-value = 2.3 × 10-2, se = 0.02), previously associated with CES, HS6ST1 (rs4662630, IS-COV beta = -0.04, p-value = 1.3 × 10-3, se = 0.01), TMEM132E (rs12941838 IS-COV beta = 0.05, p-value = 3.6 × 10-4, se = 0.01), and RFFL (rs797989 IS-COV beta = 0.03, p-value = 1.0 × 10-2, se = 0.01). A statistically significant PRS was observed for LAA. Our results suggest that IS-COV cases are genetically similar to LAA and CES subtypes. Larger cohorts are needed to assess if the genetic factors in IS-COV cases are shared with the general population or specific to viral infection.

TMEM132A regulates mouse hindgut morphogenesis and caudal development

Caudal developmental defects, including caudal regression, caudal dysgenesis and sirenomelia, are devastating conditions affecting the skeletal, nervous, digestive, reproductive and excretory systems. Defects in mesodermal migration and blood supply to the caudal region have been identified as possible causes of caudal developmental defects, but neither satisfactorily explains the structural malformations in all three germ layers. Here, we describe caudal developmental defects in transmembrane protein 132a (Tmem132a) mutant mice, including skeletal, posterior neural tube closure, genitourinary tract and hindgut defects. We show that, in Tmem132a mutant embryos, visceral endoderm fails to be excluded from the medial region of early hindgut, leading directly to the loss or malformation of cloaca-derived genitourinary and gastrointestinal structures, and indirectly to the neural tube and kidney/ureter defects. We find that TMEM132A mediates intercellular interaction, and physically interacts with planar cell polarity (PCP) regulators CELSR1 and FZD6. Genetically, Tmem132a regulates neural tube closure synergistically with another PCP regulator Vangl2. In summary, we have identified Tmem132a as a new regulator of PCP, and hindgut malformation as the underlying cause of developmental defects in multiple caudal structures.

Proteogenomic analysis of human cerebrospinal fluid identifies neurologically relevant regulation and informs causal proteins for Alzheimer's disease

The integration of quantitative trait loci (QTL) with disease genome-wide association studies (GWAS) has proven successful at prioritizing candidate genes at disease-associated loci. QTL mapping has mainly been focused on multi-tissue expression QTL or plasma protein QTL (pQTL). Here we generated the largest-to-date cerebrospinal fluid (CSF) pQTL atlas by analyzing 7,028 proteins in 3,107 samples. We identified 3,373 independent study-wide associations for 1,961 proteins, including 2,448 novel pQTLs of which 1,585 are unique to CSF, demonstrating unique genetic regulation of the CSF proteome. In addition to the established chr6p22.2-21.32 HLA region, we identified pleiotropic regions on chr3q28 near OSTN and chr19q13.32 near APOE that were enriched for neuron-specificity and neurological development. We also integrated this pQTL atlas with the latest Alzheimer's disease (AD) GWAS through PWAS, colocalization and Mendelian Randomization and identified 42 putative causal proteins for AD, 15 of which have drugs available. Finally, we developed a proteomics-based risk score for AD that outperforms genetics-based polygenic risk scores. These findings will be instrumental to further understand the biology and identify causal and druggable proteins for brain and neurological traits.

DNA Methylation and Prospects for Predicting the Therapeutic Effect of Neoadjuvant Chemotherapy for Triple-Negative and Luminal B Breast Cancer

Despite advances in the diagnosis and treatment of breast cancer (BC), the main cause of deaths is resistance to existing therapies. An approach to improve the effectiveness of therapy in patients with aggressive BC subtypes is neoadjuvant chemotherapy (NACT). Yet, the response to NACT for aggressive subtypes is less than 65% according to large clinical trials. An obvious fact is the lack of biomarkers predicting the therapeutic effect of NACT. In a search for epigenetic markers, we performed genome-wide differential methylation screening by XmaI-RRBS in cohorts of NACT responders and nonresponders, for triple-negative (TN) and luminal B tumors. The predictive potential of the most discriminative loci was further assessed in independent cohorts by methylation-sensitive restriction enzyme quantitative PCR (MSRE-qPCR), a promising method for the implementation of DNA methylation markers in diagnostic laboratories. The selected most informative individual markers were combined into panels demonstrating cvAUC = 0.83 (TMEM132D and MYO15B markers panel) for TN tumors and cvAUC = 0.76 (TTC34, LTBR and CLEC14A) for luminal B tumors. The combination of methylation markers with clinical features that correlate with NACT effect (clinical stage for TN and lymph node status for luminal B tumors) produces better classifiers, with cvAUC = 0.87 for TN tumors and cvAUC = 0.83 for luminal B tumors. Thus, clinical characteristics predictive of NACT response are independently additive to the epigenetic classifier and in combination improve prediction.

GRP75 as a functional element of cholix transcytosis

Cholix (Chx) is secreted by non-pandemic strains of Vibrio cholerae in the intestinal lumen. For this exotoxin to induce cell death in non-polarized cells in the intestinal lamina propria, it must traverse the epithelium in the fully intact form. We identified host cell elements in polarized enterocytes associated with Chx endocytosis and apical to basal (A→B) vesicular transcytosis. This pathway overcomes endogenous mechanisms of apical vesicle recycling and lysosomal targeting by interacting with several host cell proteins that include the 75 kDa glucose-regulated protein (GRP75). Apical endocytosis of Chx appears to involve the single membrane spanning protein TMEM132A, and interaction with furin before it engages GRP75 in apical vesicular structures. Sorting within these apical vesicles results in Chx being trafficked to the basal region of cells in association with the Lectin, Mannose Binding 1 protein LMAN1. In this location, Chx interacts with the basement membrane-specific heparan sulfate proteoglycan perlecan in recycling endosomes prior to its release from this basal vesicular compartment to enter the underlying lamina propria. While the furin and LMAN1 elements of this Chx transcytosis pathway undergo cellular redistribution that are reflective of the polarity shifts noted for coatamer complexes COPI and COPII, GRP75 and perlecan fail to show these dramatic rearrangements. Together, these data define essential steps in the A→B transcytosis pathway accessed by Chx to reach the intestinal lamina propria where it can engage and intoxicate certain non-polarized cells.

Protein disulfide isomerase A3 activity promotes extracellular accumulation of proteins relevant to basal breast cancer outcomes in human MDA-MB-A231 breast cancer cells

Metastasis and recurrence of breast cancer remain major causes of patient mortality, and there is an ongoing need to identify new therapeutic targets relevant to tumor invasion. Protein disulfide isomerase A3 (PDIA3) is a disulfide oxidoreductase and isomerase of the endoplasmic reticulum that has known extracellular substrates and has been correlated with aggressive breast cancers. We show that either prior PDIA3 inhibition by the disulfide isomerase inhibitor 16F16 or depletion of heparin-binding proteins strongly reduces the activity of conditioned medium (CM) of MDA-MB-231 human breast cancer cells to support promigratory cell spreading and F-actin organization by newly adherent MDA-MB-231 cells. Quantitative proteomics to investigate effects of 16F16 inhibition on heparin-binding proteins in the CM of MDA-MB-231 cells identified 80 proteins reproducibly decreased at least twofold (at q ≤ 0.05) after 16F16 treatment. By Gene Ontology analysis, many of these have roles in extracellular matrix (ECM) structure and function and cell adhesion; ribosomal proteins that also correlate with extracellular vesicles were also identified. Protein-protein interaction analysis showed that many of the extracellular proteins have known network interactions with each other. The predominant types of disulfide-bonded domains in the extracellular proteins contained β-hairpin folds, with the knottin fold the most common. From human breast cancer data sets, the extracellular proteins were found to correlate specifically with the basal subtype of breast cancer and their high expression in tumors correlated with reduced distant metastasis-free survival. These data provide new evidence that PDIA3 may be a relevant therapeutic target to alter properties of the ECM-associated microenvironment in basal breast cancer.

TMEM132A ensures mouse caudal neural tube closure and regulates integrin-based mesodermal migration

Coordinated migration of the mesoderm is essential for accurate organization of the body plan during embryogenesis. However, little is known about how mesoderm migration influences posterior neural tube closure in mammals. Here, we show that spinal neural tube closure and lateral migration of the caudal paraxial mesoderm depend on transmembrane protein 132A (TMEM132A), a single-pass type I transmembrane protein, the function of which is not fully understood. Our study in Tmem132a-null mice and cell models demonstrates that TMEM132A regulates several integrins and downstream integrin pathway activation as well as cell migration behaviors. Our data also implicates mesoderm migration in elevation of the caudal neural folds and successful closure of the caudal neural tube. These results suggest a requirement for paraxial mesodermal cell migration during spinal neural tube closure, disruption of which may lead to spina bifida.

Selecting Genetic Variants and Interactions Associated with Amyotrophic Lateral Sclerosis: A Group LASSO Approach

Amyotrophic lateral sclerosis (ALS) is a multi-system neurodegenerative disease that affects both upper and lower motor neurons, resulting from a combination of genetic, environmental, and lifestyle factors. Usually, the association between single-nucleotide polymorphisms (SNPs) and this disease is tested individually, which leads to the testing of multiple hypotheses. In addition, this classical approach does not support the detection of interaction-dependent SNPs. We applied a two-step procedure to select SNPs and pairwise interactions associated with ALS. SNP data from 276 ALS patients and 268 controls were analyzed by a two-step group LASSO in 2000 iterations. In the first step, we fitted a group LASSO model to a bootstrap sample and a random subset of predictors (25%) from the original data set aiming to screen for important SNPs and, in the second step, we fitted a hierarchical group LASSO model to evaluate pairwise interactions. An in silico analysis was performed on a set of variables, which were prioritized according to their bootstrap selection frequency. We identified seven SNPs (rs16984239, rs10459680, rs1436918, rs1037666, rs4552942, rs10773543, and rs2241493) and two pairwise interactions (rs16984239:rs2118657 and rs16984239:rs3172469) potentially involved in nervous system conservation and function. These results may contribute to the understanding of ALS pathogenesis, its diagnosis, and therapeutic strategy improvement.

Differential expression of the Tmem132 family genes in the developing mouse nervous system

The family of novel transmembrane proteins (TMEM) 132 have been associated with multiple neurological disorders and cancers in humans, but have hardly been studied in vivo. Here we report the expression patterns of the five Tmem132 genes (a, b, c, d and e) in developing mouse nervous system with RNA in situ hybridization in wholemount embryos and tissue sections. Our results reveal differential and partially overlapping expression of multiple Tmem132 family members in both the central and peripheral nervous system, suggesting potential partial redundancy among them.

Finding New Ways How to Control BACE1

Recently, all applications of BACE1 inhibitors failed as therapeutical targets for Alzheimer´s disease (AD) due to severe side effects. Therefore, alternative ways for treatment development are a hot research topic. The present analysis investigates BACE1 protein-protein interaction networks and attempts to solve the absence of complete knowledge about pathways involving BACE1. A bioinformatics analysis matched the functions of the non-substrate interaction network with Voltage-gated potassium channels, which also appear as top priority protein nodes. Targeting BACE1 interactions with PS1 and GGA-s, blocking of BACE1 access to APP by BRI3 and RTN-s, activation of Wnt signaling and upregulation of β-catenin, and brain delivery of the extracellular domain of p75^NTR, are the main alternatives to the use of BACE 1 inhibitors highlighted by the analysis. The pathway enrichment analysis also emphasized substrates and substrate candidates with essential biological functions, which cleavage must remain controlled. They include ephrin receptors, ROBO1, ROBO2, CNTN-s, CASPR-s, CD147, CypB, TTR, APLP1/APLP2, NRXN-s, and PTPR-s. The analysis of the interaction subnetwork of BACE1 functionally related to inflammation identified a connection to three cardiomyopathies, which supports the hypothesis of the common molecular mechanisms with AD. A lot of potential shows the regulation of BACE1 activity through post-translational modifications. The interaction network of BACE1 and its phosphorylation enzyme CSNK1D functionally match the Circadian clock, p53, and Hedgehog signaling pathways. The regulation of BACE1 glycosylation could be achieved through N-acetylglucosamine transferases, α-(1→6)-fucosyltransferase, β-galactoside α-(2→6)-sialyltransferases, galactosyltransferases, and mannosidases suggested by the interaction network analysis of BACE1-MGAT3. The present analysis proposes possibilities for the alternative control of AD pathology.

Transcriptomic characterization of tissues from patients and subsequent pathway analyses reveal biological pathways that are implicated in spastic ataxia

BACKGROUND

Spastic ataxias (SAs) encompass a group of rare and severe neurodegenerative diseases, characterized by an overlap between ataxia and spastic paraplegia clinical features. They have been associated with pathogenic variants in a number of genes, including GBA2. This gene codes for the non-lysososomal β-glucosylceramidase, which is involved in sphingolipid metabolism through its catalytic role in the degradation of glucosylceramide. However, the mechanism by which GBA2 variants lead to the development of SA is still unclear.

METHODS

In this work, we perform next-generation RNA-sequencing (RNA-seq), in an attempt to discover differentially expressed genes (DEGs) in lymphoblastoid, fibroblast cell lines and induced pluripotent stem cell-derived neurons derived from patients with SA, homozygous for the GBA2 c.1780G > C missense variant. We further exploit DEGs in pathway analyses in order to elucidate candidate molecular mechanisms that are implicated in the development of the GBA2 gene-associated SA.

RESULTS

Our data reveal a total of 5217 genes with significantly altered expression between patient and control tested tissues. Furthermore, the most significant extracted pathways are presented and discussed for their possible role in the pathogenesis of the disease. Among them are the oxidative stress, neuroinflammation, sphingolipid signaling and metabolism, PI3K-Akt and MAPK signaling pathways.

CONCLUSIONS

Overall, our work examines for the first time the transcriptome profiles of GBA2-associated SA patients and suggests pathways and pathway synergies that could possibly have a role in SA pathogenesis. Lastly, it provides a list of DEGs and pathways that could be further validated towards the discovery of disease biomarkers.

Roles, molecular mechanisms, and signaling pathways of TMEMs in neurological diseases

Transmembrane protein family members (TMEMs) span the entire lipid bilayer and act as channels that allow the transport of specific substances through biofilms. The functions of most TMEMs are unexplored. Numerous studies have shown that TMEMs are involved in the pathophysiological processes of various nervous system diseases, but the specific mechanisms of TMEMs in the pathogenesis of diseases remain unclear. In this review, we discuss the expression, physiological functions, and molecular mechanisms of TMEMs in brain tumors, psychiatric disorders, abnormal motor activity, cobblestone lissencephaly, neuropathic pain, traumatic brain injury, and other disorders of the nervous system. Additionally, we propose that TMEMs may be used as prognostic markers and potential therapeutic targets in patients with various neurological diseases.

17q12 Recurrent Deletions and Duplications: Description of a Case Series with Neuropsychiatric Phenotype

Syndromic neurodevelopmental disorders are usually investigated through genetics technologies, within which array comparative genomic hybridization (Array-CGH) is still considered the first-tier clinical diagnostic test. Among recurrent syndromic imbalances, 17q12 deletions and duplications are characterized by neurodevelopmental disorders associated with visceral developmental disorders, although expressive variability is common. Here we describe a case series of 12 patients with 17q12 chromosomal imbalances, in order to expand the phenotypic characterization of these recurrent syndromes whose diagnosis is often underestimated, especially if only mild traits are present. Gene content and genotype-phenotype correlations have been discussed, with special regard to neuropsychiatric features, whose impact often requires etiologic analysis.

Positive selection of skeleton-related genes during duck domestication revealed by whole genome sequencing

BACKGROUND

Domestication alters several phenotypic, neurological, and physiological traits in domestic animals compared to those in their wild ancestors. Domestic ducks originated from mallards, and some studies have shown that spot-billed ducks may have also made minor genetic contributions to domestication. Compared with the two ancestral species, domestic ducks generally differ in body size and bone morphology. In this study, we performed both genomic and transcriptomic analyses to identify candidate genes for elucidating the genetic mechanisms underlying phenotypic variation.

METHODS

In this study, the duck genome data from eight domestic breeds and two wild species were collected to study the genetic changes during domestication. And the transcriptome data of different tissues from wild ducks and seven domestic ducks were used to reveal the expression difference between wild and domestic ducks.

RESULTS

Using fixation index (Fst) algorithm and transcriptome data, we found that the genes related to skeletal development had high Fst values in wild and domestic breeds, and the differentially expressed genes were mainly enriched in the ossification pathway. Our data strongly suggest that the skeletal systems of domestic ducks were changed to adapt to artificial selection for larger sizes. In addition, by combining the genome and transcriptome data, we found that some Fst candidate genes exhibited different expression patterns, and these genes were found to be involved in digestive, immune, and metabolic functions.

CONCLUSIONS

A wide range of phenotypic differences exists between domestic and wild ducks. Through both genome and transcriptome analyses, we found that genes related to the skeletal system in domestic ducks were strongly selected. Our findings provide new insight into duck domestication and selection effects during the domestication.

HOXB8 Counteracts MAPK/ERK Oncogenic Signaling in a Chicken Embryo Model of Neoplasia

HOX transcription factors are members of an evolutionarily conserved family of proteins required for the establishment of the anteroposterior body axis during bilaterian development. Although they are often deregulated in cancers, the molecular mechanisms by which they act as oncogenes or tumor suppressor genes are only partially understood. Since the MAPK/ERK signaling pathway is deregulated in most cancers, we aimed at apprehending if and how the Hox proteins interact with ERK oncogenicity. Using an in vivo neoplasia model in the chicken embryo consisting in the overactivation of the ERK1/2 kinases in the trunk neural tube, we analyzed the consequences of the HOXB8 gain of function at the morphological and transcriptional levels. We found that HOXB8 acts as a tumor suppressor, counteracting ERK-induced neoplasia. The HOXB8 tumor suppressor function relies on a large reversion of the oncogenic transcriptome induced by ERK. In addition to showing that the HOXB8 protein controls the transcriptional responsiveness to ERK oncogenic signaling, our study identified new downstream targets of ERK oncogenic activation in an in vivo context that could provide clues for therapeutic strategies.

A Latent Gaussian Copula Model for Mixed Data Analysis in Brain Imaging Genetics

Recent advances in imaging genetics make it possible to combine different types of data including medical images like functional magnetic resonance imaging (fMRI) and genetic data like single nucleotide polymorphisms (SNPs) for comprehensive diagnosis of mental disorders. Understanding complex interactions among these heterogeneous data may give rise to a new perspective, while at the same time demand statistical models for their integration. Various graphical models have been proposed for the study of interaction or association networks with continuous, binary, and count data as well as the mixture of them. However, limited efforts have been made for the multinomial case, for instance, SNP data. Our goal is therefore to fill the void by developing a graphical model for the integration of fMRI image and SNP data, which can provide deeper understanding of the unknown neurogenetic mechanism. In this article, we propose a latent Gaussian copula model for mixed data containing multinomial components. We assume that the discrete variable is obtained by discretizing a latent (unobserved) continuous variable and then create a semi-rank based estimator of the graph structure. The simulation results demonstrate that the proposed latent correlation has more steady and accurate performance than several existing methods in detecting graph structure. When applying to a real schizophrenia data consisting of SNP array and fMRI image collected by the Mind Clinical Imaging Consortium (MCIC), the proposed method reveals a set of distinct SNP-brain associations, which are verified to be biologically significant. The proposed model is statistically promising in handling mixed types of data including multinomial components, which can find widespread applications. To promote reproducible research, the R code is available at https://github.com/Aiying0512/LGCM.

Expression Profile of New Marker Genes Involved in Differentiation of Canine Adipose-Derived Stem Cells into Osteoblasts

Next-generation sequencing (RNAseq) analysis of gene expression changes during the long-term in vitro culture and osteogenic differentiation of ASCs remains to be important, as the analysis provides important clues toward employing stem cells as a therapeutic intervention. In this study, the cells were isolated from adipose tissue obtained during routine surgical procedures and subjected to 14-day in vitro culture and differentiation. The mRNA transcript levels were evaluated using the Illumina platform, resulting in the detection of 19,856 gene transcripts. The most differentially expressed genes (fold change >|2|, adjusted p value < 0.05), between day 1, day 14 and differentiated cell cultures were extracted and subjected to bioinformatical analysis based on the R programming language. The results of this study provide molecular insight into the processes that occur during long-term in vitro culture and osteogenic differentiation of ASCs, allowing the re-evaluation of the roles of some genes in MSC progression towards a range of lineages. The results improve the knowledge of the molecular mechanisms associated with long-term in vitro culture and differentiation of ASCs, as well as providing a point of reference for potential in vivo and clinical studies regarding these cells' application in regenerative medicine.

Whole-genome resequencing provides insights into the population structure and domestication signatures of ducks in eastern China

Background

Duck is an ancient domesticated animal with high economic value, used for its meat, eggs, and feathers. However, the origin of indigenous Chinese ducks remains elusive. To address this question, we performed whole-genome resequencing to first explore the genetic relationship among variants of these domestic ducks with their potential wild ancestors in eastern China, as well as understand how the their genomes were shaped by different natural and artificial selective pressures.

Results

Here, we report the resequencing of 60 ducks from Chinese spot-billed ducks (Anas zonorhyncha), mallards (Anas platyrhnchos), Fenghua ducks, Shaoxing ducks, Shanma ducks and Cherry Valley Pekin ducks of eastern China (ten from each population) at an average effective sequencing depth of ~ 6× per individual. The results of population and demographic analysis revealed a deep phylogenetic split between wild (Chinese spot-billed ducks and mallards) and domestic ducks. By applying selective sweep analysis, we identified that several candidate genes, important pathways and GO categories associated with artificial selection were functionally related to cellular adhesion, type 2 diabetes, lipid metabolism, the cell cycle, liver cell proliferation, and muscle functioning in domestic ducks.

Conclusion

Genetic structure analysis showed a close genetic relationship of Chinese spot-billed ducks and mallards, which supported that Chinese spot-billed ducks contributed to the breeding of domestic ducks. During the long history of artificial selection, domestic ducks have developed a complex biological adaptation to captivity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07710-2.

Genome-Wide Association Studies Reveal Susceptibility Loci for Noninfectious Claw Lesions in Holstein Dairy Cattle

Sole ulcers (SUs) and white line disease (WLD) are two common noninfectious claw lesions (NICL) that arise due to a compromised horn production and are frequent causes of lameness in dairy cattle, imposing welfare and profitability concerns. Low to moderate heritability estimates of SU and WLD susceptibility indicate that genetic selection could reduce their prevalence. To identify the susceptibility loci for SU, WLD, SU and/or WLD, and any type of noninfectious claw lesion, genome-wide association studies (GWAS) were performed using generalized linear mixed model (GLMM) regression, chunk-based association testing (CBAT), and a random forest (RF) approach. Cows from five commercial dairies in California were classified as controls having no lameness records and ≥6 years old (n = 102) or cases having SU (n = 152), WLD (n = 117), SU and/or WLD (SU + WLD, n = 198), or any type of noninfectious claw lesion (n = 217). The top single nucleotide polymorphisms (SNPs) were defined as those passing the Bonferroni-corrected suggestive and significance thresholds in the GLMM analysis or those that a validated RF model considered important. Effects of the top SNPs were quantified using Bayesian estimation. Linkage disequilibrium (LD) blocks defined by the top SNPs were explored for candidate genes and previously identified, functionally relevant quantitative trait loci. The GLMM and CBAT approaches revealed the same regions of association on BTA8 for SU and BTA13 common to WLD, SU + WLD, and NICL. These SNPs had effects significantly different from zero, and the LD blocks they defined explained a significant amount of phenotypic variance for each dataset (6.1-8.1%, p < 0.05), indicating the small but notable contribution of these regions to susceptibility. These regions contained candidate genes involved in wound healing, skin lesions, bone growth and mineralization, adipose tissue, and keratinization. The LD block defined by the most significant SNP on BTA8 for SU included a SNP previously associated with SU. The RF models were overfitted, indicating that the SNP effects were very small, thereby preventing meaningful interpretation of SNPs and any downstream analyses. These findings suggested that variants associated with various physiological systems may contribute to susceptibility for NICL, demonstrating the complexity of genetic predisposition.

The C. elegans homolog of human panic-disorder risk gene TMEM132D orchestrates neuronal morphogenesis through the WAVE-regulatory complex

TMEM132D is a human gene identified with multiple risk alleles for panic disorders, anxiety and major depressive disorders. Defining a conserved family of transmembrane proteins, TMEM132D and its homologs are still of unknown molecular functions. By generating loss-of-function mutants of the sole TMEM132 ortholog in C. elegans, we identify abnormal morphologic phenotypes in the dopaminergic PDE neurons. Using a yeast two-hybrid screen, we find that NAP1 directly interacts with the cytoplasmic domain of human TMEM132D, and mutations in C. elegans tmem-132 that disrupt interaction with NAP1 cause similar morphologic defects in the PDE neurons. NAP1 is a component of the WAVE regulatory complex (WRC) that controls F-actin cytoskeletal dynamics. Decreasing activity of WRC rescues the PDE defects in tmem-132 mutants, whereas gain-of-function of TMEM132D in mammalian cells inhibits WRC, leading to decreased abundance of select WRC components, impaired actin nucleation and cell motility. We propose that metazoan TMEM132 family proteins play evolutionarily conserved roles in regulating NAP1 protein homologs to restrict inappropriate WRC activity, cytoskeletal and morphologic changes in the cell.

Association of CSF proteins with tau and amyloid β levels in asymptomatic 70-year-olds

Background

Increased knowledge of the evolution of molecular changes in neurodegenerative disorders such as Alzheimer’s disease (AD) is important for the understanding of disease pathophysiology and also crucial to be able to identify and validate disease biomarkers. While several biological changes that occur early in the disease development have already been recognized, the need for further characterization of the pathophysiological mechanisms behind AD still remains.

Methods

In this study, we investigated cerebrospinal fluid (CSF) levels of 104 proteins in 307 asymptomatic 70-year-olds from the H70 Gothenburg Birth Cohort Studies using a multiplexed antibody- and bead-based technology.

Results

The protein levels were first correlated with the core AD CSF biomarker concentrations of total tau, phospho-tau and amyloid beta (Aβ42) in all individuals. Sixty-three proteins showed significant correlations to either total tau, phospho-tau or Aβ42. Thereafter, individuals were divided based on CSF Aβ42/Aβ40 ratio and Clinical Dementia Rating (CDR) score to determine if early changes in pathology and cognition had an effect on the correlations. We compared the associations of the analysed proteins with CSF markers between groups and found 33 proteins displaying significantly different associations for amyloid-positive individuals and amyloid-negative individuals, as defined by the CSF Aβ42/Aβ40 ratio. No differences in the associations could be seen for individuals divided by CDR score.

Conclusions

We identified a series of transmembrane proteins, proteins associated with or anchored to the plasma membrane, and proteins involved in or connected to synaptic vesicle transport to be associated with CSF biomarkers of amyloid and tau pathology in AD. Further studies are needed to explore these proteins’ role in AD pathophysiology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13195-021-00789-5.

A Ten-N6-Methyladenosine (m6A)-Modified Gene Signature Based on a Risk Score System Predicts Patient Prognosis in Rectum Adenocarcinoma

OBJECTIVES

The study aims to analyze the expression of N6-methyladenosine (m6A)-modified genes in rectum adenocarcinoma (READ) and identify reliable prognostic biomarkers to predict the prognosis of READ.

MATERIALS AND METHODS

RNA sequence data of READ and corresponding clinical survival data were obtained from The Cancer Genome Atlas (TCGA) database. N6-methyladenosine (m6A)-modified genes in READ were downloaded from the "m6Avar" database. Differentially expressed m6A-modified genes in READ stratified by different clinicopathological characteristics were identified using the "limma" package in R. Protein-protein interaction (PPI) network and co-expression analysis of differentially expressed genes (DEGs) were performed using "STRING" and Cytoscape, respectively. Principal component analysis (PCA) was done using R. In addition, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were used to functionally annotate the differentially expressed genes in different subgroups. Univariate Cox regression analyses were conducted to identify the powerful independent prognostic factors in READ associated with overall survival (OS). A robust likelihood-based survival model was built using the "rbsurv" package to screen for survival-associated signature genes. The Support Vector Machine (SVM) was used to predict the prognosis of READ through the risk score of survival-associated signature genes. Correlation analysis were carried out using GraphPad prism 8.

RESULTS

We screened 974 differentially expressed m6A-modified genes among four types of READ samples. Two READ subgroups (group 1 and group 2) were identified by K means clustering according to the expression of DEGs. The two subgroups were significantly different in overall survival and pathological stages. Next, 118 differentially expressed genes between the two subgroups were screened and the expression of 112 genes was found to be related to the prognosis of READ. Next, a panel of 10 survival-associated signature genes including adamtsl1, csmd2, fam13c, fam184a, klhl4, olfml2b, pdzd4, sec14l5, setbp1, tmem132b was constructed. The signature performed very well for prognosis prediction, time-dependent receiver-operating characteristic (ROC) analysis displaying an area under the curve (AUC) of 0.863, 0.8721, and 0.8752 for 3-year survival rate, prognostic status, and pathological stage prediction, respectively. Correlation analysis showed that the expression levels of the 10 m6A-modified genes were positively correlated with that of m6A demethylase FTO and ALKBH5.

CONCLUSION

This study identified potential m6A-modified genes that may be involved in the pathophysiology of READ and constructed a novel gene expression panel for READ risk stratification and prognosis prediction.

Medication for Acromegaly Reduces Expression of MUC16, MACC1 and GRHL2 in Pituitary Neuroendocrine Tumour Tissue

Acromegaly is a disease mainly caused by pituitary neuroendocrine tumor (PitNET) overproducing growth hormone. First-line medication for this condition is the use of somatostatin analogs (SSAs), that decrease tumor mass and induce antiproliferative effects on PitNET cells. Dopamine agonists (DAs) can also be used if SSA treatment is not effective. This study aimed to determine differences in transcriptome signatures induced by SSA/DA therapy in PitNET tissue. We selected tumor tissue from twelve patients with somatotropinomas, with half of the patients receiving SSA/DA treatment before surgery and the other half treatment naive. Transcriptome sequencing was then carried out to identify differentially expressed genes (DEGs) and their protein–protein interactions, using pathway analyses. We found 34 upregulated and six downregulated DEGs in patients with SSA/DA treatment. Three tumor development promoting factors MUC16, MACC1, and GRHL2, were significantly downregulated in therapy administered PitNET tissue; this finding was supported by functional studies in GH3 cells. Protein–protein interactions and pathway analyses revealed extracellular matrix involvement in the antiproliferative effects of this type of the drug treatment, with pronounced alterations in collagen regulation. Here, we have demonstrated that somatotropinomas can be distinguished based on their transcriptional profiles following SSA/DA therapy, and SSA/DA treatment does indeed cause changes in gene expression. Treatment with SSA/DA significantly downregulated several factors involved in tumorigenesis, including MUC16, MACC1, and GRHL2. Genes that were upregulated, however, did not have a direct influence on antiproliferative function in the PitNET cells. These findings suggested that SSA/DA treatment acted in a tumor suppressive manner and furthermore, collagen related interactions and pathways were enriched, implicating extracellular matrix involvement in this anti-tumor effect of drug treatment.

Genetic architecture underlying changes in carotenoid accumulation during the evolution of the blind Mexican cavefish, Astyanax mexicanus

Carotenoids are lipid-soluble yellow to orange pigments produced by plants, bacteria, and fungi. They are consumed by animals and metabolized to produce molecules essential for gene regulation, vision, and pigmentation. Cave animals represent an interesting opportunity to understand how carotenoid utilization evolves. Caves are devoid of light, eliminating primary production of energy through photosynthesis and, therefore, limiting carotenoid availability. Moreover, the selective pressures that favor carotenoid-based traits, like pigmentation and vision, are relaxed. Astyanax mexicanus is a species of fish with multiple river-adapted (surface) and cave-adapted populations (i.e., Tinaja, Pachón, Molino). Cavefish exhibit regressive features, such as loss of eyes and melanin pigment, and constructive traits, like increased sensory neuromasts and starvation resistance. Here, we show that, unlike surface fish, Tinaja and Pachón cavefish accumulate carotenoids in the visceral adipose tissue. Carotenoid accumulation is not observed in Molino cavefish, indicating that it is not an obligatory consequence of eye loss. We used quantitative trait loci mapping and RNA sequencing to investigate genetic changes associated with carotenoid accumulation. Our findings suggest that multiple stages of carotenoid processing may be altered in cavefish, including absorption and transport of lipids, cleavage of carotenoids into unpigmented molecules, and differential development of intestinal cell types involved in carotenoid assimilation. Our study establishes A. mexicanus as a model to study the genetic basis of natural variation in carotenoid accumulation and how it impacts physiology.

Hair cell α9α10 nicotinic acetylcholine receptor functional expression regulated by ligand binding and deafness gene products

Auditory hair cells receive olivocochlear efferent innervation, which refines tonotopic mapping, improves sound discrimination, and mitigates acoustic trauma. The olivocochlear synapse involves α9α10 nicotinic acetylcholine receptors (nAChRs), which assemble in hair cells only coincident with cholinergic innervation and do not express in recombinant mammalian cell lines. Here, genome-wide screening determined that assembly and surface expression of α9α10 require ligand binding. Ion channel function additionally demands an auxiliary subunit, which can be transmembrane inner ear (TMIE) or TMEM132e. Both of these single-pass transmembrane proteins are enriched in hair cells and underlie nonsyndromic human deafness. Inner hair cells from TMIE mutant mice show altered postsynaptic α9α10 function and retain α9α10-mediated transmission beyond the second postnatal week associated with abnormally persistent cholinergic innervation. Collectively, this study provides a mechanism to link cholinergic input with α9α10 assembly, identifies unexpected functions for human deafness genes TMIE/TMEM132e, and enables drug discovery for this elusive nAChR implicated in prevalent auditory disorders.

Altered levels of CSF proteins in patients with FTD, presymptomatic mutation carriers and non-carriers

BACKGROUND

The clinical presentations of frontotemporal dementia (FTD) are diverse and overlap with other neurological disorders. There are, as of today, no biomarkers in clinical practice for diagnosing the disorders. Here, we aimed to find protein markers in cerebrospinal fluid (CSF) from patients with FTD, presymptomatic mutation carriers and non-carriers.

METHODS

Antibody suspension bead arrays were used to analyse 328 proteins in CSF from patients with behavioural variant FTD (bvFTD, n = 16) and progressive primary aphasia (PPA, n = 13), as well as presymptomatic mutation carriers (PMC, n = 16) and non-carriers (NC, n = 8). A total of 492 antibodies were used to measure protein levels by direct labelling of the CSF samples. The findings were further examined in an independent cohort including 13 FTD patients, 79 patients with Alzheimer's disease and 18 healthy controls.

RESULTS

We found significantly altered protein levels in CSF from FTD patients compared to unaffected individuals (PMC and NC) for 26 proteins. The analysis show patterns of separation between unaffected individuals and FTD patients, especially for those with a clinical diagnosis of bvFTD. The most statistically significant differences in protein levels were found for VGF, TN-R, NPTXR, TMEM132D, PDYN and NF-M. Patients with FTD were found to have higher levels of TN-R and NF-M, and lower levels of VGF, NPTXR, TMEM132D and PDYN, compared to unaffected individuals. The main findings were reproduced in the independent cohort.

CONCLUSION

In this pilot study, we show a separation of FTD patients from unaffected individuals based on protein levels in CSF. Further investigation is required to explore the CSF profiles in larger cohorts, but the results presented here has the potential to enable future clinical utilization of these potential biomarkers within FTD.

Association between herpes simplex virus 1 exposure and the risk of depression in UK Biobank

BACKGROUND

Herpes simplex virus-1 (HSV-1) infection is reported to be associated with depression. But limited efforts were made to investigate the relationship between HSV-1 infection and the risk of depression, especially from the genetic perspective.

METHODS

In UK Biobank cohort, linear and logistic regression analyses were first performed to test the association of HSV-1 seropositivity/antibody with depression, including depression status (N = 2951) and Patient Health Questionnaire (PHQ) score (N = 2839). Using individual genotypic and phenotypic data from the UK Biobank, genome-wide environmental interaction study (GWEIS) was then conducted by PLINK2.0 to evaluate gene × HSV-1 interacting effect on the risk of depression. Finally, gene set enrichment analysis was conducted to identify the biological pathways involved in the observed gene × HSV-1 interaction for depression.

RESULT

In UK Biobank cohort, significant associations were observed between depression status and HSV-1 (odds ratio [OR] = 1.09; 95% confidence interval [CI], 1.02-1.16; P = 2.40 × 10^-2 for HSV-1 antibody and OR = 1.28; 95% CI, 1.12-1.47, P = 2.59 × 10^-3 for HSV-1 seropositivity). GWEIS revealed four significant gene × HSV-1 interaction signals for PHQ score (all P < 5.0 × 10^-8 ) and the leading loci was SULF2 (rs6094791, P = 8.60 × 10^-9 ). Pathway analyses identified 21 pathways for PHQ score and 19 for depression status, including multiple neural development- and immune-related ones, such as KEGG_NEUROACTIVE_LIGAND_RECEPTOR_INTERACTION (false discovery rate [FDR] = 3.18 × 10^-2 ) for depression and LU_AGING_BRAIN_UP (FDR = 4.21 × 10^-2 ) for PHQ score.

CONCLUSION

Our results suggested that HSV-1 was associated with the risk of depression, which was modulated by the several genes that were related to the nerve development or immune function.

A single-cell transcriptomic atlas of primate pancreatic islet aging

Aging-related degeneration of pancreatic islet cells contributes to impaired glucose tolerance and diabetes. Endocrine cells age heterogeneously, complicating the efforts to unravel the molecular drivers underlying endocrine aging. To overcome these obstacles, we undertook single-cell RNA sequencing of pancreatic islet cells obtained from young and aged non-diabetic cynomolgus monkeys. Despite sex differences and increased transcriptional variations, aged β-cells showed increased unfolded protein response (UPR) along with the accumulation of protein aggregates. We observed transcriptomic dysregulation of UPR components linked to canonical ATF6 and IRE1 signaling pathways, comprising adaptive UPR during pancreatic aging. Notably, we found aging-related β-cell-specific upregulation of HSP90B1, an endoplasmic reticulum-located chaperone, impeded high glucose-induced insulin secretion. Our work decodes aging-associated transcriptomic changes that underlie pancreatic islet functional decay at single-cell resolution and indicates that targeting UPR components may prevent loss of proteostasis, suggesting an avenue to delaying β-cell aging and preventing aging-related diabetes.

Transcriptome analysis of neural progenitor cells derived from Lowe syndrome induced pluripotent stem cells: identification of candidate genes for the neurodevelopmental and eye manifestations

BACKGROUND

Lowe syndrome (LS) is caused by loss-of-function mutations in the X-linked gene OCRL, which codes for an inositol polyphosphate 5-phosphatase that plays a key role in endosome recycling, clathrin-coated pit formation, and actin polymerization. It is characterized by congenital cataracts, intellectual and developmental disability, and renal proximal tubular dysfunction. Patients are also at high risk for developing glaucoma and seizures. We recently developed induced pluripotent stem cell (iPSC) lines from three patients with LS who have hypomorphic variants affecting the 3' end of the gene, and their neurotypical brothers to serve as controls.

METHODS

In this study, we used RNA sequencing (RNA-seq) to obtain transcriptome profiles in LS and control neural progenitor cells (NPCs).

RESULTS

In a comparison of the patient and control NPCs (n = 3), we found 16 differentially expressed genes (DEGs) at the multiple test adjusted p value (padj) < 0.1, with nine at padj < 0.05. Using nominal p value < 0.05, 319 DEGs were detected. The relatively small number of DEGs could be due to the fact that OCRL is not a transcription factor per se, although it could have secondary effects on gene expression through several different mechanisms. Although the number of DEGs passing multiple test correction was small, those that were found are quite consistent with some of the known molecular effects of OCRL protein, and the clinical manifestations of LS. Furthermore, using gene set enrichment analysis (GSEA), we found that genes increased expression in the patient NPCs showed enrichments of several gene ontology (GO) terms (false discovery rate < 0.25): telencephalon development, pallium development, NPC proliferation, and cortex development, which are consistent with a condition characterized by intellectual disabilities and psychiatric manifestations. In addition, a significant enrichment among the nominal DEGs for genes implicated in autism spectrum disorder (ASD) was found (e.g., AFF2, DNER, DPP6, DPP10, RELN, CACNA1C), as well as several that are strong candidate genes for the development of eye problems found in LS, including glaucoma. The most notable example is EFEMP1, a well-known candidate gene for glaucoma and other eye pathologies.

CONCLUSION

Overall, the RNA-seq findings present several candidate genes that could help explain the underlying basis for the neurodevelopmental and eye problems seen in boys with LS.

Further evidence of involvement of TMEM132E in autosomal recessive nonsyndromic hearing impairment

Autosomal-recessive (AR) nonsyndromic hearing impairment (NSHI) displays a high degree of genetic heterogeneity with >100 genes identified. Recently, TMEM132E, which is highly expressed in inner hair cells, was suggested as a novel ARNSHI gene for DFNB99. A missense variant c.1259G>A: p.(Arg420Gln) in TMEM132E was identified that segregated with ARNSHI in a single Chinese family with two affected members. In the present study, a family of Pakistani origin with prelingual profound sensorineural hearing impairment displaying AR mode of inheritance was investigated via exome and Sanger sequencing. Compound heterozygous variants c.382G>T: p.(Ala128Ser) and c.2204C>T: p.(Pro735Leu) in TMEM132E were observed in affected but not in unaffected family members. TMEM132E variants identified in this and the previously reported ARNSHI family are located in the extracellular domain. In conclusion, we present a second ARNSHI family with TMEM132E variants which strengthens the evidence of the involvement of this gene in the etiology of ARNSHI.

Roadmap of DNA methylation in breast cancer identifies novel prognostic biomarkers

BACKGROUND

Breast cancer is a highly heterogeneous disease resulting in diverse clinical behaviours and therapeutic responses. DNA methylation is a major epigenetic alteration that is commonly perturbed in cancers. The aim of this study is to characterize the relationship between DNA methylation and aberrant gene expression in breast cancer.

METHODS

We analysed DNA methylation and gene expression profiles from breast cancer tissue and matched normal tissue in The Cancer Genome Atlas (TCGA). Genome-wide differential methylation analysis and methylation-gene expression correlation was performed. Gene expression changes were subsequently validated in the METABRIC dataset. The Oncoscore tool was used to identify genes that had previously been associated with cancer in the literature. A subset of genes that had not previously been studied in cancer was chosen for further analysis.

RESULTS

We identified 368 CpGs that were differentially methylated between tumor and normal breast tissue (∆β > 0.4). Hypermethylated CpGs were overrepresented in tumor tissue and were found predominantly (56%) in upstream promoter regions. Conversely, hypomethylated CpG sites were found primarily in the gene body (66%). Expression analysis revealed that 209 of the differentially-methylated CpGs were located in 169 genes that were differently expressed between normal and breast tumor tissue. Methylation-expression correlations were predominantly negative (70%) for promoter CpG sites and positive (74%) for gene body CpG sites. Among these differentially-methylated and differentially-expressed genes, we identified 7 that had not previously been studied in any form of cancer. Three of these, TDRD10, PRAC2 and TMEM132C, contained CpG sites that showed diagnostic and prognostic value in breast cancer, particularly in estrogen-receptor (ER)-positive samples. A pan-cancer analysis confirmed differential expression of these genes together with diagnostic and prognostic value of their respective CpG sites in multiple cancer types.

CONCLUSION

We have identified 368 DNA methylation changes that characterize breast cancer tumor tissue, of which 209 are associated with genes that are differentially-expressed in the same samples. Novel DNA methylation markers were identified, of which cg12374721 (PRAC2), cg18081940 (TDRD10) and cg04475027 (TMEM132C) show promise as diagnostic and prognostic markers in breast cancer as well as other cancer types.