Search

Pubertal timing varies considerably and is associated with later health outcomes. We performed multi-ancestry genetic analyses on ~800,000 women, identifying 1,080 signals for age at menarche. Collectively, these explained 11% of trait variance in an independent sample. Women at the top and bottom 1% of polygenic risk exhibited ~11 and ~14-fold higher risks of delayed and precocious puberty, respectively. We identified several genes harboring rare loss-of-function variants in ~200,000 women, including variants in ZNF483, which abolished the impact of polygenic risk. Variant-to-gene mapping approaches and mouse gonadotropin-releasing hormone neuron RNA sequencing implicated 665 genes, including an uncharacterized G-protein-coupled receptor, GPR83, which amplified the signaling of MC3R, a key nutritional sensor. Shared signals with menopause timing at genes involved in DNA damage response suggest that the ovarian reserve might signal centrally to trigger puberty. We also highlight body size-dependent and independent mechanisms that potentially link reproductive timing to later life disease.

Genomic time series from experimental evolution studies and ancient DNA datasets offer us a chance to directly observe the interplay of various evolutionary forces. We show how the genome-wide variance in allele frequency change between two time points can be decomposed into the contributions of gene flow, genetic drift, and linked selection. In closed populations, the contribution of linked selection is identifiable because it creates covariances between time intervals, and genetic drift does not. However, repeated gene flow between populations can also produce directionality in allele frequency change, creating covariances. We show how to accurately separate the fraction of variance in allele frequency change due to admixture and linked selection in a population receiving gene flow. We use two human ancient DNA datasets, spanning around 5,000 y, as time transects to quantify the contributions to the genome-wide variance in allele frequency change. We find that a large fraction of genome-wide change is due to gene flow. In both cases, after correcting for known major gene flow events, we do not observe a signal of genome-wide linked selection. Thus despite the known role of selection in shaping long-term polymorphism levels, and an increasing number of examples of strong selection on single loci and polygenic scores from ancient DNA, it appears to be gene flow and drift, and not selection, that are the main determinants of recent genome-wide allele frequency change. Our approach should be applicable to the growing number of contemporary and ancient temporal population genomics datasets.

Idiopathic epilepsy (IE) has been known to be inherited in the Belgian Tervuren for many decades. Risk genotypes for IE in this breed have recently been identified on Canis familiaris chromosomes (CFA) 14 and 37. In the current study, the allele frequencies of these loci were analyzed to determine whether dog breeders had employed a purposeful selection against IE, leading to a reduction in risk-associated allele frequency within the breed over time. The allele frequencies of two generational groupings of Belgian Tervuren with and without IE were compared. Allele frequencies for risk-associated alleles on CFA14 were unchanged between 1985 and 2015, whereas those on CFA37 increased during that time in the control population (p < 0.05). In contrast, dogs with IE showed a decrease (p < 0.05) in the IE risk-associated allele frequency at the CFA37 locus. Seizure prevalence in the Belgian Tervuren appears to be increasing. These results suggest that, despite awareness that IE is inherited, selection against IE has not been successful.

Understanding and treating human diseases require valid animal models. Leveraging the genetic diversity in rhesus macaque populations across eight primate centers in the United States, we conduct targeted-sequencing on 1845 individuals for 374 genes linked to inherited human retinal and neurodevelopmental diseases. We identify over 47,000 single nucleotide variants, a substantial proportion of which are shared with human populations. By combining rhesus and human allele frequencies with established variant prediction methods, we develop a machine learning-based score that outperforms established methods in predicting missense variant pathogenicity. Remarkably, we find a marked number of loss-of-function variants and putative deleterious variants, which may lead to the development of rhesus disease models. Through phenotyping of macaques carrying a pathogenic OPA1:p.A8S variant, we identify a genetic model of autosomal dominant optic atrophy. Finally, we present a public website housing variant and genotype data from over two thousand rhesus macaques.

Genome‐wide technologies open up new possibilities to clarify questions on genetic structure and phylogeographic history of taxa previously studied with microsatellite loci and mitochondrial sequences. Here, we used 736 individual red deer (Cervus elaphus) samples genotyped at 35,701 single nucleotide polymorphism loci (SNPs) to assess the population structure of the species throughout Europe. The results identified 28 populations, with higher degrees of genetic distinction in peripheral compared to mainland populations. Iberian red deer show high genetic differentiation, with lineages in Western and Central Iberia maintaining their distinctiveness, which supports separate refugial ranges within Iberia along with little recent connection between Iberian and the remaining Western European populations. The Norwegian population exhibited the lowest variability and the largest allele frequency differences from mainland European populations, compatible with a history of bottlenecks and drift during post‐glacial colonization from southern refugia. Scottish populations showed high genetic distance from the mainland but high levels of diversity. Hybrid zones were found between Eastern and Western European lineages in Central Europe as well as in the Pyrenees, where red deer from France are in close contact with Iberian red deer. Anthropogenic restocking has promoted the Pyrenean contact zone, admixture events in populations on the Isle of Rum and in the Netherlands, and at least partly the admixture of the two main lineages in central‐eastern Europe. Our analysis enabled detailed resolution of population structure of a large mammal widely distributed throughout Europe and contributes to resolving the evolutionary history, which can also inform conservation and management policies. [ABSTRACT FROM AUTHOR]

Denitrification is a crucial process in the global nitrogen cycle, in which two functionally equivalent genes, nirS and nirK, catalyse the critical reaction and are usually used as marker genes. The nirK gene can function independently, whereas nirS requires additional genes to encode nitrite reductase and is more sensitive to environmental factors than nirK. However, the ecological differentiation mechanisms of those denitrifying microbial communities and their adaptation strategies to environmental stresses remain unclear. Here, we conducted metagenomic analysis for sediments and bioreactor samples from Lake Donghu, China. We found that nirS‐type denitrifying communities had a significantly lower horizontal gene transfer frequency than that of nirK‐type denitrifying communities, and nirS gene phylogeny was more congruent with taxonomy than that of nirK gene. Metabolic reconstruction of metagenome‐assembled genomes further revealed that nirS‐type denitrifying communities have robust metabolic systems for energy conservation, enabling them to survive under environmental stresses. Nevertheless, nirK‐type denitrifying communities seemed to adapt to oxygen‐limited environments with the ability to utilize various carbon and nitrogen compounds. Thus, this study provides novel insights into the ecological differentiation mechanism of nirS and nirK‐type denitrifying communities, as well as the regulation of the global nitrogen cycle and greenhouse gas emissions. [ABSTRACT FROM AUTHOR]

AbstractObesity is a common public health problem associated with serious, life-threatening complications. MicroRNAs (miRs) have modulating roles in the immune and inflammatory systems. Therefore, this study aimed to analyze the relationship between miR-146a and morbid obesity (MO) in a Turkish population. In this study, a total of 258 subjects (110 patients with MO and 148 controls) were genotyped by the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method to analyze miR-146a rs2910164. Then, we examined the patients as males and females separately. The results of the analyses were evaluated for statistical significance. There was a significant difference in genotype and allele frequencies of miR-146a rs2910164 between patients with MO and control individuals. miR-146a rs2910164 CC genotype and C allele were shown to increase in the MO patients’ group compared to the control group (p = 0.000, p = 0.000, respectively). Also, the C allele was higher in both female and male patients compared to controls (p = 0.000, p = 0.000, respectively). High differences were also observed when the patients and the controls were compared according to CC versus GG + GC and GG versus GC + CC (p = 0.000, p = 0.000, respectively). A significant difference was found between the female/male patients and the female/male controls in terms of GG + GC versus CC (p = 0.000, p = 0.000, respectively). To the best of our knowledge, this is the first study to investigate the relationship between this variant and MO in Turkey. Our results showed that miR-146a rs2910164 is a valuable biomarker that can be used to distinguish between patients with MO and the healthy population. The findings can be extended by increasing the sample sizes with diverse ethnicities. [ABSTRACT FROM AUTHOR]

Single nucleotide polymorphism (SNP) interactions are the key to improving polygenic risk scores. Previous studies reported several significant SNP–SNP interaction pairs that shared a common SNP to form a cluster, but some identified pairs might be false positives. This study aims to identify factors associated with the cluster effect of false positivity and develop strategies to enhance the accuracy of SNP–SNP interactions. The results showed the cluster effect is a major cause of false-positive findings of SNP–SNP interactions. This cluster effect is due to high correlations between a causal pair and null pairs in a cluster. The clusters with a hub SNP with a significant main effect and a large minor allele frequency (MAF) tended to have a higher false-positive rate. In addition, peripheral null SNPs in a cluster with a small MAF tended to enhance false positivity. We also demonstrated that using the modified significance criterion based on the 3 p-value rules and the bootstrap approach (3pRule + bootstrap) can reduce false positivity and maintain high true positivity. In addition, our results also showed that a pair without a significant main effect tends to have weak or no interaction. This study identified the cluster effect and suggested using the 3pRule + bootstrap approach to enhance SNP–SNP interaction detection accuracy. [ABSTRACT FROM AUTHOR]

Introduction: Spillover events of Mycoplasma ovipneumoniae have devastating effects on the wild sheep populations. Multilocus sequence typing (MLST) is used to monitor spillover events and the spread of M. ovipneumoniae between the sheep populations. Most studies involving the typing of M. ovipneumoniae have used Sanger sequencing. However, this technology is time-consuming, expensive, and is not well suited to efficient batch sample processing. Methods: Our study aimed to develop and validate an MLST workflow for typing of M. ovipneumoniae using Nanopore Rapid Barcoding sequencing and multiplex polymerase chain reaction (PCR). We compare the workflow with Nanopore Native Barcoding library preparation and Illumina MiSeq amplicon protocols to determine the most accurate and cost-effective method for sequencing multiplex amplicons. A multiplex PCR was optimized for four housekeeping genes of M. ovipneumoniae using archived DNA samples (N = 68) from nasal swabs. Results: Sequences recovered from Nanopore Rapid Barcoding correctly identified all MLST types with the shortest total workflow time and lowest cost per sample when compared with Nanopore Native Barcoding and Illumina MiSeq methods. Discussion: Our proposed workflow is a convenient and effective method for strain typing of M. ovipneumoniae and can be applied to other bacterial MLST schemes. The workflow is suitable for diagnostic settings, where reduced hands-on time, cost, and multiplexing capabilities are important. [ABSTRACT FROM AUTHOR]

Cultivated potato, Solanum tuberosum L., is considered an autotetraploid with 12 chromosomes with four homologous phases. However, recent evidence found that, due to frequent large phase deletions in the genome, gene ploidy is not constant across the genome. The elite cultivar "Otava" was found to have an average gene copy number of 3.2 across all loci. Breeding programs for elite potato cultivars rely increasingly on genomic prediction tools for selection breeding and elucidation of quantitative trait loci underpinning trait genetic variance. These are typically based on anonymous single nucleotide polymorphism (SNP) markers, which are usually called from, for example, SNP array or sequencing data using a tetraploid model. In this study, we analyzed the impact of using whole genome markers genotyped as either tetraploid or observed allele frequencies from genotype-by-sequencing data on single-trait additive genomic best linear unbiased prediction (GBLUP) genomic prediction (GP) models and single-marker regression genome-wide association studies of potato to evaluate the implications of capturing varying ploidy on the statistical models employed in genomic breeding. A panel of 762 offspring of a diallel cross of 18 parents of elite breeding material was used for modeling. These were genotyped by sequencing and phenotyped for five key performance traits: chipping quality, length/width ratio, senescence, dry matter content, and yield. We also estimated the read coverage required to confidently discriminate between a heterozygous triploid and tetraploid state from simulated data. It was found that using a tetraploid model neither impaired nor improved genomic predictions compared to using the observed allele frequencies that account for true marker ploidy. In genome-wide associations studies (GWAS), very minor variations of both signal amplitude and number of SNPs supporting both minor and major quantitative trait loci (QTLs) were observed between the two data sets. However, all major QTLs were reproducible using both data sets. [ABSTRACT FROM AUTHOR]

Background: Ventricular septal defect (VSD) is the most common congenital heart disease. Although a small number of genes associated with VSD have been found, the genetic factors of VSD remain unclear. In this study, we evaluated the association of 10 candidate single nucleotide polymorphisms (SNPs) with isolated VSD in a population from Southwest China. Methods: Based on the results of 34 congenital heart disease whole-exome sequencing and 1000 Genomes databases, 10 candidate SNPs were selected. A total of 618 samples were collected from the population of Southwest China, including 285 VSD samples and 333 normal samples. Ten SNPs in the case group and the control group were identified by SNaPshot genotyping. The chi-square (χ2) test was used to evaluate the relationship between VSD and each candidate SNP. The SNPs that had significant P value in the initial stage were further analysed using linkage disequilibrium, and haplotypes were assessed in 34 congenital heart disease whole-exome sequencing samples using Haploview software. The bins of SNPs that were in very strong linkage disequilibrium were further used to predict haplotypes by Arlequin software. ViennaRNA v2.5.1 predicted the haplotype mRNA secondary structure. We evaluated the correlation between mRNA secondary structure changes and ventricular septal defects. Results: The χ2 results showed that the allele frequency of FLT4 rs383985 (P = 0.040) was different between the control group and the case group (P < 0.05). FLT4 rs3736061 (r2 = 1), rs3736062 (r2 = 0.84), rs3736063 (r2 = 0.84) and FLT4 rs383985 were in high linkage disequilibrium (r2 > 0.8). Among them, rs3736061 and rs3736062 SNPs in the FLT4 gene led to synonymous variations of amino acids, but predicting the secondary structure of mRNA might change the secondary structure of mRNA and reduce the free energy. Conclusions: These findings suggest a possible molecular pathogenesis associated with isolated VSD, which warrants investigation in future studies. [ABSTRACT FROM AUTHOR]

Genetic variability persists across diverse populations, and it may impact the characterization of heritable diseases in different ancestral groups. Cystinosis is a metabolic disease caused by pathogenic variants in the CTNS gene causing the cellular accumulation of cystine. We attempted to assess the currently poorly characterized prevalence of cystinosis by employing a population genetics methodology. However, we encountered a significant challenge due to genetic variations across different populations, and the consideration of potential disparities in access to healthcare made our results inconclusive. Pathogenic CTNS variants were identified in a representative global population cohort using The Human Gene Mutation Database (HGMD) and the 1000 Genomes (1 KG) database. The c.124G>A (p.Val42Ile) variant was reported to be pathogenic based on an observation in the white population presenting with atypical phenotypes, but it would be reclassified as benign in the African ancestral group if applying the ACMG allele frequency guideline due to its high allele frequency specifically in this population. Inclusion or exclusion of this c.124G>A (p.Val42Ile) variant results in a significant change in estimated disease prevalence, which can impact the diagnosis and treatment of affected patients with a broad range of phenotypic presentations. This observation led us to postulate that pathogenic manifestations of the disease may be underdiagnosed due to variable expressivity and systemic inequities in access to care, specifically in the African subpopulation. We call for a more cautious and inclusive approach to achieve more equitable care across diverse populations. [ABSTRACT FROM AUTHOR]

In marine species with large populations and high dispersal potential, large-scale genetic differences and clinal trends in allele frequency can provide insight into the evolutionary processes that shape diversity. Lumpfish, Cyclopterus lumpus , is found throughout the North Atlantic and has traditionally been harvested for roe and more recently used as a cleaner fish in salmon aquaculture. We used a 70 K SNP array to evaluate trans-Atlantic differentiation, genetic structuring, and clinal variation across the North Atlantic. Basin-scale structuring between the Northeast and Northwest Atlantic was significant, with enrichment for loci associated with developmental/mitochondrial function. We identified a putative structural variant on chromosome 2, likely contributing to differentiation between Northeast and Northwest Atlantic Lumpfish, and consistent with post-glacial trans-Atlantic secondary contact. Redundancy analysis identified climate associations both in the Northeast (N  = 1269 loci) and Northwest (N  = 1637 loci), with 103 shared loci between them. Clinal patterns in allele frequencies were observed in some loci (15%—Northwest and 5%—Northeast) of which 708 loci were shared and involved with growth, developmental processes, and locomotion. The combined evidence of trans-Atlantic differentiation, environmental associations, and clinal loci, suggests that both regional and large-scale potentially-adaptive population structuring is present across the North Atlantic. [ABSTRACT FROM AUTHOR]

In aquaculture, sterile triploids are commonly used for production as sterility gives them potential gains in growth, yields, and quality. However, they cannot be reproduced, and DNA parentage assignment to their diploid or tetraploid parents is required to estimate breeding values for triploid phenotypes. No publicly available software has the ability to assign triploids to their parents. Here, we updated the R package APIS to support triploids induced from diploid parents. First, we created new exclusion and likelihood tables that account for the double allelic contribution of the dam and the recombination that can occur during female meiosis. As the effective recombination rate of each marker with the centromere is usually unknown, we set it at 0.5 and found that this value maximizes the assignment rate even for markers with high or low recombination rates. The number of markers needed for a high true assignment rate did not strongly depend on the proportion of missing parental genotypes. The assignment power was however affected by the quality of the markers (minor allele frequency, call rate). Altogether, 96–192 SNPs were required to have a high parentage assignment rate in a real rainbow trout dataset of 1,232 triploid progenies from 288 parents. The likelihood approach was more efficient than exclusion when the power of the marker set was limiting. When more markers were used, exclusion was more advantageous, with sensitivity reaching unity, very low false discovery rate (<0.01), and excellent specificity (0.96–0.99). Thus, APIS provides an efficient solution to assign triploids to their diploid parents. [ABSTRACT FROM AUTHOR]

Rice landraces from the fragile riverine ecosystems surrounding Majuli islands in Assam are pivotal for local livelihoods, as well as for their potential in developing flood tolerant varieties. We characterized 87 landraces from Majuli, Dhemaji, and North Lakhimpur districts in Assam by assessing 36 agro‐morphological traits, submergence tolerance, and microsatellite diversity. Multivariate analyses revealed that grain and leaf characteristics, plant height, and tillering are important traits differentiating the cultivar types. Sixty‐six SSR markers generated 271 alleles with an average of 4.1 alleles per marker with population structuring in two clusters having an allele frequency divergence of 0.124. The Bao or deep‐water landraces displayed the highest phenotypic and genetic diversity, with potential genetic links to the aus and wild species. Sali (winter rice) and bora (sticky rice) landraces were predominantly indica type. Nearly 47% of the germplasm exhibited tolerance to two‐week submergence at seedling stage. Surveying SUB1A and SNORKEL genes revealed the presence of SUB1A‐1 and SNORKEL genes in 64% and 57% of the tolerant accessions, respectively. This study highlights the untapped potential of rice landraces in a fragile ecological region, providing insights into their agronomical traits, submergence tolerance, and genetic diversity. The findings lay the foundation for improved germplasm management and developing flood‐resilient rice varieties, crucial for sustainable agriculture in riverine ecosystems. [ABSTRACT FROM AUTHOR]

Blood cell-free DNA (cfDNA) can be a new reliable tool for detecting epidermal growth factor receptor (EGFR) mutations in non-small cell lung cancer (NSCLC) patients. However, the currently reported cfDNA assays have a limited role in detecting drug-resistant mutations due to their deficiencies in sensitivity, stability, or mutation detection rate.We developed an Archaeoglobus fulgidus-derived flap endonuclease (Afu FEN)-based DNA-enhanced amplification system of mutated cfDNA by designing a pair of hairpin probes to anneal with wild-type cfDNA to form two 5′-flaps, allowing for the specific cleavage of wild-type cfDNA by Afu FEN. When the dominant wild-type somatic cfDNA fragments were cleaved by structure-recognition-specific Afu FEN, the proportion of mutated cfDNA in the reaction system was greatly enriched. As the amount of mutated cfDNA in the system was further increased by PCR amplification, the mutation status could be easily detected through first-generation sequencing.In a mixture of synthetic wild-type and T790M EGFR DNA fragments, our new assay still could detect T790M mutation at the fg level with remarkably high sensitivity. We also tested its performance in detecting low variant allele frequency (VAF) mutations in clinical samples from NSCLC patients. The plasma cfDNA samples with low VAF (0.1 and 0.5 %) could be easily detected by DNA-enhanced amplification.This system with enhanced amplification of mutated cfDNA is an effective tool used for the early screening and individualized targeted therapy of NSCLC by providing a rapid, sensitive, and economical way for the detection of drug-resistant mutations in tumors. [ABSTRACT FROM AUTHOR]

• Baseline plasma inflammatory markers were associated with acamprosate treatment outcomes. • Lower baseline IL17RB expression was associated with alcohol relapse. • IL17RB genetic variants might contribute to acamprosate treatment outcomes. Acamprosate is a Food and Drug Administration (FDA) approved medication for the treatment of alcohol use disorder (AUD). However, only a subset of patients achieves optimal treatment outcomes. Currently, no biological measures are utilized to predict response to acamprosate treatment. We applied our established pharmaco-omics informed genomics strategy to identify potential biomarkers associated with acamprosate treatment response. Specifically, our previous open-label acamprosate clinical trial recruited 442 patients with AUD who were treated with acamprosate for three months. We first performed proteomics using baseline plasma samples to identify potential biomarkers associated with acamprosate treatment outcomes. Next, we applied our established "proteomics-informed genome-wide association study (GWAS)" research strategy, and identified 12 proteins, including interleukin-17 receptor B (IL17RB), associated with acamprosate treatment response.​ A GWAS for IL17RB concentrations identified several genome-wide significant signals. Specifically, the top hit single nucleotide polymorphism (SNP) rs6801605 with a minor allele frequency of 38% in the European American population mapped 4 kilobase (Kb) upstream of IL17RB , and intron 1 of the choline dehydrogenase (CHDH) gene on chromosome 3 (p: 4.8E-20). The variant genotype (AA) for the SNP rs6801605 was associated with lower IL17RB protein expression. In addition, we identified a series of genetic variants in IL17RB that were associated with acamprosate treatment outcomes. Furthermore, the variant genotypes for all of those IL17RB SNPs were protective for alcohol relapse. Finally, we demonstrated that the basal level of mRNA expression of IL17RB was inversely correlated with those of nuclear factor-κB (NF-κB) subunits, and a significantly higher expression of NF-κB subunits was observed in AUD patients who relapsed to alcohol use. In summary, this study illustrates that IL17RB genetic variants might contribute to acamprosate treatment outcomes. This series of studies represents an important step toward generating functional hypotheses that could be tested to gain insight into mechanisms underlying acamprosate treatment response phenotypes. (The ClinicalTrials.gov Identifier: NCT00662571). [ABSTRACT FROM AUTHOR]

We consider a single genetic locus with two alleles A 1 and A 2 in a large haploid population. The locus is subject to selection and two-way, or recurrent, mutation. Assuming the allele frequencies follow a Wright–Fisher diffusion and have reached stationarity, we describe the asymptotic behaviors of the conditional gene genealogy and the latent mutations of a sample with known allele counts, when the count n 1 of allele A 1 is fixed, and when either or both the sample size n and the selection strength | α | tend to infinity. Our study extends previous work under neutrality to the case of non-neutral rare alleles, asserting that when selection is not too strong relative to the sample size, even if it is strongly positive or strongly negative in the usual sense (α → − ∞ or α → + ∞), the number of latent mutations of the n 1 copies of allele A 1 follows the same distribution as the number of alleles in the Ewens sampling formula. On the other hand, very strong positive selection relative to the sample size leads to neutral gene genealogies with a single ancient latent mutation. We also demonstrate robustness of our asymptotic results against changing population sizes, when one of | α | or n is large. [ABSTRACT FROM AUTHOR]

Preclinical evidence implicating cannabinoid receptor 2 (CB2) in various diseases has led researchers to question whether CB2 genetics influence aetiology or progression. Associations between conditions and genetic loci are often studied via single nucleotide polymorphism (SNP) prevalence in case versus control populations. In the CNR2 coding exon, ~36 SNPs have high overall population prevalence (minor allele frequencies [MAF] ~37%), including non‐synonymous SNP (ns‐SNP) rs2501432 encoding CB2 63Q/R. Interspersed are ~27 lower frequency SNPs, four being ns‐SNPs. CNR2 introns also harbour numerous SNPs. This review summarises CB2 ns‐SNP molecular pharmacology and evaluates evidence from ~70 studies investigating CB2 genetic variants with proposed linkage to disease. Although CNR2 genetic variation has been associated with a wide variety of conditions, including osteoporosis, immune‐related disorders, and mental illnesses, further work is required to robustly validate CNR2 disease links and clarify specific mechanisms linking CNR2 genetic variation to disease pathophysiology and potential drug responses. [ABSTRACT FROM AUTHOR]

Gastric cancer (GC) is one of the most heterogeneous tumors. However, research on normal tissue adjacent to the tumor (NAT) is very limited. We performed multi‐regional omics sequencing on 150 samples to assess the genetic basis and immune microenvironment in NAT and matched primary tumor or lymph node metastases. NATs demonstrated different mutated genes compared with GC, and NAT genomes underwent independent evolution with low variant allele frequency. Mutation profiles were predominated by aging and smoking‐associated signatures in NAT instead of signatures associated with genetic instability. Although the immune microenvironment within NATs shows substantial intra‐patient heterogeneity, the proportion of shared TCR clones among NATs is five times higher than that of tumor regions. These findings support the notion that subclonal expansion is not pronounced in NATs. We also demonstrated remarkable intra‐patient heterogeneity of GCs and revealed heterogeneity of focal amplification of CD274 (encoding PD‐L1) that leads to differential expression. Finally, we identified that monoclonal seeding is predominant in GC, which is followed by metastasis‐to‐metastasis dissemination in individual lymph nodes. These results provide novel insights into GC carcinogenesis. © 2024 The Pathological Society of Great Britain and Ireland. [ABSTRACT FROM AUTHOR]

Background: The transition from vegetative to reproductive growth is a key factor in yield maximization. Sesame (Sesamum indicum), an indeterminate short-day oilseed crop, is rapidly being introduced into new cultivation areas. Thus, decoding its flowering mechanism is necessary to facilitate adaptation to environmental conditions. In the current study, we uncover the effect of day-length on flowering and yield components using F 2 populations segregating for previously identified quantitative trait loci (Si_DTF QTL) confirming these traits. Results: Generally, day-length affected all phenotypic traits, with short-day preceding days to flowering and reducing yield components. Interestingly, the average days to flowering required for yield maximization was 50 to 55 days, regardless of day-length. In addition, we found that Si_DTF QTL is more associated with seed-yield and yield components than with days to flowering. A bulk-segregation analysis was applied to identify additional QTL differing in allele frequencies between early and late flowering under both day-length conditions. Candidate genes mining within the identified major QTL intervals revealed two flowering-related genes with different expression levels between the parental lines, indicating their contribution to sesame flowering regulation. Conclusions: Our findings demonstrate the essential role of flowering date on yield components and will serve as a basis for future sesame breeding. [ABSTRACT FROM AUTHOR]

Introduction: Genotyping large-scale gene bank collections requires an appropriate sampling strategy to represent the diversity within and between accessions. Methods: A panel of 44 common bean (Phaseolus vulgaris L.) landraces from the Alliance Bioversity and The Alliance of Bioversity International and the International Center for Tropical Agriculture (CIAT) gene bank was genotyped with DArTseq using three sampling strategies: a single plant per accession, 25 individual plants per accession jointly analyzed after genotyping (in silico-pool), and by pooling tissue from 25 individual plants per accession (seq-pool). Sampling strategies were compared to assess the technical aspects of the samples, the marker information content, and the genetic composition of the panel. Results: The seq-pool strategy resulted inmore consistent DNA libraries for quality and call rate, although with fewer polymorphic markers (6,142 single-nucleotide polymorphisms) than the in silico-pool (14,074) or the single plant sets (6,555). Estimates of allele frequencies by seq-pool and in silico-pool genotyping were consistent, but the results suggest that the difference between pools depends on population heterogeneity. Principal coordinate analysis, hierarchical clustering, and the estimation of admixture coefficients derived froma single plant, in silico-pool, and seq-pool successfully identified the well-known structure of Andean and Mesoamerican gene pools of P. vulgaris across all datasets. Conclusion: In conclusion, seq-pool proved to be a viable approach for characterizing common bean germplasm compared to genotyping individual plants separately by balancing genotyping effort and costs. This study provides insights and serves as a valuable guide for gene bank researchers embarking on genotyping initiatives to characterize their collections. It aids curators in effectively managing the collections and facilitates marker-trait association studies, enabling the identification of candidate markers for key traits. [ABSTRACT FROM AUTHOR]

Objective: This study aimed to introduce a lytic bacteriophage against Staphylococcus saprophyticus from wastewater in Gorgan, northern Iran. Results: The vB_SsapS-46 phage was isolated from urban wastewater and formed round and clear plaques on bacterial culture. It was visualized by electron microscopy and had a large head (approximately 106 nm) and a long tail (approximately 150 nm), indicating that it belongs to the Siphoviridae family. The host range of vB_SsapS-46 was determined using a spot test on 35 S. saprophyticus clinical isolates, and it was able to lyse 12 of the 35 clinical isolates (34%). Finally, the relationship between phage sensitivity and adherence genes was assessed, revealing no significant correlation between phage sensitivity and the frequency of adherence genes. The vB_SsapS-46 phage can be used alone or in a mixture in future studies to control urinary tract infections caused by this bacterium, especially in the elimination of drug-resistant pathogens. [ABSTRACT FROM AUTHOR]

Background: A member of the Interleukin-1 superfamily of cytokines, interleukin-18 (IL-18) is essential to the etiology and progression of inflammatory bowel disease (IBD), a chronic inflammatory illness that affects the digestive system. This study investigated the possible association between two genetic variations, IL-18 rs187238 and IL-18 rs1946518, and IBD in Iraqi patients. Methods: We evaluated the association of two SNPs of the IL-18 gene at rs187238 and rs1946518 in 54 IBD patients with 19 Crohn's disease (CD), 35 ulcerative colitis (UC), and 46 healthy controls using PCR-RFLP and PCR-AS techniques for detecting IL-18 rs187238 and IL-18 rs1946518, respectively, by extracting genomic DNA from blood samples. Results: Our findings indicated no statistically significant variations between the IL-18 rs187238 genotypes and incidences of CD and UC (P = 0.189 and 0.59, respectively). However, the allele frequency showed a significant difference with CD (P = 0.049) but did not show a significant association with UC (P = 0.887). There was no significant association between the genotype and allele frequency of IL-18 rs1946518C/A and CD risk (P = 0.171 and 0.053, respectively). However, there was a significant association between the genotype and allele frequency of IL-18 rs1946518C/A and the risk of developing UC (P = 0.028 and 0.002, respectively). Conclusion: The study revealed statistically significant distinctions between the genetic and allelic frequencies of IL-18 rs1946518 and the probability of developing UC. Nonetheless, there were no significant distinctions between them and CD. According to the research, there were no major differences between IL-18 rs187238 and the two diseases. The frequency of the C allele is connected to CD. [ABSTRACT FROM AUTHOR]

Background: Yaks are a vital livestock in the Qinghai-Tibetan Plateau area for providing food products, maintaining sustainable ecosystems, and promoting cultural heritage. Because of uncontrolled mating, it is impossible to estimate inbreeding level of yak populations using the pedigree-based approaches. With the aims to accurately evaluate inbreeding level of two Chinese yak populations (Maiwa and Jiulong), we obtained genome-wide single nucleotide polymorphisms (SNPs) by DNA sequencing and calculated five SNP-by-SNP estimators ( , , , , and ), as well as two segment-based estimators of runs of homozygosity (ROH, ) and homozygous-by-descent (HBD, ). Functional implications were analyzed for the positional candidate genes located within the related genomic regions. Results: A total of 151,675 and 190,955 high-quality SNPs were obtained from 71 Maiwa and 30 Jiulong yaks, respectively. Jiulong had greater genetic diversity than Maiwa in terms of allele frequency and nucleotide diversity. The two populations could be genetically distinguished by principal component analysis, with the mean differentiation index (Fst) of 0.0054. The greater genomic inbreeding levels of Maiwa yaks were consistently supported by all five SNP-by-SNP estimators. Based on simple proportion of homozygous SNPs (), a lower inbreeding level was indicated by three successfully sequenced old leather samples that may represent historical Maiwa yaks about five generations ago. There were 3304 ROH detected among all samples, with mean and median length of 1.97 Mb and 1.0 Mb, respectively. A total of 94 HBD segments were found among all samples, whereas 92 of them belonged to the shortest class with the mean length of 10.9 Kb. Based on the estimates of and , however, there was no difference in inbreeding level between Maiwa and Jiulong yaks. Within the genomic regions with the significant Fst or enriched by ROH, we found several candidate genes and pathways that have been reported to be related to diverse production traits in farm animals. Conclusions: We successfully evaluated the genomic inbreeding level of two Chinese yak populations. Although different estimators resulted in inconsistent conclusions on their genomic inbreeding levels, our results may be helpful to implement the genetic conservation and utilization programs for the two yak populations. [ABSTRACT FROM AUTHOR]

Detection of specific gene mutations in cell-free DNA (cfDNA) serves as a valuable cancer biomarker and is increasingly being explored as an appealing alternative to tissue-based methods. However, the lack of available reference materials poses challenges in accurately evaluating the performance of different assays. In this study, we present the development of a comprehensive reference material panel for cfDNA detection, encompassing nine hotspot mutations in KRAS/BRAF/EGFR/PIK3CA at three variant allele frequencies (VAFs), ranging from 0.33 to 23.9%. To mimic cfDNA, these reference materials were generated by enzymatically digesting cell-line DNA into approximately 154-bp to 173-bp fragments using a laboratory-developed reaction system. The VAFs for each variation were precisely determined through validated digital PCR assays with high accuracy. Furthermore, the reliability and applicability of this panel were confirmed through two independent NGS assays, yielding concordant results. Collectively, our findings suggest that this novel reference material panel holds great potential for validation, evaluation, and quality control processes associated with liquid biopsy assays. [ABSTRACT FROM AUTHOR]

Introduction: Natural hybridization is common and plays a crucial role in driving biodiversity in nature. Despite its significance, the understanding of hybridization in ferns remains inadequate. Therefore, it is imperative to study fern hybridization to gain a more comprehensive understanding of fern biodiversity. Our study delves into the role of hybridization in shaping fern species, employing Microlepia matthewii as a case study to investigate its origins of hybridization. Methods: We performed double digest Genotyping-by-sequencing (dd-GBS) on M. matthewii and its potential parent species, identifying nuclear and chloroplast SNPs. Initially, nuclear SNPs were employed to construct the three cluster analysis: phylogenetic tree, principal component analysis, and population structure analysis. Subsequently, to confirm whether the observed genetic mixture pattern resulted from hybridization, we utilized two methods: ABBABABA statistical values in the D-suite program and gene frequency covariance in the Treemix software to detect gene flow. Finally, we employed chloroplast SNPs to construct a phylogenetic tree, tracing the maternal origin. Results and discussion: The analysis of the nuclear SNP cluster revealed that M. matthewii possesses a genetic composition that is a combination of M. hancei and M. calvescens. Furthermore, the analysis provided strong evidence of significant gene flow signatures from the parental species to the hybrid, as indicated by the two gene flow analyses. The samples of M. matthewii cluster separately with M. hancei or M. calvescens on the chloroplast systematic tree. However, the parentage ratio significantly differs from 1:1, suggesting that M. matthewii is a bidirectional and asymmetrical hybrid offspring of M. hancei and M. calvescens. [ABSTRACT FROM AUTHOR]

Objective: Fructose-1,6-bisphosphatase deficiency (FBP1D) is a rare inborn error due to mutations in the FBP1 gene. The genetic spectrum of FBP1D in China is unknown, also nonspecific manifestations confuse disease diagnosis. We systematically estimated the FBP1D prevalence in Chinese and explored genotype-phenotype association. Methods: We collected 101 FBP1 variants from our cohort and public resources, and manually curated pathogenicity of these variants. Ninety-seven pathogenic or likely pathogenic variants were used in our cohort to estimate Chinese FBP1D prevalence by three Methods: 1) carrier frequency, 2) permutation and combination, 3) Bayesian framework. Allele frequencies (AFs) of these variants in our cohort, China Metabolic Analytics Project (ChinaMAP) and gnomAD were compared to reveal the different hotspots in Chinese and other populations. Clinical and genetic information of 122 FBP1D patients from our cohort and published literature were collected to analyze the genotype-phenotypes association. Phenotypes of 68 hereditary fructose intolerance (HFI) patients from our previous study were used to compare the phenotypic differences between these two fructose metabolism diseases. Results: The estimated Chinese FBP1D prevalence was 1/1,310,034. In the Chinese population, c.490G>A and c.355G>A had significantly higher AFs than in the non-Finland European population, and c.841G>A had significantly lower AF value than in the South Asian population (all p values < 0.05). The genotypephenotype association analyses showed that patients carrying homozygous c.841G>A were more likely to present increased urinary glycerol, carrying two CNVs (especially homozygous exon1 deletion) were often with hepatic steatosis, carrying compound heterozygous variants were usually with lethargy, and carrying homozygous variants were usually with ketosis and hepatic steatosis (all p values < 0.05). By comparing to phenotypes of HFI patients, FBP1D patients were more likely to present hypoglycemia, metabolic acidosis, and seizures (all p-value < 0.05). Conclusion: The prevalence of FBP1D in the Chinese population is extremely low. Genetic sequencing could effectively help to diagnose FBP1D. [ABSTRACT FROM AUTHOR]

Background: Heritability partitioning approaches estimate the contribution of different functional classes, such as coding or regulatory variants, to the genetic variance. This information allows a better understanding of the genetic architecture of complex traits, including complex diseases, but can also help improve the accuracy of genomic selection in livestock species. However, methods have mainly been tested on human genomic data, whereas livestock populations have specific characteristics, such as high levels of relatedness, small effective population size or long-range levels of linkage disequilibrium. Results: Here, we used data from 14,762 cows, imputed at the whole-genome sequence level for 11,537,240 variants, to simulate traits in a typical livestock population and evaluate the accuracy of two state-of-the-art heritability partitioning methods, GREML and a Bayesian mixture model. In simulations where a single functional class had increased contribution to heritability, we observed that the estimators were unbiased but had low precision. When causal variants were enriched in variants with low (< 0.05) or high (> 0.20) minor allele frequency or low (below 1st quartile) or high (above 3rd quartile) linkage disequilibrium scores, it was necessary to partition the genetic variance into multiple classes defined on the basis of allele frequencies or LD scores to obtain unbiased results. When multiple functional classes had variable contributions to heritability, estimators showed higher levels of variation and confounding between certain categories was observed. In addition, estimators from small categories were particularly imprecise. However, the estimates and their ranking were still informative about the contribution of the classes. We also demonstrated that using methods that estimate the contribution of a single category at a time, a commonly used approach, results in an overestimation. Finally, we applied the methods to phenotypes for muscular development and height and estimated that, on average, variants in open chromatin regions had a higher contribution to the genetic variance (> 45%), while variants in coding regions had the strongest individual effects (> 25-fold enrichment on average). Conversely, variants in intergenic or intronic regions showed lower levels of enrichment (0.2 and 0.6-fold on average, respectively). Conclusions: Heritability partitioning approaches should be used cautiously in livestock populations, in particular for small categories. Two-component approaches that fit only one functional category at a time lead to biased estimators and should not be used. [ABSTRACT FROM AUTHOR]

Understanding the contribution of gene-environment interactions (GxE) to complex trait variation can provide insights into disease mechanisms, explain sources of heritability, and improve genetic risk prediction. While large biobanks with genetic and deep phenotypic data hold promise for obtaining novel insights into GxE, our understanding of GxE architecture in complex traits remains limited. We introduce a method to estimate the proportion of trait variance explained by GxE (GxE heritability) and additive genetic effects (additive heritability) across the genome and within specific genomic annotations. We show that our method is accurate in simulations and computationally efficient for biobank-scale datasets. We applied our method to common array SNPs (MAF ≥ 1 %), fifty quantitative traits, and four environmental variables (smoking, sex, age, and statin usage) in unrelated white British individuals in the UK Biobank. We found 68 trait-E pairs with significant genome-wide GxE heritability (p < 0.05 / 200) with a ratio of GxE to additive heritability of ≈ 6.8 % on average. Analyzing ≈ 8 million imputed SNPs (MAF ≥ 0.1 %), we documented an approximate 28 % increase in genome-wide GxE heritability compared to array SNPs. We partitioned GxE heritability across minor allele frequency (MAF) and local linkage disequilibrium (LD) values, revealing that, like additive allelic effects, GxE allelic effects tend to increase with decreasing MAF and LD. Analyzing GxE heritability near genes highly expressed in specific tissues, we find significant brain-specific enrichment for body mass index (BMI) and basal metabolic rate in the context of smoking and adipose-specific enrichment for waist-hip ratio (WHR) in the context of sex. Pazokitoroudi et al. introduced a scalable method to estimate trait variance explained by gene-environment interactions across the genome and within specific genomic annotations. Application of the method in the UK Biobank uncovered significant genome-wide GxE heritability and enrichment of GxE heritability within population and functional genomic annotations. [ABSTRACT FROM AUTHOR]

Blackleg, caused by Leptosphaeria maculans, is a major oilseed rape disease present in many traditional growing regions worldwide and is also becoming increasingly important in Tunisia. Cultivation of oilseed rape varieties harbouring major blackleg resistance genes (R‐genes) in an integrated approach represents the most effective means of controlling this fungal threat; knowledge of the pathogen population structure is a prerequisite for its efficient use. Leptosphaeria spp. isolates (145) were collected and characterized during 2018–2020 from oilseed rape stubble originating from seven northern Tunisian regions. As in most other growing areas, L. maculans and L. biglobosa occurred together on diseased plants. Most isolates (78.6%, n = 114) were L. maculans. For those, a combination of differential set phenotyping and PCR‐based avirulence gene analyses was used to perform race typing. Overall, 30 races were identified, carrying at least two and up to eight of the 11 avirulence genes tested. The most predominant race (14.9%) carried AvrLm7, AvrLepR1 and AvrLepR2 and was detected each year and in six out of the seven regions sampled, indicating the effectiveness of major R‐genes Rlm7 and LepR1 in Tunisia. Rlm2 was overcome by all Tunisian L. maculans isolates, but most were avirulent on Rlm6 and LepR2 (>50%), suggesting that using these R‐genes would provide protection in Tunisia. Both mating types were identified, with MAT1‐2 being twice as frequent as MAT1‐1. These results endorse the need for regular race surveillance of the local pathogen population at temporal and regional scales to develop appropriate blackleg management strategies. [ABSTRACT FROM AUTHOR]

Motivation Genetic variants present differential effects on humans according to various environmental exposures, the so-called "gene–environment interactions" (GxE). Many diseases can be diagnosed with multiple traits, such as obesity, diabetes, and dyslipidemia. I developed a multivariate scale test (MST) for detecting the GxE of a disease with several continuous traits. Given a significant MST result, I continued to search for which trait and which E enriched the GxE signals. Simulation studies were performed to compare MST with the univariate scale test (UST). Results MST can gain more power than UST because of (1) integrating more traits with GxE information and (2) the less harsh penalty on multiple testing. However, if only few traits account for GxE, MST may lose power due to aggregating non-informative traits into the test statistic. As an example, MST was applied to a discovery set of 93 708 Taiwan Biobank (TWB) individuals and a replication set of 25 200 TWB individuals. From among 2 570 487 SNPs with minor allele frequencies ≥5%, MST identified 18 independent variance quantitative trait loci (P  < 2.4E−9 in the discovery cohort and P  < 2.8E−5 in the replication cohort) and 41 GxE signals (P  < .00027) based on eight trait domains (including 29 traits). Availability and implementation https://github.com/WanYuLin/Multivariate-scale-test-MST- [ABSTRACT FROM AUTHOR]

INTRODUCTION: Systemic Lupus Erythematosus (SLE) is a systematic autoimmune disorder characterized by the production of autoantibodies against nuclear antigens and inflammation initiation. We aimed to examine the correlation between IL-18, IL-37, and PD-1, and the potential link between polymorphisms in the PD-1 gene located in intron-4 and the susceptibility to SLE. MATERIALS AND METHODS: This cross-sectional study included 43 SLE and 53 healthy individuals. Blood samples were obtained and underwent biochemical examination. The polymorphisms were screened by amplifying the intron-4 of the PD-1 gene using particular primers and then verified through sequencing. RESULTS: Our findings demonstrated statistically significant positive correlations between IL-18, IL-37, and PD-1, while the AUC of the ROC curve is 0.985, 0.968, and 0.940, and cut-off concentration is ≥132.87, ≥62.98, and ≥169.02, respectively. Moreover, two separate SNPs (rs6705653 and rs41386349) were discovered within intron-4 of the PD-1 gene. The genotype AA of the +7499 (G/A) SNP was significantly related with an increased risk of SLE (OR=3.11, 95%CI=1.52–5.94, p-value=0.031). Additionally, the A allele was identified as a risk allele (OR=1.59, 95%CI=1.09–2.31, p-value=0.043). Nevertheless, our study didn’t find any noteworthy connection between the allele and genotype of the +7209 (C/T) polymorphism region of the PD-1 gene frequencies and the susceptibility to SLE. CONCLUSION: IL-18, IL-37, and PD-1 may play significant roles in SLE immune responses and processes. Furthermore, the sequencing examination of intron-4 within the PD-1 gene demonstrated a noteworthy correlation between the A allele and the AA genotypes of PD-1 +7499 (G/A) SNP presence with the increased SLE susceptibility. [ABSTRACT FROM AUTHOR]

Genetic adaptation is the change of a population toward a phenotype that best fits the present ecological conditions of the environment it inhabits. As environmental conditions change, allele frequencies shift, resulting in different populations of the same species possessing genetic variation and divergent phenotypes. Cooperatively breeding common mole-rats (Cryptomys hottentotus hottentotus) inhabit environments along an aridity gradient in South Africa, which provides an opportunity for local genetic adaptations to occur. Using one mitochondrial gene (cytochrome b) and 3,540 SNP loci across the whole genome, we determined the phylogenetic relationship, population structure and genetic diversity of five populations of C. h. hottentotus located along an aridity gradient. Mitochondrial data identified population-specific clades that were less distinct in the two mesic populations, potentially indicating historical or recent gene flow, or the retention of ancestral haplotypes. Arid and semi-arid populations formed a distinct cluster from the non-arid populations. Genetic diversity and gene flow were higher in arid-dwelling individuals, suggesting greater connectivity and interactions between colonies in arid regions in comparison to mesic ones. Using an Aridity Index, we determined that isolation by environment, rather than isolation by geographical distance, best explains the genetic distance between the populations. Further analyses using target loci may determine if there are differing underlying genetic adaptations among populations of C. h. hottentotus. These analyses could help unravel population differences in response to environmental factors within a subspecies of bathyergid mole-rat and determine the adaptive capacity of this small nonmigratory subterranean rodent species in response to aridification in the face of climate change. [ABSTRACT FROM AUTHOR]

Mangrove sediments host a diverse array of microbial populations and are characterized by high heterogeneity along their vertical depths. However, the genetic diversity within these populations is largely unknown, hindering our understanding of their adaptive evolution across the sediment depths. To elucidate their genetic diversity, we utilized metagenome sequencing to identify 16 high‐frequency microbial populations comprised of two archaea and 14 bacteria from mangrove sediment cores (0–100 cm, with 10 depths) in Qi'ao Island, China. Our analysis of the genome‐wide genetic variation revealed extensive nucleotide diversity in the microbial populations. The genes involved in the transport and the energy metabolism displayed a high nucleotide diversity (HND; 0.0045–0.0195; an indicator of shared minor alleles with the microbial populations). By tracking the processes of homologous recombination, we found that each microbial population was subjected to different purification selection levels at different depths (44.12% genes). This selection resulted in significant differences in synonymous/non‐synonymous mutation ratio between 0–20 and 20–100 cm layers, indicating the adaptive evolutionary process of microbial populations. Furthermore, our assessment of differentiation in the allele frequencies between these two layers showed that the functional genes involved in the metabolic processes of amino acids or cofactors were highly differential in more than half of them. Together, we showed that the nucleotide diversity of microbial populations was shaped by homologous recombination and gene‐specific selection, finally resulting in the stratified differentiation occurring between 0–20 and 20–100 cm. These results enhance our cognition of the microbial adaptation mechanisms to vertical environmental changes during the sedimentation process of coastal blue carbon ecosystems. [ABSTRACT FROM AUTHOR]

Crop wild relatives (CWR) provide a valuable resource for improving crops. They possess desirable traits that confer resilience to various environmental stresses. To fully utilize crop wild relatives in breeding and conservation programs, it is important to understand the genetic basis of their adaptation. Landscape genomics associates environments with genomic variation and allows for examining the genetic basis of adaptation. Our study examined the differences in allele frequency of 15,416 single nucleotide polymorphisms (SNPs) generated through genotyping by sequencing approach among 153 accessions of 15 wild eggplant relatives and two cultivated species from Africa, the principal hotspot of these wild relatives. We also explored the correlation between these variations and the bioclimatic and soil conditions at their collection sites, providing a comprehensive understanding of the genetic signals of environmental adaptation in African wild eggplant. Redundancy analysis (RDA) results showed that the environmental variation explained 6% while the geographical distances among the collection sites explained 15% of the genomic variation in the eggplant wild relative populations when controlling for population structure. Our findings indicate that even though environmental factors are not the main driver of selection in eggplant wild relatives, it is influential in shaping the genomic variation over time. The selected environmental variables and candidate SNPs effectively revealed grouping patterns according to the environmental characteristics of sampling sites. Using four genotype–environment association methods, we detected 396 candidate SNPs (2.5% of the initial SNPs) associated with eight environmental factors. Some of these SNPs signal genes involved in pathways that help adapt to environmental stresses such as drought, heat, cold, salinity, pests, and diseases. These candidate SNPs will be useful for marker‐assisted improvement and characterizing the germplasm of this crop for developing climate‐resilient eggplant varieties. The study provides a model for applying landscape genomics to other crops' wild relatives. [ABSTRACT FROM AUTHOR]

The Merovingian period (5th to 8th cc AD) was a time of demographic, socioeconomic, cultural, and political realignment in Western Europe. Here, we report the whole-genome shotgun sequence data of 30 human skeletal remains from a coastal Late Merovingian site of Koksijde (675 to 750 AD), alongside 18 remains from two Early to Late Medieval sites in present-day Flanders, Belgium. We find two distinct ancestries, one shared with Early Medieval England and the Netherlands, while the other, minor component, reflecting likely continental Gaulish ancestry. Kinship analyses identified no large pedigrees characteristic to elite burials revealing instead a high modularity of distant relationships among individuals of the main ancestry group. In contrast, individuals with >90% Gaulish ancestry had no kinship links among sampled individuals. Evidence for population structure and major differences in the extent of Gaulish ancestry in the main group, including in a mother-daughter pair, suggests ongoing admixture in the community at the time of their burial. The isotopic and genetic evidence combined supports a model by which the burials, representing an established coastal nonelite community, had incorporated migrants from inland populations. The main group of burials at Koksijde shows an abundance of >5 cM long shared allelic intervals with the High Medieval site nearby, implying long-term continuity and suggesting that similarly to Britain, the Early Medieval ancestry shifts left a significant and long-lasting impact on the genetic makeup of the Flemish population. We find substantial allele frequency differences between the two ancestry groups in pigmentation and diet-associated variants, including those linked with lactase persistence, likely reflecting ancestry change rather than local adaptation. [ABSTRACT FROM AUTHOR]

Aim: Short tandem repeats (STRs) are repetitive DNA sequences and highly mutable in various human disorders. While the involvement of STRs in various genetic disorders has been extensively studied, their role in autism spectrum disorder (ASD) remains largely unexplored. In this study, we aimed to investigate genetic association of STR expansions with ASD using whole genome sequencing (WGS) and identify risk loci associated with ASD phenotypes. Methods: We analyzed WGS data of 634 ASD families and performed genome‐wide evaluation for 12,929 STR loci. We found rare STR expansions that exceeded normal repeat lengths in autism cases compared to unaffected controls. By integrating single cell RNA and ATAC sequencing datasets of human postmortem brains, we prioritized STR loci in genes specifically expressed in cortical development stages. A deep learning method was used to predict functionality of ASD‐associated STR loci. Results: In ASD cases, rare STR expansions predominantly occurred in early cortical layer‐specific genes involved in neurodevelopment, highlighting the cellular specificity of STR‐associated genes in ASD risk. Leveraging deep learning prediction models, we demonstrated that these STR expansions disrupted the regulatory activity of enhancers and promoters, suggesting a potential mechanism through which they contribute to ASD pathogenesis. We found that individuals with ASD‐associated STR expansions exhibited more severe ASD phenotypes and diminished adaptability compared to non‐carriers. Conclusion: Short tandem repeat expansions in cortical layer‐specific genes are associated with ASD and could potentially be a risk genetic factor for ASD. Our study is the first to show evidence of STR expansion associated with ASD in an under‐investigated population. [ABSTRACT FROM AUTHOR]

Oleaceae, a eudicot family with great species diversity, has attracted much attention from botanists because it contains many plants with important economic, medicinal, and ornamental values. However, the history of polyploidization and ancestral genome reshuffling of Oleaceae remains unclear. Here, we clarified an Oleaceae‐common hexaploidization (OCH) event occurring at ~53–61 million years ago (Ma) common in all Oleaceae plants and an Oleaceae‐recent tetraploidization (ORT) event occurring at ~18–21 Ma shared by the lineages of Syringa, Olea, Osmanthus, and Fraxinus. We found that high‐frequency polyploidization events drove the frequency of gene loss in Oleaceae genomes and extended the size of regions containing adjacent gene loss, thereby promoting the degree of genome fragmentation. We revealed that biased fractionation between the OCH‐ and ORT‐produced subgenomes is likely attributed to the origin of allopolyploidization in the OCH and ORT events. Significantly, through paleochromosome rearrangement comparisons, we proposed a "two‐step" genome duplication model for OCH and determined the duplicated orders of OCH tripled genome. We reconstructed 11 protochromosomes of the most recent ancestral Oleaceae karyotype (AOK) and elucidated the trajectories of immense paleochromosome reorganization of Oleaceae species from ancestral eudicot karyotype. Notably, we tracked the diversification history of secondary metabolite synthesis genes in the Oleaceae and explored the effects of paleogenome evolution on specialized metabolite synthesis. Our findings provide new insights into the polyploidization and paleogenomic evolution of Oleaceae and have important scientific significance for understanding the genetic basis of species and secondary metabolic diversity in Oleaceae. [ABSTRACT FROM AUTHOR]

In dairy cattle production, it is important to understand how inbreeding affects production, fertility, and health traits. However, there is still limited use of genomic information to estimate inbreeding, despite advancements in genotyping technologies. To address this gap, we investigated the effect of inbreeding on German Holstein dairy cattle using both pedigree-based and genomic-based inbreeding estimators. We employed one method based on pedigree information (F ped) together with 6 genomic-based methods, including 3 genome-wide complex trait analysis software estimators (F hat1 , F hat2 , F hat3), VanRaden's first method (F VR1 , with observed allele frequencies, and F VR0.5 , when allele frequencies are set to 0.5), and one based on runs of homozygosity (F roh). Data from 24,489 cows with both phenotypes and genotypes were used, with a pedigree including 232,780 animals born between 1970 and 2018. We analyzed the effects of inbreeding depression on production, fertility, and health traits separately, using single-trait linear animal models as well as threshold models to account for the binary nature of the health traits. For the health traits, we transformed solutions from the liability scale to a probability scale for easier interpretation. Our results showed that the mean inbreeding coefficients from all estimators ranged from −0.003 to 0.243, with negative values observed for most genomic-based methods. We found out that a 1% increase in inbreeding caused a depression ranging from 25.94 kg (F hat1) to 40.62 kg (F hat3), 1.18 kg (F hat2) to 1.70 kg (F hat3), 0.90 kg (F hat2) to 1.45 kg (F roh and F hat3), 0.19 (F ped) to 0.34 d (F hat3) for 305-d milk yield, fat, protein, and calving interval, respectively. The health traits showed very slight gradual changes when inbreeding was increased steadily from 0% to 50%, with digital dermatitis showing a rather contrasting trend to that of mastitis, which increased the more an animal was inbred. Overall, our study highlights the importance of considering both pedigree-based and genomic-based inbreeding estimators when assessing the impact on inbreeding, emphasizing that not all inbreeding is harmful. [ABSTRACT FROM AUTHOR]

Interest in recovering DNA from the surface of ammunition evidence for genotyping has increased over the past few years. Numerous studies have examined a variety of methods to maximize DNA recovery from these types of challenging samples, but successful DNA profiling has been inconsistent. Low amounts of DNA and PCR inhibition due to metal ions have been suggested as the leading causes of poor results; however, no study quantitatively examined the presence of metal ions at various stages of the DNA analysis workflow from DNA collection through to amplification. In this study, the effectiveness of six different DNA collection and purification methods commonly used by forensic laboratories to process brass ammunition for DNA evidence was investigated. The amount of copper, zinc, and other metals co-recovered from fired and unfired brass casings during DNA collection (using numerous soaking, swabbing, and direct PCR protocols) was quantified via Inductively Coupled Plasma - Optical Emission Spectrometry (ICP-OES). This same panel of metals was subsequently quantified after DNA lysis and purification steps. Results demonstrated that low amounts of DNA, DNA damage, and degradation are more detrimental to STR typing results than PCR inhibition, as metal ions were successfully removed by all DNA purification methods tested. Furthermore, the use of metal ion chelators increased the amount of DNA recovered and number of reportable STR alleles. This research informs the forensic community on the most effective way to collect and process trace amounts of biological material from brass ammunition and similar evidence. [ABSTRACT FROM AUTHOR]

Genome-wide association studies (GWAS) are often underpowered due to small effect sizes of common single nucleotide polymorphisms (SNPs) on phenotypes and extreme multiple testing thresholds. The most common approach for increasing statistical power is to increase sample size. We propose an alternative strategy of redefining case-control outcomes into ordinal case-subthreshold-asymptomatic variables. While maintaining the clinical case threshold, we subdivide controls into two groups: individuals who are symptomatic but do not meet the clinical criteria for diagnosis (subthreshold) and individuals who are effectively asymptomatic. We conducted a simulation study to examine the impact of effect size, minor allele frequency, population prevalence, and the prevalence of the subthreshold group on statistical power to detect genetic associations in three scenarios: a standard case-control, an ordinal, and a case-asymptomatic control analysis. Our results suggest the ordinal model consistently provides the greatest statistical power while the case-control model the least. Power in the case-asymptomatic control model reflects the case-control or ordinal model depending on the population prevalence and size of the subthreshold category. We then analyzed a major depression phenotype from the UK Biobank to corroborate our simulation results. Overall, the ordinal model improves statistical power in GWAS consistent with increasing the sample size by approximately 10%. [ABSTRACT FROM AUTHOR]

The Mediterranean fever (MEFV) gene codes for protein pyrin, which is among the modulators of inflammasome activity in innate immune cells. It was suggested that there is a relation between MEFV variations and inflammatory diseases. The aim of this study was to investigate MEFV gene variations in the patients with primary dysmenorrhea. The prevalence of common MEFV gene variations (M694V, M680I, V726A, E148Q and R202Q) was investigated in 145 young women with primary dysmenorrhea and 135 unrelated healthy controls. MEFV gene variations were genotyped using PCR-based RFLP assay. Number of childbirth and marriage were significantly lower in the study group than the controls, respectvely (p < 0.001, p = 0.001). Family history was statistically higher in the patient group (p < 0.001). In total, MEFV genotype and allele frequencies were significantly higher in patients than controls, respectively (p = 0.008 and p = 0.005, respectively). It was found that MEFV gene E148Q allele was more common in patient group (p = 0.039). MEFV R202Q A allele was higher in the patients than the controls (p = 0.045). A significant association was observed when the patients were compared with the controls according to R202Q variant AA versus GG+GA genotypes (p=0.020). Our findings suggest that MEFV variations may be a risk factor for patients with dysmenorrhea in a Turkish cohort. There are very few studies in the literature regarding the relationship between pathological variants of MEFV and dysmenorrhea disease. The common MEFV mutations/variants were evaluated in primary dysmenorrhea patients. Family history was statistically higher in the patient group (p <.001). MEFV gene variations were found 52 (35.9%) in patients and 29 (21.5%2) in controls. MEFV gene allele frequency was significantly higher in-patient group than control (p =.005). [ABSTRACT FROM AUTHOR]

There is growing interest in uncovering genetic kinship patterns in past societies using low‐coverage palaeogenomes. Here, we benchmark four tools for kinship estimation with such data: lcMLkin, NgsRelate, KIN, and READ, which differ in their input, IBD estimation methods, and statistical approaches. We used pedigree and ancient genome sequence simulations to evaluate these tools when only a limited number (1 to 50 K, with minor allele frequency ≥0.01) of shared SNPs are available. The performance of all four tools was comparable using ≥20 K SNPs. We found that first‐degree related pairs can be accurately classified even with 1 K SNPs, with 85% F1 scores using READ and 96% using NgsRelate or lcMLkin. Distinguishing third‐degree relatives from unrelated pairs or second‐degree relatives was also possible with high accuracy (F1 > 90%) with 5 K SNPs using NgsRelate and lcMLkin, while READ and KIN showed lower success (69 and 79% respectively). Meanwhile, noise in population allele frequencies and inbreeding (first‐cousin mating) led to deviations in kinship coefficients, with different sensitivities across tools. We conclude that using multiple tools in parallel might be an effective approach to achieve robust estimates on ultra‐low‐coverage genomes. [ABSTRACT FROM AUTHOR]

Two important characteristics of metapopulations are extinction–(re)colonization dynamics and gene flow between subpopulations. These processes can cause strong shifts in genome-wide allele frequencies that are generally not observed in "classical" (large, stable, and panmictic) populations. Subpopulations founded by one or a few individuals, the so-called propagule model, are initially expected to show intermediate allele frequencies at polymorphic sites until natural selection and genetic drift drive allele frequencies toward a mutation–selection–drift equilibrium characterized by a negative exponential-like distribution of the site frequency spectrum. We followed changes in site frequency spectrum distribution in a natural metapopulation of the cyclically parthenogenetic pond-dwelling microcrustacean Daphnia magna using biannual pool-seq samples collected over a 5-yr period from 118 ponds occupied by subpopulations of known age. As expected under the propagule model, site frequency spectra in newly founded subpopulations trended toward intermediate allele frequencies and shifted toward right-skewed distributions as the populations aged. Immigration and subsequent hybrid vigor altered this dynamic. We show that the analysis of site frequency spectrum dynamics is a powerful approach to understand evolution in metapopulations. It allowed us to disentangle evolutionary processes occurring in a natural metapopulation, where many subpopulations evolve in parallel. Thereby, stochastic processes like founder and immigration events lead to a pattern of subpopulation divergence, while genetic drift leads to converging site frequency spectrum distributions in the persisting subpopulations. The observed processes are well explained by the propagule model and highlight that metapopulations evolve differently from classical populations. [ABSTRACT FROM AUTHOR]

Simple Summary: The BMPR1B and T/Brachyury genes are considered major genetic factors influencing sheep reproduction and tail bone number. In this study, we conducted an association analysis between nine mutations of the BMPR1B gene and the litter sizes of short-tailed Gobi sheep and larger-tailed Ujimqin sheep. We also compared the frequency of the favorable alleles of four mutations among Gobi short tail sheep, Ujimqin sheep, and Mongolia sheep. In addition, we confirmed that the frequency of tail-bone-number-related T alleles at position c.363G>T in the T gene was significantly higher in short-tailed Gobi sheep than larger-tailed Ujimqin sheep. The significant deposition of tail fat in sheep has a profound impact on the economic benefits of animal husbandry. Furthermore, increasing the litter size is a crucial means of enhancing economic benefits. The BMPR1B and T/Brachyury genes are considered major functional genes that could affect sheep litter size and tail bone number, respectively. In this study, we employed direct sequencing to identify specific mutations of the BMPR1B gene in Gobi short tail sheep and carried out genotyping using MassARRAY technology for each variant of both the BMPR1B and T genes. Significant associations were demonstrated between the c.687G>A mutation of BMPR1B and the litter size in both the Gobi short tail sheep and Ujimqin sheep breeds. Meanwhile, the g.30058882_30058873GCAGATTAAAIndel mutation was significantly associated with the litter size in Gobi short tail sheep. These findings may provide valuable genetic markers for expanding sheep litter size. In addition, we also confirmed that the frequency of tail-bone-number-related T alleles was significantly higher in Gobi short tail sheep than in longer-tailed Ujimqin sheep. [ABSTRACT FROM AUTHOR]

Objective Methods Results Conclusion Our objective was to evaluate whether there is an enrichment of rare variants in familial hemophagocytic lymphohistiocytosis (HLH)‐associated genes among patients with systemic juvenile idiopathic arthritis (sJIA) with or without macrophage activation syndrome (MAS).Targeted sequencing of HLH genes (LYST, PRF1, RAB27A, STX11, STXBP2, UNC13D) was performed in patients with sJIA from an established cohort. Sequence data from control participants were obtained in silico (database of Genotypes and Phenotypes: phs000280.v8.p2). Rare variant association testing (RVT) was performed with sequence kernel association test package. Significance was defined as P < 0.05 after 100,000 permutations.Sequencing data from 524 sJIA cases were jointly called and harmonized with exome‐derived target data from 3,000 controls. Quality control operations produced a set of 480 cases and 2,924 ancestrally matched control participants. RVT of cases and controls revealed a significant association with rare protein‐altering variants (minor allele frequency [MAF] < 0.01) of STXBP2 (P = 0.020) and ultrarare variants (MAF < 0.001) of STXBP2 (P = 0.006) and UNC13D (P = 0.046). A subanalysis of 32 cases with known MAS and 90 without revealed a significant difference in the distribution of rare UNC13D variants (P = 0.0047) between the groups. Additionally, patients with sJIA more often carried two or more HLH variants than did controls (P = 0.007), driven largely by digenic combinations involving LYST.We identified an enrichment of rare HLH variants in patients with sJIA compared with controls, driven by STXBP2 and UNC13D. Biallelic variation in HLH genes was associated with sJIA, driven by LYST. Only UNC13D displayed enrichment in patients with MAS. This suggests that HLH variants may contribute to the pathophysiology of sJIA, even without MAS. [ABSTRACT FROM AUTHOR]

Background: High blood pressure, also known as hypertension (HTN), is a complicated disorder that is controlled by a complex network of physiological processes. Untreated hypertension is associated with increased death incidence, rise the need for understanding the genetic basis affecting hypertension susceptibility and development. The current study sought to identify the genetic association between twelve single nucleotide polymorphisms (SNPs) within seven candidate genes (NOS3, NOS1AP, REN, PLA2G4A, TCF7L, ADRB1, and PTPRD) Methods: The current study included 200 Jordanian individuals diagnosed with hypertension, compared to 224 healthy controls. Whole blood samples were drawn from each individual for DNA isolation and genotyping. The SNPStats tool was used to assess haplotype, genotype, and allele frequencies by the mean of chi-square (χ2). Results: Except for rs10739150 of PTPRD (P = 0.0003), the genotypic and allelic distribution of the SNP was identical between patients and controls. The prevalence of the G/G genotype in healthy controls (45.5%) was lower than in hypertension patients (64.3%), suggesting that it might be a risk factor for the disease. PTPRD TTC genetic haplotypes were strongly linked with hypertension (P = 0.003, OR = 4.03). Conclusion: This study provides a comprehensive understanding of the involvement of rs10739150 within the PTPRD gene in hypertension. This new knowledge could potentially transform the way we approach hypertension diagnosis, providing an accurate diagnostic tool for classifying individuals who are at a higher risk of developing this condition. [ABSTRACT FROM AUTHOR]

Introduction: Identifying new biomarkers for predicting immune checkpoint inhibitors (ICIs) response in non-small cell lung cancer (NSCLC) is crucial. We aimed to assess the variant allele frequency (VAF)-related profile as a novel biomarker for NSCLC personalized therapy. Methods: We utilized genomic data of 915 NSCLC patients via cBioPortal and a local cohort of 23 patients for model construction and mutational analysis. Genomic, transcriptomic data from 952 TCGA NSCLC patients, and immunofluorescence (IF) assessment with the local cohort supported mechanism analysis. Results: Utilizing the random forest algorithm, a 15-gene VAF-related model was established, differentiating patients with durable clinical benefit (DCB) from no durable benefit (NDB). The model demonstrated robust performance, with ROC-AUC values of 0.905, 0.737, and 0.711 across training (n = 313), internal validation (n = 133), and external validation (n = 157) cohorts. Stratification by the model into high- and low-score groups correlated significantly with both progression-free survival (PFS) (training: P < 0.0001, internal validation: P < 0.0001, external validation: P = 0.0066) and overall survival (OS) (n = 341) (P < 0.0001). Notably, the stratification system was independent of PD-L1 (P < 0.0001) and TMB (P < 0.0001). High-score patients exhibited an increased DCB ratio and longer PFS across both PD-L1 and TMB subgroups. Additionally, the high-score group appeared influenced by tobacco exposure, with activated DNA damage response pathways. Whereas, immune/inflammation-related pathways were enriched in the low-score group. Tumor immune microenvironment analyses revealed higher proportions of exhausted/effector memory CD8 + T cells in the high-score group. Conclusions: The mutational VAF profile is a promising biomarker for ICI therapy in NSCLC, with enhanced therapeutic stratification and management as a supplement to PD-L1 or TMB. [ABSTRACT FROM AUTHOR]

Central African rainforests are predicted to be disproportionately affected by future climate change. How species will cope with these changes is unclear, but rapid environmental changes will likely impose strong selection pressures. Here we examined environmental drivers of genomic variation in the central African puddle frog (Phrynobatrachus auritus) to identify areas of elevated environmentally-associated turnover. We also compared current and future climate models to pinpoint areas of high genomic vulnerability where allele frequencies will have to shift the most in order to keep pace with future climate change. Neither physical landscape barriers nor the effects of past Pleistocene refugia influenced genomic differentiation. Alternatively, geographic distance and seasonal aspects of precipitation are the most important drivers of SNP allele frequency variation. Patterns of genomic differentiation coincided with key ecological gradients across the forest-savanna ecotone, montane areas, and a coastal to interior rainfall gradient. Areas of greatest vulnerability were found in the lower Sanaga basin, the southeastern region of Cameroon, and southwest Gabon. In contrast with past conservation efforts that have focused on hotspots of species richness or endemism, our findings highlight the importance of maintaining environmentally heterogeneous landscapes to preserve genomic variation and ongoing evolutionary processes in the face of climate change. [ABSTRACT FROM AUTHOR]