Genome-wide Single Nucleotide Polymorphism Research Articles

Intraspecific divergence in a coastal plant, Euphorbia jolkinii, at a major biogeographic boundary in East Asia.

Quaternary climatic fluctuations and long-distance seed dispersal across the sea are critical factors affecting the distribution of coastal plants, but the spatiotemporal nature of population expansion and distribution change of East Asian coastal plants during this period are rarely examined. To explore this process, we investigated the genome-wide phylogenetic patterns of Euphorbia jolkinii Boiss. (Euphorbiaceae), which grows widely on littoral areas of Japan, Korea, and Taiwan. We used plastome sequences and genome-wide single nucleotide polymorphisms in samples across the species range to reveal phylogeographic patterns and spatiotemporal distributional changes. We conducted ecological niche modeling for the present and the last glacial maximum (LGM). Genetic differentiation was observed between the northern and southern populations of E. jolkinii, separated by the major biogeographic boundary, the Tokara Gap. These two groups of populations differentiated during the glacial period and subsequently intermingled in the intermorainic areas of the central Ryukyu Islands after the LGM. Ecological niche models suggested that the potential range of E. jolkinii was restricted to southern Kyushu; however, it was widespread in the southern Ryukyu Islands and Taiwan during the LGM. This study provides evidence of genetic differentiation among coastal plant populations separated by the prominent biogeographical boundary. Although coastal plants are typically expected to maintain population connectivity through sea-drifted seed dispersal, our findings suggest that genetic differences may arise because of a combination of limited gene flow and changes in climate during the glacial period.

Novel insights into rosacea's role in cancer risk: A Mendelian randomization approach.

Chronic inflammation has been shown to promote cancer progression. Rosacea is indeed a long-term inflammatory skin condition and had been reported to link with increased risk for several types of malignancies, but evidence for causality is lacking. To systematically estimate the causal relationship between rosacea and several types of cancer, including cutaneous malignant melanoma (CMM), cutaneous squamous cell carcinoma (cSCC), basal cell carcinoma (BCC), actinic keratosis (AK), thyroid cancer, breast cancer, glioma and hepatic cancer, as well as explore the potential underlying pathogenesis. We conducted a bidirectional two-sample Mendelian randomization study to probe the potential causal relationships between rosacea and several types of cancer. Instrumental variables were established using genome-wide significant single nucleotide polymorphisms associated with rosacea and cancers. The assessment of causality was carried out through multiple methods, and the robustness of the results was evaluated via sensitivity analyses. There was no significant indication of causal effects of rosacea on CMM (pivw=0.71), cSCC (pivw=0.45), BCC (pivw=0.90), AK (pivw=0.73), thyroid cancer (pivw=0.59), glioma (pivw=0.15), and hepatic cancer (pivw=0.07), but the genetic risk of rosacea was associated with an increased susceptibility to human epidermal growth factor receptor (HER)-negative malignant neoplasm of breast (odds ratio [OR], 1.10; 95% confidence interval [CI], 1.02-1.18; pivw=0.01). TANK (TRAF family member associated nuclear factor kappa B (NFKB) activator) was identified as a common protective gene for both rosacea (OR, 0.90; 95% CI, 0.82-0.99; pivw=0.048) and HER-negative malignant neoplasm of the breast (OR, 0.86; 95% CI, 0.75-0.98; pivw=0.032), which was primarily enriched in the negative regulation of NF-κB signal transduction and may contribute to the genetic links between rosacea and this subtype of breast cancer. Our findings provide suggestive evidence for causal links between rosacea and HER-negative malignant neoplasm of the breast risk.

Impact of gene-by-trauma interaction in MDD-related multimorbidity clusters

BackgroundMajor depressive disorder (MDD) is considerably heterogeneous in terms of comorbidities, which may hamper the disentanglement of its biological mechanism. In a previous study, we classified the lifetime trajectories of MDD-related multimorbidities into seven distinct clusters, each characterized by unique genetic and environmental risk-factor profiles. The current objective was to investigate genome-wide gene-by-environment (G × E) interactions with childhood trauma burden, within the context of these clusters. MethodsWe analyzed 77,519 participants and 6,266,189 single-nucleotide polymorphisms (SNPs) of the UK Biobank database. Childhood trauma burden was assessed using the Childhood Trauma Screener (CTS). For each cluster, Plink 2.0 was used to calculate SNP × CTS interaction effects on the participants' cluster membership probabilities. We especially focused on the effects of 31 candidate genes and associated SNPs selected from previous G × E studies for childhood maltreatment's association with depression. ResultsAt SNP-level, only the high-multimorbidity Cluster 6 revealed a genome-wide significant SNP rs145772219. At gene-level, MPST and PRH2 were genome-wide significant for the low-multimorbidity Clusters 1 and 3, respectively. Regarding candidate SNPs for G × E interactions, individual SNP results could be replicated for specific clusters. The candidate genes CREB1, DBH, and MTHFR (Cluster 5) as well as TPH1 (Cluster 6) survived multiple testing correction. LimitationsCTS is a short retrospective self-reported measurement. Clusters could be influenced by genetics of individual disorders. ConclusionsThe first G × E GWAS for MDD-related multimorbidity trajectories successfully replicated findings from previous G × E studies related to depression, and revealed risk clusters for the contribution of childhood trauma.

Causal association of metformin treatment with diverse cardiovascular diseases: a Mendelian randomization analysis.

The cardiovascular effects of metformin continue to be a subject of debate within the medical community. The Mendelian randomization (MR) study used data from genome-wide association studies (GWAS) to explore the causal association with six diseases that are associated with bimatoprost treatment and myocardial infarction, chronic heart failure, atrial fibrillation, hypertrophic and dilated cardiomyopathy, and valvular disease. Genome-wide significant single nucleotide polymorphisms (SNPs), that are associated with metformin use were selected as the instrumental variables. To determine the causal relationship between metformin use and various cardiovascular diseases, MR analysis was conducted, employing methods such as Instrumental Variable Weighting (IVW). The IVW analysis demonstrated a positive association between metformin treatment and the risk of myocardial infarction (OR = 22.67, 95% CI 3.22-34.01; P = 0.002). Conversely, metformin treatment exhibited a negative association with the risk of developing valvular disease (OR = 0.98, 95% CI 0.95-1.00; P = 0.046) and hypertrophic cardiomyopathy (OR = 0.01, 95% CI 0.00-0.22; P = 0.016). Multiple test correction found that metformin treatment was causally associated with the risk of both hypertrophic cardiomyopathy (PFDR = 0.048) and myocardial infarction (PFDR = 0.012). The analysis revealed limited heterogeneity in the individual results, absence of pleiotropy evidence, and indications of stability in the findings. The MR study discovered from a genetic standpoint that metformin may lower the risk of hypertrophic cardiomyopathy and valvular heart disease, yet it could elevate the risk of myocardial infarction.

Genetic history of esophageal cancer group in southwestern China revealed by Y‐chromosome STRs and genomic evolutionary connection analysis

Genetic and environmental factors play crucial roles in the development of esophageal cancer (EC) and contribute uniquely or cooperatively to human cancer susceptibility. Sichuan is located in the interior of southwestern China, and the northern part of Sichuan is one of the regions with a high occurrence of EC. However, the factors influencing the high incidence rate of EC in the Sichuan Han Chinese population and its corresponding genetic background and origins are still poorly understood. Here, we utilized genome-wide single nucleotide polymorphisms (SNPs) and Y-chromosome short tandem repeats (Y-STRs) to characterize the genetic structure, connection, and origin of cancer groups and general populations. We generated Y-STR-based haplotype data from 214 Sichuan individuals, including the Han Chinese EC population and a control group of Han Chinese individuals. Our results, obtained from Y-STR-based population statistical methods (analysis of molecular variance (AMOVA), principal component analysis (PCA), and phylogenetic analysis), demonstrated that there was a genetic substructure difference between the EC population in the high-incidence area of northern Sichuan Province and the control population. Additionally, there was a strong genetic relationship between the EC population in the northern Sichuan high-incidence area and those at high risk in both the Fujian and Chaoshan areas. In addition, we obtained high-density SNP data from saliva samples of 60 healthy Han Chinese individuals from three high-prevalence areas of EC in China: Sichuan Nanchong, Fujian Quanzhou, and Henan Xinxiang. As inferred from the allele frequency of SNPs and sharing patterns of haplotype segments, the evolutionary history and admixture events suggested that the Han population from Nanchong in northern Sichuan Province shared a close genetic relationship with the Han populations from Xinxiang in Henan Province and Quanzhou in Fujian Province, both of which are regions with a high prevalence of EC. Our study illuminated the genetic profile and connection of the Northern Sichuan Han population and enriched the genomic resources and features of the Han Chinese populations in China, especially for the Y-STR genetic data of the Han Chinese EC population. Populations living in different regions with high incidences of EC may share similar genetic backgrounds, which offers new insights for further exploring the genetic mechanisms underlying EC.

PM2.5 associates with blood pressure: a Mendelian randomization analysis

ABSTRACT The relationship between fine particulate matter (PM2.5) and blood pressure (BP) is a controversial issue. We conducted a two-sample Mendelian randomization (MR) analysis and identified 58 genome-wide significant single-nucleotide polymorphisms associated with PM2.5 as instrument variables. Inverse-variance weighted (IVW) was used as the primary analysis approach. MR-Egger, weighted median, simple model, and weighted model methods were selected for quality control. We found a significant negative causal association of higher genetically predicted PM2.5 levels with lower systolic BP (SBP), while no causal relationship was identified between PM2.5 and diastolic BP (DBP). For each 1 standard deviation increase in genetically predicted PM2.5 levels, the beta value (95% CI) of SBP was -0.14 (-0.25, -0.03) for IVW (p=0.02), and -0.13 (-0.22, -0.04) for weighted median (p=0.005). Increased PM2.5 concentrations can lead to decreased SBP levels. Our findings provided novel insights into the controversial topic on the causal relationship between PM2.5 and BP.

Genome-wide single nucleotide polymorphisms reveal the genetic diversity and population structure of Creole goats from northern Peru

Goat farming constitutes a significant source of income for farmers in northern Peru. There is currently an absence of information about the genetics of Peruvian Creole goats that would enable us to understand their origins and genetic spread. The objective of this study was to estimate the genetic diversity of Creole goats from northern Peru using SNP markers. This study involved the collection of 192 male Creole goats from three key goat production geographical departments in northern Peru. These goat samples were genotyped using the GGPGoat70k SNP panel. To explore the genetic influence of other breeds on Peruvian Creole goats, our dataset was combined with previously published SNP genotypes. External data set includes multiple breeds genotypes sampled from Argentina, Brazil, Spain, and Alpine breed from Italy, France, and Switzerland. After quality control 52,832 autosomal SNPs were used to assess genetic diversity in the Peruvian goats. For the population structure analysis of the merged data 20,513 common SNPs were used. Estimations for expected heterozygosity (He), observed heterozygosity (Ho), and inbreeding coefficient (FIS) were computed for the Peruvian groups. AMOVA, principal component analysis and ADMIXTURE were conducted to evaluate the population structure in the two data sets, Peru and merged. The results revealed a considerable genetic diversity, with Ho values ranging from 0.40 to 0.41 for the Peruvian sampling groups, and inbreeding coefficient was notably low for Peruvian goat. The population structure analysis demonstrated a distinction (p < 0.05) from other breeds. These findings suggest a level of genetic differentiation of the Peruvian goat population among other breeds, although further research is needed considering samples from other Peruvian areas. We expect this study will contribute to define genetic management strategies to prevent the loss of genetic diversity in Peruvian goat populations and for upcoming advancements in this field.

Comprehensive genome assembly reveals genetic diversity and carcass consumption insights in critically endangered Asian king vultures

The Asian king vulture (AKV), a vital forest scavenger, is facing globally critical endangerment. This study aimed to construct a reference genome to unveil the mechanisms underlying its scavenger abilities and to assess the genetic relatedness of the captive population in Thailand. A reference genome of a female AKV was assembled from sequencing reads obtained from both PacBio long-read and MGI short-read sequencing platforms. Comparative genomics with New World vultures (NWVs) and other birds in the Family Accipitridae revealed unique gene families in AKV associated with retroviral genome integration and feather keratin, contrasting with NWVs’ genes related to olfactory reception. Expanded gene families in AKV were linked to inflammatory response, iron regulation and spermatogenesis. Positively selected genes included those associated with anti-apoptosis, immune response and muscle cell development, shedding light on adaptations for carcass consumption and high-altitude soaring. Using restriction site-associated DNA sequencing (RADseq)-based genome-wide single nucleotide polymorphisms (SNPs), genetic relatedness and inbreeding status of five captive AKVs were determined, revealing high genomic inbreeding in two females. In conclusion, the AKV reference genome was established, providing insights into its unique characteristics. Additionally, the potential of RADseq-based genome-wide SNPs for selecting AKV breeders was demonstrated.

Tasmanian devil (Sarcophilus harrisii) gene flow and source-sink dynamics

Increased access to genetic data has substantially improved how we manage threatened species. The Tasmanian devil (Sarcophilus harrisii) is listed as endangered due to the ongoing threat of a highly contagious cancer, devil facial tumour disease (DFTD), causing more than 80% population reductions. To assist future management interventions (e.g. releases into wild sites) we expanded upon previous studies of gene flow for the devil by assessing more recent and broad-scale patterns. We use genome-wide single nucleotide polymorphisms generated via DArTSeq across 21 devil sites to delineate source-sink dynamics across the species’ range. Our findings revealed gene flow is stronger on the northeast and central regions of Tasmania, with high rates of bidirectional gene flow among central sites. The northwest exhibits weaker connectivity relative to other regions of Tasmania, while gene flow appears to be non-existent between the southwest and other areas. Northeast coastal sites tend to serve as ‘sources’ for inland central sites, whereas gene flow appears restricted to the coastline in the northwest. These results are consistent with genetic structure of devil sites and spatial spread of DFTD, which has yet to arrive in the southwest region of Tasmania. Southwest isolation is probably due to mountain ranges and lack of roadways. Interestingly, some waterbodies did not appear to restrict devil movement among sites. Conversely, areas of high elevation act as apparent barriers, as evidenced by limited gene flow observed between eastern and western sites. Integrating source-sink dynamics into conservation management planning will be crucial in developing effective strategies to safeguard the Tasmanian devil and other threatened species facing similar threats (i.e. disease, habitat loss).

Low genetic variability and high isolation of a post-harvest South American pinniped population as revealed by genome-wide single nucleotide polymorphisms

Abstract Otaria flavescens has been one of the most heavily exploited pinnipeds during the last 200 years, with depletion of ~90% in some colonies. After the prohibition on sealing in South America, populations stabilized except for the Uruguayan population, which showed a constant decrease. The underlying causes of its decline are unknown. This study used genome-wide single nucleotide polymorphisms to assess the variability and connectivity of some of the most overexploited sea lion colonies in the Atlantic Ocean. Our results revealed low allelic richness, nucleotide diversity, and heterozygosity in the Uruguayan population and evidence of complete isolation from the Argentinean populations under study. In contrast, the Patagonian populations showed a high degree of connectivity, which could explain their recovery and high levels of diversity at present. Our study highlights the critical situation of the Uruguayan sea lion population, emphasizing the need for maintaining continuous conservation efforts in the region.

Shared Genetic Determinants of Schizophrenia and Autism Spectrum Disorder Implicate Opposite Risk Patterns: A Genome-Wide Analysis of Common Variants.

The synaptic pruning hypothesis posits that schizophrenia (SCZ) and autism spectrum disorder (ASD) may represent opposite ends of neurodevelopmental disorders: individuals with ASD exhibit an overabundance of synapses and connections while SCZ was characterized by excessive pruning of synapses and a reduction. Given the strong genetic predisposition of both disorders, we propose a shared genetic component, with certain loci having differential regulatory impacts. Genome-Wide single nucleotide polymorphism (SNP) data of European descent from SCZ (N cases = 53386, N controls = 77258) and ASD (N cases = 18381, N controls = 27969) were analyzed. We used genetic correlation, bivariate causal mixture model, conditional false discovery rate method, colocalization, Transcriptome-Wide Association Study (TWAS), and Phenome-Wide Association Study (PheWAS) to investigate the genetic overlap and gene expression pattern. We found a positive genetic correlation between SCZ and ASD (rg = .26, SE = 0.01, P = 7.87e-14), with 11 genomic loci jointly influencing both conditions (conjFDR <0.05). Functional analysis highlights a significant enrichment of shared genes during early to mid-fetal developmental stages. A notable genetic region on chromosome 17q21.31 (lead SNP rs2696609) showed strong evidence of colocalization (PP.H4.abf = 0.85). This SNP rs2696609 is linked to many imaging-derived brain phenotypes. TWAS indicated opposing gene expression patterns (primarily pseudogenes and long noncoding RNAs [lncRNAs]) for ASD and SCZ in the 17q21.31 region and some genes (LRRC37A4P, LINC02210, and DND1P1) exhibit considerable variation in the cerebellum across the lifespan. Our findings support a shared genetic basis for SCZ and ASD. A common genetic variant, rs2696609, located in the Chr17q21.31 locus, may exert differential risk regulation on SCZ and ASD by altering brain structure. Future studies should focus on the role of pseudogenes, lncRNAs, and cerebellum in synaptic pruning and neurodevelopmental disorders.

Highly structured populations of copepods at risk to deep-sea mining: Integration of genomic data with demogenetic and biophysical modelling.

Copepoda is the most abundant taxon in deep-sea hydrothermal vents, where hard substrate is available. Despite the increasing interest in seafloor massive sulphides exploitation, there have been no population genomic studies conducted on vent meiofauna, which are known to contribute over 50% to metazoan biodiversity at vents. To bridge this knowledge gap, restriction-site-associated DNA sequencing, specifically 2b-RADseq, was used to retrieve thousands of genome-wide single-nucleotide polymorphisms (SNPs) from abundant populations of the vent-obligate copepod Stygiopontius lauensis from the Lau Basin. SNPs were used to investigate population structure, demographic histories and genotype-environment associations at a basin scale. Genetic analyses also helped to evaluate the suitability of tailored larval dispersal models and the parameterization of life-history traits that better fit the population patterns observed in the genomic dataset for the target organism. Highly structured populations were observed on both spatial and temporal scales, with divergence of populations between the north, mid, and south of the basin estimated to have occurred after the creation of the major transform fault dividing the Australian and the Niuafo'ou tectonic plate (350 kya), with relatively recent secondary contact events (<20 kya). Larval dispersal models were able to predict the high levels of structure and the highly asymmetric northward low-level gene flow observed in the genomic data. These results differ from most studies conducted on megafauna in the region, elucidating the need to incorporate smaller size when considering site prospecting for deep-sea exploitation of seafloor massive sulphides, and the creation of area-based management tools to protect areas at risk of local extinction, should mining occur.

Genome-wide association study identifies a novel candidate gene for egg production traits in chickens.

Qingyuan partridge chicken is a renowned indigenous yellow broiler breed in China. Egg production traits are important economic traits for chickens. With the decreasing cost of whole genome resequencing, identifying candidate genes with more precision has become possible. In order to identify molecular markers and candidate genes associated with egg production traits, we conducted genome-wide association studies based on the resequencing data of 287 female Qingyuan partridge chickens. For each hen, age at first egg and egg laying rate were recorded and calculated, respectively. With a univariate linear mixed model, we detected one genome-wide significant single nucleotide polymorphism (SNP) and three chromosome-wide significant SNPs associated with egg laying rate. MTA2 is highly likely to be a functional gene for egg laying rate. Our study identifies MTA2 as the first time to be associated with egg laying rate. Findings in our study will advance our understanding of the genetic basis of egg production and have the potential to improve the efficiency of genomic selection in chickens.

Impact of Family History and Germline Genetic Risk Single Nucleotide Polymorphisms on Long-Term Outcomes of Favorable-Risk Prostate Cancer.

Family history and germline genetic risk single nucleotide polymorphisms (SNPs) have been separately shown to stratify lifetime risk of prostate cancer. Here, we evaluate the combined prognostic value of family history of prostate and other related cancers and germline risk SNPs among patients with favorable-risk prostate cancer. A total of 1367 participants from the prospective Health Professionals Follow-up Study diagnosed with low- or favorable intermediate-risk prostate cancer from 1986 to 2017 underwent genome-wide SNP genotyping. Multivariable Cox regression was used to estimate the association between family history, specific germline risk variants, and a 269 SNP polygenic risk score with prostate cancer‒specific death. Family history of prostate, breast, and/or pancreatic cancer was observed in 489 (36%) participants. With median follow-up from diagnosis of 14.9 years, participants with favorable-risk prostate cancer with a positive family history had a significantly higher risk of prostate cancer‒specific death (HR 1.95, 95% CI 1.15-3.32, P = .014) compared to those without any family history. The rs2735839 (19q13) risk allele was associated with prostate cancer‒specific death (HR 1.81 per risk allele, 95% CI 1.04-3.17, P = .037), whereas the polygenic risk score was not. Combined family history and rs2735839 risk allele were each associated with an additive risk of prostate cancer‒specific death (HR 1.78 per risk factor, 95% CI 1.25-2.53, P = .001). Family history of prostate, breast, or pancreatic cancer and/or a 19q13 germline risk allele are associated with an elevated risk of prostate cancer‒specific death among favorable-risk patients. These findings have implications for how family history and germline genetic risk SNPs should be factored into clinical decision-making around favorable-risk prostate cancer.

Heritability and predictive ability for the heterobothriosis resistance and growth performance in the tiger pufferfish Takifugu rubripes fed standard or low fishmeal diets

Heterobothriosis resistance and high growth performance are major concerns in the tiger pufferfish (Takifugu rubripes) farming. Our previous small-scale study suggested the feasibility of genomic selection (GS) for these traits but also identified an undesired genetic correlation between them, which could hamper simultaneous genetic improvement. Here, we sought to validate these initial findings through larger samples and crosses and also explore the potential of GS for the individuals fed a normal fishmeal diet and a low fishmeal (LFM) diet. We produced three populations: one was fed a normal diet (PopFM: 1100 individuals) and two were fed an LFM diet (PopLFM1: 993 individuals; PopLFM2: 1017 individuals); all three populations were artificially infected with Heterobothrium okamotoi, a monogenean parasite causing heterobothriosis. We recorded standard lengths of the fish after artificial infection (SLEAI) and the number of parasites on the gill cavity wall (HC) and genotyped genome-wide single nucleotide polymorphisms (SNPs) using genotyping by random amplicon sequencing, direct (GRAS-Di) technology for each sample: PopFM, 12 K SNPs; PopLFM1, 19 K SNPs; PopLFM2, 21 K SNPs. Heritability and prediction ability for each trait were estimated for each population using genomic best linear unbiased prediction (GBLUP), Bayes C, and Bayesian ridge regression methods. SNP effects on each trait were estimated by Bayes C, Bayesian ridge regression, and genome-wide association analyses (GWAS) using the same data sets. We observed moderate to high levels of heritability for HC (0.16–0.57) and SLEAI (0.50–0.77), with predictive ability of 0.25–0.55 for HC and of 0.51–0.73 for SLEAI. These values were slightly larger in LFM diet groups. As reported previously, a weak undesired genetic correlation was observed between the traits, and the correlation was more pronounced in the LFM groups. In conclusion, our results indicate that genetic gain in growth and parasite resistance traits can be expected for the tiger pufferfish through GS, especially for individuals reared on an LFM diet. However, close attention should be paid to the extent of the undesirable genetic correlation between the two traits that may be amplified in the populations fed an LFM diet.

The role of historical biogeography in shaping colour morph diversity in the common wall lizard.

The maintenance of polymorphisms often depends on multiple selective forces, but less is known on the role of stochastic or historical processes in maintaining variation. The common wall lizard (Podarcis muralis) is a colour polymorphic species in which local colour morph frequencies are thought to be modulated by natural and sexual selection. Here, we used genome-wide single-nucleotide polymorphism data to investigate the relationships between morph composition and population biogeography at a regional scale, by comparing morph composition with patterns of genetic variation of 54 populations sampled across the Pyrenees. We found that genetic divergence was explained by geographic distance but not by environmental features. Differences in morph composition were associated with genetic and environmental differentiation, as well as differences in sex ratio. Thus, variation in colour morph frequencies could have arisen via historical events and/or differences in the permeability to gene flow, possibly shaped by the complex topography and environment. In agreement with this hypothesis, colour morph diversity was positively correlated with genetic diversity and rates of gene flow and inversely correlated with the likelihood of the occurrence of bottlenecks. Concurrently, we did not find conclusive evidence for selection in the two colour loci. As an illustration of these effects, we observed that populations with higher proportions of the rarer yellow and yellow-orange morphs had higher genetic diversity. Our results suggest that processes involving a decay in overall genetic diversity, such as reduced gene flow and/or bottleneck events have an important role in shaping population-specific morph composition via non-selective processes.

Effects of genetically predicted posttraumatic stress disorder on autoimmune phenotypes

Observational studies suggest that posttraumatic stress disorder (PTSD) increases risk for various autoimmune diseases. Insights into shared biology and causal relationships between these diseases may inform intervention approaches to PTSD and co-morbid autoimmune conditions. We investigated the shared genetic contributions and causal relationships between PTSD, 18 autoimmune diseases, and 3 immune/inflammatory biomarkers. Univariate MiXeR was used to contrast the genetic architectures of phenotypes. Genetic correlations were estimated using linkage disequilibrium score regression. Bi-directional, two-sample Mendelian randomization (MR) was performed using independent, genome-wide significant single nucleotide polymorphisms; inverse variance weighted and weighted median MR estimates were evaluated. Sensitivity analyses for uncorrelated (MR PRESSO) and correlated horizontal pleiotropy (CAUSE) were also performed. PTSD was considerably more polygenic (10,863 influential variants) than autoimmune diseases (median 255 influential variants). However, PTSD evidenced significant genetic correlation with nine autoimmune diseases and three inflammatory biomarkers. PTSD had putative causal effects on autoimmune thyroid disease (p = 0.00009) and C-reactive protein (CRP) (p = 4.3 × 10−7). Inferences were not substantially altered by sensitivity analyses. Additionally, the PTSD-autoimmune thyroid disease association remained significant in multivariable MR analysis adjusted for genetically predicted inflammatory biomarkers as potential mechanistic pathway variables. No autoimmune disease had a significant causal effect on PTSD (all p values > 0.05). Although causal effect models were supported for associations of PTSD with CRP, shared pleiotropy was adequate to explain a putative causal effect of CRP on PTSD (p = 0.18). In summary, our results suggest a significant genetic overlap between PTSD, autoimmune diseases, and biomarkers of inflammation. PTSD has a putative causal effect on autoimmune thyroid disease, consistent with existing epidemiologic evidence. A previously reported causal effect of CRP on PTSD is potentially confounded by shared genetics. Together, results highlight the nuanced links between PTSD, autoimmune disorders, and associated inflammatory signatures, and suggest the importance of targeting related pathways to protect against disease and disability.

Large-scale genome-wide SNP analysis reveals the rugged (and ragged) landscape of global ancestry, phylogeny, and demographic history in chicken breeds.

The worldwide chicken gene pool encompasses a remarkable, but shrinking, number of divergently selected breeds of diverse origin. This study was a large-scale genome-wide analysis of the landscape of the complex molecular architecture, genetic variability, and detailed structure among 49 populations. These populations represent a significant sample of the world's chicken breeds from Europe (Russia, Czech Republic, France, Spain, UK, etc.), Asia (China), North America (USA), and Oceania (Australia). Based on the results of breed genotyping using the Illumina 60K single nucleotide polymorphism (SNP) chip, a bioinformatic analysis was carried out. This included the calculation of heterozygosity/homozygosity statistics, inbreeding coefficients, and effective population size. It also included assessment of linkage disequilibrium and construction of phylogenetic trees. Using multidimensional scaling, principal component analysis, and ADMIXTURE-assisted global ancestry analysis, we explored the genetic structure of populations and subpopulations in each breed. An overall 49-population phylogeny analysis was also performed, and a refined evolutionary model of chicken breed formation was proposed, which included egg, meat, dual-purpose types, and ambiguous breeds. Such a large-scale survey of genetic resources in poultry farming using modern genomic methods is of great interest both from the viewpoint of a general understanding of the genetics of the domestic chicken and for the further development of genomic technologies and approaches in poultry breeding. In general, whole genome SNP genotyping of promising chicken breeds from the worldwide gene pool will promote the further development of modern genomic science as applied to poultry.

Genome-wide association studies of egg production traits by whole genome sequencing of Laiwu Black chicken

Compared to high-yield commercial laying hens, Chinese indigenous chicken breeds have poor egg laying capacity due to the lack of intensive selection. However, as these breeds have not undergone systematic selection, it is possible that there is a greater abundance of genetic variations related to egg laying traits. In this study, we assessed 5 egg number (EN) traits at different stages of the egg-laying period: EN1 (from the first egg to 23 wk), EN2 (from 23 to 35 wk), EN3 (from 35 to 48 wk), EN4 (from the first egg to 35 wk), and EN5 (from the first egg to 48 wk). To investigate the molecular mechanisms underlying egg number traits in a Chinese local chicken breed, we conducted a genome-wide association study (GWAS) using data from whole-genome sequencing (WGS) of 399 Laiwu Black chickens. We obtained a total of 3.01 Tb of raw data with an average depth of 7.07 × per individual. A total of 86 genome-wide suggestive or significant single-nucleotide polymorphisms (SNP) contained within a set of 45 corresponding candidate genes were identified and found to be associated with stages EN1-EN5. The genes vitellogenin 2 (VTG2), lipase maturation factor 1 (LMF1), calcium voltage-gated channel auxiliary subunit alpha2delta 3 (CACNA2D3), poly(A) binding protein cytoplasmic 1 (PABPC1), programmed cell death 11 (PDCD11) and family with sequence similarity 213 member A (FAM213A) can be considered as the candidate genes associated with egg number traits, due to their reported association with animal reproduction traits. Noteworthy, results suggests that VTG2 and PDCD11 are not only involved in the regulation of EN3, but also in the regulation of EN5, implies that VTG2 and PDCD11 have a significant influence on egg production traits. Our study offers valuable genomic insights into the molecular genetic mechanisms that govern egg number traits in a Chinese indigenous egg-laying chicken breed. These findings have the potential to enhance the egg-laying performance of chickens.

Phylogenetic relationships of Quercus species (Fagaceae) in Vietnam based on multiplexed inter simple sequence repeat genotyping by sequencing

Quercus is the largest genus of the family Fagaceae in the world and the second-largest in Vietnam, with approximately 50 species. However, the phylogeny study of the Quercus genus in Vietnam has yet to be thoroughly explored. In this study, we utilize the genome-wide single-nucleotide polymorphisms (SNPs) data obtained through Multiplexed Inter Simple Sequence Repeat Genotyping by sequencing (MIG-seq) to explore the phylogenetic relationships among Quercus species in Vietnam. The results of this study reveal that all Quercus species in Vietnam belong to subgenus Cerris and the phylogenetic analysis strongly supports the recognition of two infrageneric sections: Quercus and Ilex section for the Vietnamese Quercus. These results also confirm that the monophyly of Quercus in Vietnam is strongly supported by both morphological and molecular data. The results of this study also align entirely with previous research, indicating that Quercus species in Asia fall under the subgenus Cerris, comprising two sections: section Ilex and section Cyclobalanopsis.