ALBERT

Description: There are still several unanswered questions regarding about ancient events in the Lithuanian population. The Lithuanians, as the subject of this study, are of great interest as they represent a partially isolated population maintaining an ancient genetic composition and show genetic uniqueness in European comparisons. To elucidate the genetic relationships between the Lithuania and North-Eastern European and West Siberian populations, we analyzed the population structure, effective population size, and recent positive selection from genome-wide single nucleotide polymorphism (SNP) data. We identified the close genetic proximity of Lithuanians to neighboring populations (Latvians, Estonians, Belarusians) and in part with West and South Slavs (Poles, Slovaks, and Slovenians), however, with particular genetic distinctiveness. The estimated long-term Ne values ranged from ~5900 in the Estonian population to ~2400 in the South Russian population. The divergence times between the Lithuanian and study populations ranged from 240 to 12,871 YBP. We also found evidence of selection in 24 regions, 21 of which have not been discovered in previous analyses of selection. Undoubtedly, the genetic diversity analysis of geographically specific regions may provide new insights into microevolutionary processes affecting local human populations.

Description: Wolfram syndrome (WS) 1 is a rare monogenic neurodegenerative disorder caused by mutations in the gene encoding WFS1. Knowledge of the pathophysiology of WS is incomplete and to date, there is no treatment available. Here, we describe early deviations in the renin-angiotensin-aldosterone system (RAAS) and bradykinin pathway (kallikrein kinin system, KKS) observed in a rat model of WS (Wfs1 KO) and the modulative effect of glucagon-like peptide-1 receptor agonist liraglutide (LIR) and anti-epileptic drug valproate (VPA), which have been proven effective in delaying WS progression in WS animal models. We found that the expression of key receptors of the RAAS and KKS, Agtr2 and Bdkrb1, were drastically downregulated both in vitro and in vivo at an early stage in a rat model of WS. Moreover, in Wfs1, KO serum aldosterone levels were substantially decreased and bradykinin levels increased compared to WT animals. Neither treatment nor their combination affected the gene expression levels seen in the Wfs1 KO animals. However, all the treatments elevated serum aldosterone and decreased bradykinin in the Wfs1 KO rats, as well as increasing angiotensin II levels independent of genotype. Altogether, our results indicate that Wfs1 deficiency might disturb the normal functioning of RAAS and KKS and that LIR and VPA have the ability to modulate these systems.

Description: The identification of human remains is challenging mostly due to the bad condition of the remains and the available background information that is sometimes limited. The current case report is related to the identification of an unknown soldier from the Estonian War of Independence (1918–1920). The case includes an anthropological study of the remains, examinations of documents found with the exhumed remains, and kinship estimations based on archival documents, and DNA analyses. As the preliminary data pointed to remains of male origin, Y-chromosomal STR (short tandem repeat) analyses of 22 Y-STR loci were used to analyze the exhumed teeth. Reference samples from individuals from two paternal lineages were collected based on archival documents. Y-chromosomal STR results for the tooth samples were consistent with a patrilineal relationship to only one reference sample out of two proposed paternal lineages. Based on the provided pedigrees in the consistent case, the Y-STR results are approximately four million times more likely if the tooth sample originated from an individual related along the paternal line to the matching reference sample, than if the tooth sample originated from another person in the general population. Special considerations have to be met when limited evidence is available.

Description: The detection of circulating microRNA (miRNA)-based biomarkers represents an innovative, non-invasive method for the early detection of cancer. However, the low concentration of miRNAs released in body fluids and the difficult identification of the tumor site have limited their clinical use as effective cancer biomarkers. To evaluate if ultrasound treatment could amplify the release of extracellular cancer biomarkers, we treated a panel of prostate cancer (PCa) cell lines with an ultrasound-based prototype and profiled the release of miRNAs in the extracellular space, with the aim of identifying novel miRNA-based biomarkers that could be used for PCa diagnosis and the monitoring of tumor evolution. We provide evidence that US-mediated sonoporation amplifies the release of miRNAs from both androgen-dependent (AD) and -independent (AI) PCa cells. We identified four PCa-related miRNAs, whose levels in LNCaP and DU145 supernatants were significantly increased following ultrasound treatment: mir-629-5p, mir-374-5p, mir-194-5p, and let-7d-5p. We further analyzed a publicly available dataset of PCa, showing that the serum expression of these novel miRNAs was upregulated in PCa patients compared to controls, thus confirming their clinical relevance. Our findings highlight the potential of using ultrasound to identify novel cell-free miRNAs released from cancer cells, with the aim of developing new biomarkers with diagnostic and predictive value.

Description: Trisomy 21 (Ts21) causes alterations in skeletal development resulting in decreased bone mass, shortened stature and weaker bones in individuals with Down syndrome (DS). There is a sexual dimorphism in bone mineral density (BMD) deficits associated with DS with males displaying earlier deficits than females. The relationships between causative trisomic genes, cellular mechanisms, and influence of sex in DS skeletal abnormalities remain unknown. One hypothesis is that the low bone turnover phenotype observed in DS results from attenuated osteoblast function, contributing to impaired trabecular architecture, altered cortical geometry, and decreased mineralization. DYRK1A, found in three copies in humans with DS, Ts65Dn, and Dp1Tyb DS model mice, has been implicated in the development of postnatal skeletal phenotypes associated with DS. Reduced copy number of Dyrk1a to euploid levels from conception in an otherwise trisomic Ts65Dn mice resulted in a rescue of appendicular bone deficits, suggesting DYRK1A contributes to skeletal development and homeostasis. We hypothesized that reduction of Dyrk1a copy number in trisomic osteoblasts would improve cellular function and resultant skeletal structural anomalies in trisomic mice. Female mice with a floxed Dyrk1a gene (Ts65Dn,Dyrk1afl/wt) were mated with male Osx-Cre+ (expressed in osteoblasts beginning around E13.5) mice, resulting in reduced Dyrk1a copy number in mature osteoblasts in Ts65Dn,Dyrk1a+/+/Osx-Cre P42 male and female trisomic and euploid mice, compared with littermate controls. Male and female Ts65Dn,Dyrk1a+/+/+ (3 copies of DYRK1A in osteoblasts) and Ts65Dn,Dyrk1a+/+/Osx-Cre (2 copies of Dyrk1a in osteoblasts) displayed similar defects in both trabecular architecture and cortical geometry, with no improvements with reduced Dyrk1a in osteoblasts. This suggests that trisomic DYRK1A does not affect osteoblast function in a cell-autonomous manner at or before P42. Although male Dp1Tyb and Ts65Dn mice exhibit similar skeletal deficits at P42 in both trabecular and cortical bone compartments between euploid and trisomic mice, female Ts65Dn mice exhibit significant cortical and trabecular deficits at P42, in contrast to an absence of genotype effect in female Dp1Tyb mice in trabecular bone. Taken together, these data suggest skeletal deficits in DS mouse models and are sex and age dependent, and influenced by strain effects, but are not solely caused by the overexpression of Dyrk1a in osteoblasts. Identifying molecular and cellular mechanisms, disrupted by gene dosage imbalance, that are involved in the development of skeletal phenotypes associated with DS could help to design therapies to rescue skeletal deficiencies seen in DS.

Description: The developmentally programmed loss of a plant organ is called abscission. This process is characterized by the ultimate separation of adjacent cells in the abscission zone (AZ). The discovery of an American oil palm (Elaeis oleifera) variant that does not shed its has allowed for the study of the mechanisms of ripe fruit abscission in this species. A comparative transcriptome analysis was performed to compare the fruit AZs of the non-shedding E. oleifera variant to an individual of the same progeny that sheds its ripe fruit normally. The study provides evidence for widespread perturbation to gene expression in the AZ of the non-shedding variant, compared to the normal fruit-shedding control, and offers insight into abscission-related functions. Beyond the genes with known or suspected roles during organ abscission or indehiscence that were identified, a list of genes with hormone-related functions, including ethylene, jasmonic acid, abscisic acid, cytokinin and salicylic acid, in addition to reactive oxygen species (ROS) metabolism, transcriptional responses and signaling pathways, was compiled. The results also allowed a comparison between the ripe fruit abscission processes of the African and American oil palm species at the molecular level and revealed commonalities with environmental stress pathways.

Description: Milk fat and protein content are among key elements of milk quality, and they are attracting more attention in response to consumers′ demand for high-quality dairy products. To investigate the potential regulatory roles of DNA methylation underlying milk component yield, whole genome bisulfite sequencing was employed to profile the global DNA methylation patterns of mammary gland tissues from 17 Canada Holstein cows with various milk fat and protein contents. A total of 706, 2420 and 1645 differentially methylated CpG sites (DMCs) were found between high vs. low milk fat (HMF vs. LMF), high vs. low milk protein (HMP vs. LMP), and high vs. low milk fat and protein (HMFP vs. LMFP) groups, respectively (q value 〈 0.1). Twenty-seven, 56 and 67 genes harboring DMCs in gene regions (denoted DMC genes) were identified for HMF vs. LMF, HMP vs. LMP and HMFP vs. LMFP, respectively. DMC genes from HMP vs. LMP and HMFP vs. LMFP comparisons were significantly overrepresented in GO terms related to aerobic electron transport chain and/or mitochondrial ATP (adenosine triphosphate) synthesis coupled electron transport. A total of 83 (HMF vs. LMF), 708 (HMP vs. LMP) and 408 (HMFP vs. LMFP) DMCs were co-located with 87, 147 and 158 quantitative trait loci (QTL) for milk component and yield traits, respectively. In conclusion, the identified methylation changes are potentially involved in the regulation of milk fat and protein yields, as well as the variation in reported co-located QTLs.

Description: Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) is not only involved in carbohydrate metabolism, but also plays an important role in stress resistance. However, it has not been reported in Brassica oleracea. In this study, we performed a genome-wide identification of BoGAPDH in B. oleracea and performed cloning and expression analysis of one of the differentially expressed genes, BoGAPC. A total of 16 members of the BoGAPDH family were identified in B. oleracea, which were conserved, distributed unevenly on chromosomes and had tandem repeat genes. Most of the genes were down-regulated during self-pollination, and the highest expression was found in stigmas and sepals. Different transcriptome data showed that BoGAPDH genes were differentially expressed under stress, which was consistent with the results of qRT-PCR. We cloned and analyzed the differentially expressed gene BoGAPC and found that it was in the down-regulated mode 1 h after self-pollination, and the expression was the highest in the stigma, which was consistent with the result of GUS staining. The promoter region of the gene not only has stress response elements and plant hormone response elements, but also has a variety of specific elements for regulating floral organ development. Subcellular localization indicates that the BoGAPC protein is located in the cytoplasm and belongs to the active protein in the cytoplasm. The results of prokaryotic expression showed that the size of the BoGAPC protein was about 37 kDa, which was consistent with the expected results, indicating that the protein was induced in prokaryotic cells. The results of yeast two-hybrid and GST pull-down showed that the SRK kinase domain interacted with the BoGAPC protein. The above results suggest that the BoGAPDH family of B. oleracea plays an important role in the process of plant stress resistance, and the BoGAPC gene may be involved in the process of self-incompatibility in B. oleracea, which may respond to SI by encoding proteins directly interacting with SRK.

Description: Tumor mutational burden (TMB) is considered a potential biomarker for predicting the response and effect of immune checkpoint inhibitors (ICIs). To find specific gene mutations related to TMB and the prognosis of patients, the frequently mutated genes in gastric cancer patients from TCGA and ICGC were obtained and the correlation between gene mutation, TMB, and prognosis was analyzed. Furthermore, to clarify whether specific gene mutations can be used as predictive biomarkers of ICIs, a gene set enrichment analysis (GSEA) for immune pathways and an immune infiltration analysis were conducted. The results showed that CUB and Sushi multiple domains 1 (CSMD1) mutation (CSMD1-mut) were associated with higher TMB and better prognosis in patients. The genetic map showed that, compared with wild-type samples, the loss of chromosomes 4q, 5q, 8p, and 9p decreased and the status of microsatellite instability increased in the CSMD1-mut samples. The GSEA analysis showed that immune-related pathways were enriched in the CSMD1-mut samples. The immune infiltration analysis showed that the anti-tumor immune cells were upregulated and that the tumor-promoting immune cells were downregulated in the CSMD1-mut samples. The gene co-expression analysis showed that PD-L1 expression was higher in the CSMD1-mut samples. In summary, CSMD1-mut in gastric cancer was associated with increased TMB and favorable survival and may have potential significance in predicting the efficacy of anti-PD-L1.

Description: The additive genetic model as implemented in logistic regression has been widely used in genome-wide association studies (GWASs) for binary outcomes. Unfortunately, for many complex diseases, the underlying genetic models are generally unknown and a mis-specification of the genetic model can result in a substantial loss of power. To address this issue, the MAX3 test (the maximum of three separate test statistics) has been proposed as a robust test that performs plausibly regardless of the underlying genetic model. However, the original implementation of MAX3 utilizes the trend test so it cannot adjust for any covariates such as age and gender. This drawback has significantly limited the application of the MAX3 in GWASs, as covariates account for a considerable amount of variability in these disorders. In this paper, we extended the MAX3 and proposed the CMAX3 (covariate-adjusted MAX3) based on logistic regression. The proposed test yielded a similar robust efficiency as the original MAX3 while easily adjusting for any covariate based on the likelihood framework. The asymptotic formula to calculate the p-value of the proposed test was also developed in this paper. The simulation results showed that the proposed test performed desirably under both the null and alternative hypotheses. For the purpose of illustration, we applied the proposed test to re-analyze a case-control GWAS dataset from the Collaborative Studies on Genetics of Alcoholism (COGA). The R code to implement the proposed test is also introduced in this paper and is available for free download.