Genetic characterization of the Latvian local goat breed and genetic traits associated with somatic cell count

The Latvian local goat (LVK) breed represents the only native domestic goat breed in Latvia, but its limited population places it within the endangered category. However, the LVK breed has not yet undergone a comprehensive genetic characterization. Therefore, we completed whole genome sequencing to reveal the genetic foundation of the LVK breed while identifying genetic traits linked to the somatic cell count (SCC) levels. The study included 40 genomes of LVK goats sequenced to acquire at least 35x or 10x coverage. A Principal component analysis, a genetic distance tree, and an admixture analysis showed LVK's similarity to some European breeds, such as Finnish Landrace, Alpine, and Saanen, which aligns with the breed's history. An analysis of genome-wide heterozygosity, nucleotide diversity, and LD analysis indicated that the LVK population exhibits substantial levels of genetic diversity. LVK genome was dominated by short runs of homozygosity (ROHs, ≤ 500 kb) with a median length of 25 kb. With FROH 2.49%, average inbreeding levels were low; however, FROH ranged broadly from 0.13 to 12.2%. With the exception of one pure-blood breeding buck exhibiting FROH of 9.3% and FSNP of 8.5%, animals with at least 66% LVK ancestry showed moderate or no inbreeding. Overall, this study demonstrated that the LVK goats can be differentiated from imported breeds, although the population has a complex genetic structure. We were able to identify potential genetic traits associated with SCC levels, although the kinship of the animals and the heterogenic substructure of the population might have largely influenced the association analysis. We identified 26 genetic variants associated with SCC levels, which included the potentially relevant SNP rs662053371 in the OSBPL8 gene, indicating a potential signal linked to lipid metabolism in goats. To conclude, these findings present valuable insight into the genetic structure of the LVK breed for the conservation of local genetic resources.

Variations in the Relative Abundance of Gut Bacteria Correlate with Lipid Profiles in Healthy Adults

The gut microbiome is a versatile system regulating numerous aspects of host metabolism. Among other traits, variations in the composition of gut microbial communities are related to blood lipid patterns and hyperlipidaemia, yet inconsistent association patterns exist. This study aims to assess the relationships between the composition of the gut microbiome and variations in lipid profiles among healthy adults. This study used data and samples from 23 adult participants of a previously conducted dietary intervention study. Circulating lipid measurements and whole-metagenome sequences of the gut microbiome were derived from 180 blood and faecal samples collected from eight visits distributed across an 11-week study. Lipid-related variables explained approximately 4.5% of the variation in gut microbiome compositions, with higher effects observed for total cholesterol and high-density lipoproteins. Species from the genera Odoribacter, Anaerostipes, and Parabacteroides correlated with increased serum lipid levels, whereas probiotic species like Akkermansia muciniphila were more abundant among participants with healthier blood lipid profiles. An inverse correlation with serum cholesterol was also observed for Massilistercora timonensis, a player in regulating lipid turnover. The observed correlation patterns add to the growing evidence supporting the role of the gut microbiome as an essential regulator of host lipid metabolism.

Genomic Characterization and Initial Insight into Mastitis-Associated SNP Profiles of Local Latvian Bos taurus Breeds

Latvia has two local Bos taurus breeds-Latvian Brown (LBG) and Latvian Blue (LZG)-characterized by a good adaptation to the local climate, longevity, and high fat and protein contents in milk. Since these are desired traits in the dairy industry, this study investigated the genetic background of the LBG and LZG breeds and identified the genetic factors associated with mastitis. Blood and semen samples were acquired, and whole genome sequencing was then performed to acquire a genomic sequence with at least 35× or 10× coverage. The heterozygosity, nucleotide diversity, and LD analysis indicated that LBG and LZG cows have similar levels of genetic diversity compared to those of other breeds. An analysis of the population structure revealed that each breed clustered together, but the overall differentiation between the breeds was small. The highest genetic variance was observed in the LZG breed compared with the LBG breed. Our results show that SNP rs721295390 is associated with mastitis in the LBG breed, and SNPs rs383806754, chr29:43998719CG>C, and rs462030680 are associated with mastitis in the LZG breed. This study shows that local Latvian LBG and LZG breeds have a pronounced genetic differentiation, with each one suggesting its own mastitis-associated SNP profile.

Detection of SARS-CoV-2 RNA in wastewater and importance of population size assessment in smaller cities: An exploratory case study from two municipalities in Latvia

Wastewater-based epidemiology (WBE) has regained global importance during the COVID-19 pandemic. The mobility of people and other factors, such as precipitation and irregular inflow of industrial wastewater, are complicating the estimation of the disease prevalence through WBE, which is crucial for proper crisis management. These estimations are particularly challenging in urban areas with moderate or low numbers of inhabitants in situations where movement restrictions are not adopted (as in the case of Latvia) because residents of smaller municipalities tend to be more mobile and less strict in following the rules and measures of disease containment. Thus, population movement can influence the outcome of WBE measurements significantly and may not reflect the actual epidemiological situation in the respective area. Here, we demonstrate that by combining the data of detected SARS-CoV-2 RNA copy number, 5-hydroxyindoleacetic acid (5-HIAA) analyses in wastewater and mobile call detail records it was possible to provide an accurate assessment of the COVID-19 epidemiological situation in towns that are small (COVID-19 28-day cumulative incidence r = 0.609 and 35-day cumulative incidence r = 0.89, p < 0.05) and medium-sized towns (COVID-19 21-day cumulative incidence r = 0.997, 28-day cumulative incidence r = 0.98 and 35-day cumulative incidence r = 0.997, p < 0.05). This is the first study demonstrating WBE for monitoring COVID-19 outbreaks in Latvia. We demonstrate that the application of population size estimation measurements such as total 5-HIAA and call detail record data improve the accuracy of the WBE approach.

Baseline gut microbiome composition predicts metformin therapy short-term efficacy in newly diagnosed type 2 diabetes patients

Background

The study was conducted to investigate the effects of metformin treatment on the human gut microbiome’s taxonomic and functional profile in the Latvian population, and to evaluate the correlation of these changes with therapeutic efficacy and tolerance.

Methods

In this longitudinal observational study, stool samples for shotgun metagenomic sequencing-based analysis were collected in two cohorts. The first cohort included 35 healthy nondiabetic individuals (metformin dose 2x850mg/day) at three time-points during metformin administration. The second cohort was composed of 50 newly-diagnosed type 2 diabetes patients (metformin dose–determined by an endocrinologist) at two concordant times. Patients were defined as Responders if their HbA1c levels during three months of metformin therapy had decreased by ≥12.6 mmol/mol (1%), while in Non-responders HbA1c were decreased by <12.6 mmol/mol (1%).

Results

Metformin reduced the alpha diversity of microbiota in healthy controls (p = 0.02) but not in T2D patients. At the species level, reduction in the abundance of Clostridium bartlettii and Barnesiella intestinihominis, as well as an increase in the abundance of Parabacteroides distasonis and Oscillibacter unclassified overlapped between both study groups. A large number of group-specific changes in taxonomic and functional profiles was observed. We identified an increased abundance of Prevotella copri (FDR = 0.01) in the Non-Responders subgroup, and enrichment of Enterococcus faecium, Lactococcus lactis, Odoribacter, and Dialister at baseline in the Responders group. Various taxonomic units were associated with the observed incidence of side effects in both cohorts.

Conclusions

Metformin effects are different in T2D patients and healthy individuals. Therapy induced changes in the composition of gut microbiome revealed possible mediators of observed short-term therapeutic effects. The baseline composition of the gut microbiome may influence metformin therapy efficacy and tolerance in T2D patients and could be used as a powerful prediction tool.

Whole-blood transcriptome profiling reveals signatures of metformin and its therapeutic response

Metformin, a biguanide agent, is the first-line treatment for type 2 diabetes mellitus due to its glucose-lowering effect. Despite its wide application in the treatment of multiple health conditions, the glycemic response to metformin is highly variable, emphasizing the need for reliable biomarkers. We chose the RNA-Seq-based comparative transcriptomics approach to evaluate the systemic effect of metformin and highlight potential predictive biomarkers of metformin response in drug-naïve volunteers with type 2 diabetes in vivo. The longitudinal blood-derived transcriptome analysis revealed metformin-induced differential expression of novel and previously described genes involved in cholesterol homeostasis (SLC46A1 and LRP1), cancer development (CYP1B1, STAB1, CCR2, TMEM176B), and immune responses (CD14, CD163) after administration of metformin for three months. We demonstrate for the first time a transcriptome-based molecular discrimination between metformin responders (delta HbA1c ≥ 1% or 12.6 mmol/mol) and non-responders (delta HbA1c < 1% or 12.6 mmol/mol), that is determined by expression levels of 56 genes, explaining 13.9% of the variance in the therapeutic efficacy of the drug. Moreover, we found a significant upregulation of IRS2 gene (log₂FC 0.89) in responders compared to non-responders before the use of metformin. Finally, we provide evidence for the mitochondrial respiratory complex I as one of the factors related to the high variability of the therapeutic response to metformin in patients with type 2 diabetes mellitus.

Significantly altered peripheral blood cell DNA methylation profile as a result of immediate effect of metformin use in healthy individuals

BACKGROUND

Metformin is a widely prescribed antihyperglycemic agent that has been also associated with multiple therapeutic effects in various diseases, including several types of malignancies. There is growing evidence regarding the contribution of the epigenetic mechanisms in reaching metformin's therapeutic goals; however, the effect of metformin on human cells in vivo is not comprehensively studied. The aim of our study was to examine metformin-induced alterations of DNA methylation profiles in white blood cells of healthy volunteers, employing a longitudinal study design.

RESULTS

Twelve healthy metformin-naïve individuals where enrolled in the study. Genome-wide DNA methylation pattern was estimated at baseline, 10 h and 7 days after the start of metformin administration. The whole-genome DNA methylation analysis in total revealed 125 differentially methylated CpGs, of which 11 CpGs and their associated genes with the most consistent changes in the DNA methylation profile were selected: POFUT2, CAMKK1, EML3, KIAA1614, UPF1, MUC4, LOC727982, SIX3, ADAM8, SNORD12B, VPS8, and several differentially methylated regions as novel potential epigenetic targets of metformin. The main functions of the majority of top-ranked differentially methylated loci and their representative cell signaling pathways were linked to the well-known metformin therapy targets: regulatory processes of energy homeostasis, inflammatory responses, tumorigenesis, and neurodegenerative diseases.

CONCLUSIONS

Here we demonstrate for the first time the immediate effect of short-term metformin administration at therapeutic doses on epigenetic regulation in human white blood cells. These findings suggest the DNA methylation process as one of the mechanisms involved in the action of metformin, thereby revealing novel targets and directions of the molecular mechanisms underlying the various beneficial effects of metformin.

TRIAL REGISTRATION

EU Clinical Trials Register, 2016-001092-74. Registered 23 March 2017, https://www.clinicaltrialsregister.eu/ctr-search/trial/2016-001092-74/LV .

Association of metformin administration with gut microbiome dysbiosis in healthy volunteers

BACKGROUND

Metformin is a widely used first-line drug for treatment of type 2 diabetes. Despite its advantages, metformin has variable therapeutic effects, contraindications, and side effects. Here, for the very first time, we investigate the short-term effect of metformin on the composition of healthy human gut microbiota.

METHODS

We used an exploratory longitudinal study design in which the first sample from an individual was the control for further samples. Eighteen healthy individuals were treated with metformin (2 × 850 mg) for 7 days. Stool samples were collected at three time points: prior to administration, 24 hours and 7 days after metformin administration. Taxonomic composition of the gut microbiome was analyzed by massive parallel sequencing of 16S rRNA gene (V3 region).

RESULTS

There was a significant reduction of inner diversity of gut microbiota observed already 24 hours after metformin administration. We observed an association between the severity of gastrointestinal side effects and the increase in relative abundance of common gut opportunistic pathogen Escherichia-Shigella spp. One week long treatment with metformin was associated with a significant decrease in the families Peptostreptococcaceae and Clostridiaceae_1 and four genera within these families.

CONCLUSIONS

Our results are in line with previous findings on the capability of metformin to influence gut microbiota. However, for the first time we provide evidence that metformin has an immediate effect on the gut microbiome in humans. It is likely that this effect results from the increase in abundance of opportunistic pathogens and further triggers the occurrence of side effects associated with the observed dysbiosis. An additional randomized controlled trial would be required in order to reach definitive conclusions, as this is an exploratory study without a placebo control arm. Our findings may be further used to create approaches that improve the tolerability of metformin.

Single nucleotide polymorphisms in the intergenic region between metformin transporter OCT2 and OCT3 coding genes are associated with short-term response to metformin monotherapy in type 2 diabetes mellitus patients

OBJECTIVES

High variability in clinical response to metformin is often observed in type 2 diabetes (T2D) patients, and it highlights the need for identification of genetic components affecting the efficiency of metformin therapy. Aim of this observational study is to evaluate the role of tagSNPs (tagging single nucleotide polymorphisms) from genomic regions coding for six metformin transporter genes with respect to the short-term efficiency.

DESIGN

102 tagSNPs in 6 genes coding for metformin transporters were genotyped in the group of 102 T2D patients treated with metformin for 3 months.

METHODS

Most significant hits were analyzed in the group of 131 T2D patients from Slovakia. Pharmacokinetic study in 25 healthy nondiabetic volunteers was conducted to investigate the effects of identified polymorphisms.

RESULTS

In the discovery group of 102 patients, minor alleles of rs3119309, rs7757336 and rs2481030 were significantly nominally associated with metformin inefficiency (P = 1.9 × 10^-6 to 8.1 × 10^-6). Effects of rs2481030 and rs7757336 did not replicate in the group of 131 T2DM patients from Slovakia alone, whereas rs7757336 was significantly associated with a reduced metformin response in combined group. In pharmacokinetic study, group of individuals harboring risk alleles of rs7757336 and rs2481030 displayed significantly reduced AUC_∞ of metformin in plasma.

CONCLUSIONS

For the first time, we have identified an association between the lack of metformin response and SNPs rs3119309 and rs7757336 located in the 5' flanking region of the genes coding for Organic cation transporter 2 and rs2481030 located in the 5' flanking region of Organic cation transporter 3 that was supported by the results of a pharmacokinetic study on 25 healthy volunteers.

Genome-wide analysis reveals DNA methylation markers that vary with both age and obesity

The combination of the obesity epidemic and an aging population presents growing challenges for the healthcare system. Obesity and aging are major risk factors for a diverse number of diseases and it is of importance to understand their interaction and the underlying molecular mechanisms. Herein the authors examined the methylation levels of 27578 CpG sites in 46 samples from adult peripheral blood. The effect of obesity and aging was ascertained with general linear models. More than one hundred probes were correlated to aging, nine of which belonged to the KEGG group map04080. Additionally, 10 CpG sites had diverse methylation profiles in obese and lean individuals, one of which was the telomerase catalytic subunit (TERT). In eight of ten cases the methylation change was reverted between obese and lean individuals. One region proved to be differentially methylated with obesity (LINC00304) independent of age. This study provides evidence that obesity influences age driven epigenetic changes, which provides a molecular link between aging and obesity. This link and the identified markers may prove to be valuable biomarkers for the understanding of the molecular basis of aging, obesity and associated diseases.

The role of common and rare MC4R variants and FTO polymorphisms in extreme form of obesity

Melanocortin 4 receptor (MC4R) is an important regulator of food intake and number of studies report genetic variations influencing the risk of obesity. Here we explored the role of common genetic variation from MC4R locus comparing with SNPs from gene FTO locus, as well as the frequency and functionality of rare MC4R mutations in cohort of 380 severely obese individuals (BMI > 39 kg/m(2)) and 380 lean subjects from the Genome Database of Latvian Population (LGDB). We found correlation for two SNPs--rs11642015 and rs62048402 in the fat mass and obesity-associated protein (FTO) with obesity but no association was detected for rs17782313 located in the MC4R locus in these severely obese individuals. We sequenced the whole gene MC4R coding region in all study subjects and found five previously known heterozygous non-synonymous substitutions V103I, I121T, S127L, V166I and I251L. Expression in mammalian cells showed that the S127L, V166I and double V103I/S127L mutant receptors had significantly decreased quantity at the cell surface compared to the wild type MC4R. We carried out detailed functional analysis of V166I that demonstrated that, despite low abundance in plasma membrane, the V166I variant has lower EC50 value upon αMSH activation than the wild type receptor, while the level of AGRP inhibition was decreased, implying that V166I cause hyperactive satiety signalling. Overall, this study suggest that S127L may be the most frequent functional MC4R mutation leading to the severe obesity in general population and provides new insight into the functionality of population based variants of the MC4R.

Determination of the obesity-associated gene variants within the entire FTO gene by ultra-deep targeted sequencing in obese and lean children

BACKGROUND

The Fat mass and obesity-associated gene (FTO) was the first gene reliably associated with body mass index in genome-wide association studies on a population level. At present, the genetic variations within the FTO gene are still the common variants that have the largest influence on body mass index.

METHODS

In the current study, we amplified the entire FTO gene, in total 412 Kbp, in over 200 long-range PCR fragments from each individual, from 524 severely obese and 527 lean Swedish children, and sequenced the products as two DNA pools using massive parallel sequencing (SOLiD).

RESULTS

The sequencing achieved very high coverage (median 18 000 reads) and we detected and estimated allele frequencies for 705 single nucleotide polymorphisms (SNPs) (19 novel) and 40 indels (24 novel) using a sophisticated statistical approach to remove false-positive SNPs. We identified 19 obesity-associated SNPs within intron one of the FTO gene, and validated our findings with genotyping. Ten of the validated obesity-associated SNPs have a stronger obesity association (P<0.007) than the commonly studied rs9939609 SNP (P<0.012).

CONCLUSIONS

This study provides a comprehensive obesity-associated variation map of FTO, identifies novel lead SNPs and evaluates putative causative variants. We conclude that intron one is the only region within the FTO gene associated with obesity, and finally, we establish next generation sequencing of pooled DNA as a powerful method to investigate genetic association with complex diseases and traits.

Stronger association of common variants in TCF7L2 gene with nonobese type 2 diabetes in the Latvian population

Polymorphisms in the gene coding for transcription factor 7 like 2 (TCF7L2) are recognized as the strongest common genetic risk factors for type 2 diabetes (T2D) across multiple ethnicities. This study was conducted to evaluate an association between TCF7L2 variants and diabetes susceptibility in the population of Latvia. We genotyped 4 single nucleotide polymorphisms (SNP) rs7901695, rs7903146, rs11196205 and rs12255372 in 1 093 controls and 1 043 diabetic subjects. Association with T2D was found for 3 SNPs rs7901695, rs7903146 and rs12255372 in the whole sample (under an additive genetic model, the adjusted odds ratios (OR) were 1.26, 95% CI [1.08-1.48], P=0.003; OR=1.32, 95% CI [1.12-1.55], P=0.001 and OR=1.35, 95% CI [1.15-1.60], P=0.0004 respectively). In addition observed effects on T2D susceptibility for analysed SNPs were higher among subjects with BMI under 30 kg/m². The impact of TCF7L2 variation on T2D risk in Latvian population is compatible with that demonstrated by a range of studies conducted in various ethnic groups.

Identification of domains responsible for specific membrane transport and ligand specificity of the ACTH receptor (MC2R)

The adrenocorticotropic hormone (ACTH) receptor has highly specific membrane expression that is limited to adrenal cells; in other cell types the polypeptide fails to be transported to the cell surface. Unlike other evolutionarily related members of the melanocortin receptor family (MC1R-MC5R) that recognize different melanocortin peptides, ACTHR (MC2R) binds only ACTH. We used a mutagenesis approach involving systematic construction of chimeric ACTHR/MC4R receptors to identify the domains determining the selectivity of ACTHR membrane transport and ACTH binding. In total 15 chimeric receptors were created by replacement of selected domains of human ACTHR with the corresponding regions of human MC4R. We developed an analytical method to accurately quantify cell-membrane localization of recombinant receptors fused with enhanced green fluorescent protein by confocal fluorescence microscopy. The chimeric receptors were also tested for their ability to bind ACTH (1-24) and the melanocyte-stimulating hormone (MSH) analog, Nle4, DPhe7-alpha-MSH, and to induce a cAMP response. Our results indicate that substitution of the MC4R N-terminal segment with the homologous segment of ACTHR significantly decreased membrane transport. We also identified another signal localized in the third and fourth transmembrane regions as the main determinant of ACTHR intracellular retention. In addition, we found that the fourth and fifth transmembrane domains of the ACTHR are involved in ACTH binding selectivity. We discuss the mechanisms involved in bypassing these arrest signals via an interaction with melanocortin 2 receptor accessory protein (MRAP) and the possible mechanisms that determine the high ligand-binding specificity of ACTHR.

Association between a rare SNP in the second intron of human Agouti related protein gene and increased BMI

BACKGROUND

The agouti related protein (AGRP) is an endogenous antagonist of the melanocortin 4 receptor and is one of the most potent orexigenic factors. The aim of the present study was to assess the genetic variability of AGRP gene and investigate whether the previously reported SNP rs5030980 and the rs11575892, a SNP that so far has not been studied with respect to obesity is associated with increased body mass index (BMI).

METHODS

We determined the complete sequence of the AGRP gene and upstream promoter region in 95 patients with severe obesity (BMI > 35 kg/m2). Three polymorphisms were identified: silent mutation c.123G>A (rs34123523) in the second exon, non-synonymous mutation c.199G>A (rs5030980) and c.131-42C>T (rs11575892) located in the second intron. We further screened rs11575892 in a selected group of 1135 and rs5030980 in group of 789 participants from the Genome Database of Latvian Population and Latvian State Research Program Database.

RESULTS

The CT heterozygotes of rs11575892 had significantly higher mean BMI value (p = 0.027). After adjustment for age, gender and other significant non-genetic factors (presence of diseases), the BMI levels remained significantly higher in carriers of the rs11575892 T allele (p = 0.001). The adjusted mean BMI value of CC genotype was 27.92 +/- 1.01 kg/m2 (mean, SE) as compared to 30.97 +/- 1.03 kg/m2 for the CT genotype. No association was found between rs5030980 and BMI.

CONCLUSION

This study presents an association of rare allele of AGRP polymorphism in heterozygous state with increased BMI. The possible functional effects of this polymorphism are unclear but may relate to splicing defects.

The evolutionary history and tissue mapping of GPR123: specific CNS expression pattern predominantly in thalamic nuclei and regions containing large pyramidal cells

The Adhesion family of G protein-coupled receptors (GPCRs) includes 33 receptors and is the second largest GPCR family. Most of these proteins are still orphans and fairly little is known of their tissue distribution and evolutionary context. We report the evolutionary history of the Adhesion family protein GPR123 as well as mapping of GPR123 mRNA expression in mouse and rat using in situ hybridization and real-time PCR, respectively. GPR123 was found to be well conserved within the vertebrate lineage, especially within the transmembrane regions and in the distal part of the cytoplasmic tail, containing a potential PDZ binding domain. The real-time PCR data indicates that GPR123 is predominantly expressed in CNS. The in situ data show high expression in thalamic nuclei and regions containing large pyramidal cells like cortex layers 5 and 6 and subiculum. Moreover, we found distinct expression in amygdala, hypothalamus, inferior olive and spinal cord. The CNS specific expression, together with the high sequence conservation between the vertebrate sequences investigated, indicate that GPR123 may have an important role in the regulation of neuronal signal transduction.