Three de novo assembled wild cacao genomes from the Upper Amazon

Theobroma cacao, the chocolate tree, is indigenous to the Amazon basin, the greatest biodiversity hotspot on earth. Recent advancement in plant genomics highlights the importance of de novo sequencing of multiple reference genomes to capture the genome diversity present in different cacao populations. In this study, three high-quality chromosome-level genomes of wild cacao were constructed, de novo assembled with HiFi long reads sequencing, and scaffolded using a reference-free strategy. These genomes represent the three most important genetic clusters of cacao trees from the Upper Amazon region. The three wild cacao genomes were compared with two reference genomes of domesticated cacao. The five cacao genetic clusters were inferred to have diverged in the early and middle Pleistocene period, approximately 1.83-0.69 million years ago. The results shown here serve as an example of understanding how the Amazonian biodiversity was developed. The three wild cacao genomes provide valuable resources for studying genetic diversity and advancing genetic improvement of this species.

Reduced evolutionary potential under heat and drought stress at the southern range edge of North American Arabidopsis lyrata

The warm edges of species' distributions are vulnerable under global warming. Evidence is the recent range retraction from there found in many species. It is unclear why populations cannot easily adapt to warmer, drier, or combined hot and dry conditions and locally persist. Here, we assessed the ability to adapt to these stressors in the temperate species Arabidopsis lyrata. We grew plants from replicate seed families of a central population with high genetic diversity under a temperature and precipitation regime typical of the low-latitude margin or under hotter and/or drier conditions within naturally occurring amplitudes. We then calculated genetic variance-covariance (G-) matrices of traits depicting growth and allocation as well as selection vectors to compare the predicted adaptation potential under the different climate-stress regimes. We found that the sum of genetic variances and genetic correlations were not significantly different under stress as compared to benign conditions. But under drought and heat-drought, the predicted ability to adapt was severely constrained due to strong selection and selection pointing in a direction with less multivariate genetic variation. The much-reduced ability to adapt to dry and hot-dry conditions is likely to reduce the persistence of populations at the low-latitude margin of the species' distribution and contributes to the local extinction of the species under further warming.

Allopolyploid origin and diversification of the Hawaiian endemic mints

Island systems provide important contexts for studying processes underlying lineage migration, species diversification, and organismal extinction. The Hawaiian endemic mints (Lamiaceae family) are the second largest plant radiation on the isolated Hawaiian Islands. We generated a chromosome-scale reference genome for one Hawaiian species, Stenogyne calaminthoides, and resequenced 45 relatives, representing 34 species, to uncover the continental origins of this group and their subsequent diversification. We further resequenced 109 individuals of two Stenogyne species, and their purported hybrids, found high on the Mauna Kea volcano on the island of Hawai'i. The three distinct Hawaiian genera, Haplostachys, Phyllostegia, and Stenogyne, are nested inside a fourth genus, Stachys. We uncovered four independent polyploidy events within Stachys, including one allopolyploidy event underlying the Hawaiian mints and their direct western North American ancestors. While the Hawaiian taxa may have principally diversified by parapatry and drift in small and fragmented populations, localized admixture may have played an important role early in lineage diversification. Our genomic analyses provide a view into how organisms may have radiated on isolated island chains, settings that provided one of the principal natural laboratories for Darwin's thinking about the evolutionary process.

Genetic Variation and Composition of Two Commercial Estonian Dairy Cattle Breeds Assessed by SNP Data

The aims of this study were to assess the genomic relatedness of Estonian and selected European dairy cattle breeds and to examine the within-breed diversity of two Estonian dairy breeds using genome-wide SNP data. This study was based on a genotyped heifer population of the Estonian Red (ER) and Estonian Holstein (EH) breeds, including about 10% of all female cattle born in 2017-2020 (sample sizes n = 215 and n = 2265, respectively). The within-breed variation study focused on the level of inbreeding using the ROH-based inbreeding coefficient. The genomic relatedness analyses were carried out among two Estonian and nine European breeds from the WIDDE database. Admixture analysis revealed the heterogeneity of ER cattle with a mixed pattern showing several ancestral populations containing a relatively low proportion (1.5-37.0%) of each of the reference populations used. There was a higher FROH in EH (FROH = 0.115) than in ER (FROH = 0.044). Compared to ER, the long ROHs of EH indicated more closely related parents. The paternal origin of the genetic material used in breeding had a low effect on the inbreeding level. However, among EH, the highest genomic inbreeding was estimated in daughters of USA-born sires.

The drivers of intraspecific trait variation and their implications for future tree productivity and survival

Forests are facing unprecedented levels of stress from pest and disease outbreaks, disturbance, fragmentation, development, and a changing climate. These selective agents act to alter forest composition from regional to cellular levels. Thus, a central challenge for understanding how forests will be impacted by future change is how to integrate across scales of biology. Phenotype, or an observable trait, is the product of an individual's genes (G) and the environment in which an organism lives (E). To date, researchers have detailed how environment drives variation in tree phenotypes over long time periods (e.g., long-term ecological research sites [LTERs]) and across large spatial scales (e.g., flux network). In parallel, researchers have discovered the genes and pathways that govern phenotypes, finding high degrees of genetic control and signatures of local adaptation in many plant traits. However, the research in these two areas remain largely independent of each other, hindering our ability to generate accurate predictions of plant response to environment, an increasingly urgent need given threats to forest systems. I present the importance of both genes and environment in determining tree responses to climate stress. I highlight why the difference between G versus E in driving variation is critical for our understanding of climate responses, then propose means of accelerating research that examines G and E simultaneously by leveraging existing long-term, large-scale phenotypic data sets from ecological networks and adding newly affordable sequence (-omics) data to both drill down to find the genes and alleles influencing phenotypes and scale up to find how patterns of demography and local adaptation may influence future response to change.

Association of the neonatal Fc receptor promoter variable number of tandem repeat polymorphism with immunoglobulin response in patients with chronic inflammatory demyelinating polyneuropathy

BACKGROUND AND PURPOSE

Chronic inflammatory demyelinating polyneuropathy (CIDP) is an autoimmune disease with humoral and cellular autoimmunity causing demyelination of peripheral nerves, commonly treated with intravenous immunoglobulins (IVIg). The neonatal Fc receptor (FcRn), encoded by the FCGRT gene, prevents the degradation of immunoglobulin G (IgG) by recycling circulating IgG. A variable number of tandem repeat (VNTR) polymorphism in the promoter region of the FCGRT gene is associated with different expression levels of mRNA and protein. Thus, patients with genotypes associated with relatively low FcRn expression may show a poorer treatment response to IVIg due to increased IVIg degradation.

METHODS

VNTR genotypes were analyzed in 144 patients with CIDP. Patients' clinical data, including neurological scores and treatment data, were collected as part of the Immune-Mediated Neuropathies Biobank registry.

RESULTS

Most patients (n = 124, 86%) were VNTR 3/3 homozygotes, and 20 patients (14%) were VNTR 2/3 heterozygotes. Both VNTR 3/3 and VNTR 2/3 genotype groups showed no difference in clinical disability and immunoglobulin dosage. However, patients with a VNTR 2 allele were more likely to receive subcutaneous immunoglobulins (SCIg) than patients homozygous for the VNTR 3 allele (25% vs. 9.7%, p = 0.02) and were more likely to receive second-line therapy (75% vs. 54%, p = 0.05).

CONCLUSIONS

The VNTR 2/3 genotype is associated with the administration of SCIg, possibly reflecting a greater benefit from SCIg due to more constant immunoglobulin levels without lower IVIg levels between the treatment circles. Also, the greater need for second-line treatment in VNTR 2/3 patients could be an indirect sign of a lower response to immunoglobulins.

Causal role of immune cells in chronic obstructive pulmonary disease: Mendelian randomization study

OBJECTIVES

Innate and adaptive immunity play different roles in the pathogenesis of chronic obstructive pulmonary disease (COPD). However, previous studies on the relationship between immune cells and COPD reported inconsistent results.

METHODS

The causal connection between 731 immune cells and COPD was established using a two-sample Mendelian randomization (MR) analysis through publicly accessible genetic data. The heterogeneity and horizontal pleiotropism of the findings were confirmed using sensitivity analysis.

RESULTS

In the B-cell panel, B-cell activating factor receptor (BAFF-R) on CD20- and CD20 on IgD-CD38bright (OR (95% CI): 0.93 (0.88, 0.99) and 0.97 (0.95, 0.98), respectively) were discovered to be protective. In the cDC panel, CD62L- plasmacytoid DC AC, CD80 on monocytes and CD11c on myeloid DCs (OR (95% CI): 0.94 (0.92, 0.97), 0.97 (0.94, 0.99) and (0.97 (0.95, 0.98), respectively) exerted protective effects. However, unswitched memory AC (OR (95%CI): 1.08 (1.01,1.15)) and CD 19 on IgD- CD 27- (OR (95%CI): 1.06 (1.02,1.10)) were hazardous in the B-cell panel. However, among the 731 immune cell phenotypes, no causal relationship was found for COPD on immune cells.

CONCLUSION

This study found a potential causal relationship between immune cells in COPD, ruling out reverse causation. This study provides new avenues for studying the mechanisms of COPD.

Genetic variation drives differences in obesity-related gene regulation

Heikkinen and colleagues recently demonstrated that genetic variation, rather than dietary changes, governs gene regulation in liver. This finding highlights the impact of noncoding variants on chromatin accessibility, histone modifications, transcription factor binding, and gene expression and has implications for future research directions in understanding the genetic basis of disease.

Ecological characteristics explain neutral genetic variation of three coastal sparrow species

Eco-phylogeographic approaches to comparative population genetic analyses allow for the inclusion of intrinsic influences as drivers of intraspecific genetic structure. This insight into microevolutionary processes, including changes within a species or lineage, provides better mechanistic understanding of species-specific interactions and enables predictions of evolutionary responses to environmental change. In this study, we used single nucleotide polymorphisms (SNPs) identified from reduced representation sequencing to compare neutral population structure, isolation by distance (IBD), genetic diversity and effective population size (Ne) across three closely related and co-distributed saltmarsh sparrow species differing along a specialization gradient-Nelson's (Ammospiza nelsoni subvirgata), saltmarsh (A. caudacuta) and seaside sparrows (A. maritima maritima). Using an eco-phylogeographic lens within a conservation management context, we tested predictions about species' degree of evolutionary history and ecological specialization to tidal marshes, habitat, current distribution and population status on population genetic metrics. Population structure differed among the species consistent with their current distribution and habitat factors, rather than degree of ecological specialization: seaside sparrows were panmictic, saltmarsh sparrows showed hierarchical structure and Nelson's sparrows were differentiated into multiple, genetically distinct populations. Neutral population genetic theory and demographic/evolutionary history predicted patterns of genetic diversity and Ne rather than degree of ecological specialization. Patterns of population variation and evolutionary distinctiveness (Shapely metric) suggest different conservation measures for long-term persistence and evolutionary potential in each species. Our findings contribute to a broader understanding of the complex factors influencing genetic variation, beyond specialist-generalist status and support the role of an eco-phylogeographic approach in population and conservation genetics.

The relationship between human blood metabolites and preeclampsia-eclampsia: A Mendelian randomization study

Preeclampsia and eclampsia are serious complications of pregnancy, leading to high rates of maternal and neonatal mortality. During pregnancy, there are changes in relevant serum metabolites in women. However, it remains unclear if these serum metabolites contribute to the development of associated disorders during pregnancy. Therefore, we conducted a Mendelian randomization study to explore the causal relationship between serum metabolites and preeclampsia and eclampsia. We utilized the inverse variance weighted model as our primary analysis approach. We complemented this with sensitivity analyses, including the heterogeneity test, horizontal pleiotropy test, and leave-one-out analysis, to ensure the robustness of our findings. Furthermore, we conducted linkage disequilibrium score regression, multivariable Mendelian randomization, and metabolic pathway analysis to further explore the genetic data. The Mendelian randomization analysis has identified γ-glutamylglutamine, inosine, and isoleucine 10 metabolites that are significantly associated with preeclampsia, and γ-glutamylglutamine and phenylacetate 8 metabolites that may potentially contribute to the development of eclampsia. Notably, γ-glutamylglutamine has been found to have a causal relationship with both preeclampsia and eclampsia. In the multivariable Mendelian randomization analysis, our research findings suggest that both isoleucine and X-14304-leucylalanine directly impact preeclampsia within the context of amino acids and peptides. Moreover, our observations reveal that carbohydrates can also have a direct effect on preeclampsia. Importantly, it should be emphasized that only 3-lactate in amino acids has been shown to have a direct influence on eclampsia. This research has the potential to enhance our understanding of the biological variances related to disease status, providing a foundation for future investigations.

Characterization of Vestibular Phenotypes in Patients with Genetic Hearing Loss

Background: The vestibular phenotypes of patients with genetic hearing loss are poorly understood. Methods: we performed genetic testing including exome sequencing and vestibular function tests to investigate vestibular phenotypes and functions in patients with genetic hearing loss. Results: Among 627 patients, 143 (22.8%) had vestibular symptoms. Genetic variations were confirmed in 45 (31.5%) of the 143 patients. Nineteen deafness genes were linked with vestibular symptoms; the most frequent genes in autosomal dominant and recessive individuals were COCH and SLC26A4, respectively. Vestibular symptoms were mostly of the vertigo type, recurrent, and persisted for hours in the genetically confirmed and unconfirmed groups. Decreased vestibular function in the caloric test, video head impulse test, cervical vestibular-evoked myogenic potential, and ocular vestibular-evoked myogenic potential was observed in 42.0%, 16.3%, 57.8%, and 85.0% of the patients, respectively. The caloric test revealed a significantly higher incidence of abnormal results in autosomal recessive individuals than in autosomal dominant individuals (p = 0.011). The genes, including SLC26A4, COCH, KCNQ4, MYH9, NLRP3, EYA4, MYO7A, MYO15A, and MYH9, were heterogeneously associated with abnormalities in the vestibular function test. Conclusions: In conclusion, diverse vestibular symptoms are commonly concomitant with genetic hearing loss and are easily overlooked.

Serial passage of PDCoV in cell culture reduces its pathogenicity and its damage of gut microbiota homeostasis in piglets

Porcine deltacoronavirus (PDCoV) is an enteropathogenic coronavirus that mainly causes diarrhea in suckling piglets, and also has the potential for cross-species transmission. However, there are still no commercial vaccines available to prevent and control PDCoV infection. In this study, PDCoV strain HNZK-02 was serially propagated in vitro for up to 150 passages and the amino acid changes have mainly occurred in the S protein during serial passage which caused structure change. PDCoV HNZK-02-passage 5 (P5)-infected piglets exhibited acute and severe watery diarrhea, an obvious intestinal damage, while the piglets infected with PDCoV HNZK-02-P150 showed no obvious clinical signs, weak intestinal lesions, and lower viral loads in rectal swabs and various tissues. Compared with the PDCoV HNZK-02-P5 infection, HNZK-02-P150 infection resulted in a decrease in intestinal mucosal permeability and pro-inflammatory cytokines. Moreover, PDCoV HNZK-02-P5 infection had significantly reduced bacterial diversity and increased relative abundance of opportunistic pathogens, while PDCoV HNZK-02-P150 infection did not significantly affect the bacterial diversity, and the relative abundance of probiotics increased. Furthermore, the alterations of gut microbiota were closely related to the change of pro-inflammatory factor. Metagenomics prediction analysis demonstrated that HNZK-02-P150 modulated the tyrosine metabolism, Nucleotide-binding and oligomerization domain (NOD)-like receptor signaling pathway, and lipopolysaccharide biosynthesis, which coincided with lower inflammatory response and intestinal permeability in the piglets infected with HNZK-02-P150. In conclusion, the PDCoV HNZK-02 was successfully attenuated by serial passage in vitro, and the changes of S gene, metabolic function, and gut microbiota may contribute to the attenuation. The PDCoV HNZK-02-P150 may have the potential for developing live-attenuated vaccine.IMPORTANCEPorcine deltacoronavirus (PDCoV) is an enteropathogen causing severe diarrhea, dehydration, and death in nursing piglets, devastating great economic losses for the global swine industry, and has cross-species transmission and zoonotic potential. There are currently no approved treatments or vaccines available for PDCoV. In addition, gut microbiota has an important relationship with the development of many diseases. Here, the PDCoV virulent HNZK-02 strain was successfully attenuated by serial passage on cell cultures, and the pathogenesis and effects on the gut microbiota composition and metabolic function of the PDCoV HNZK-02-P5 and P150 strains were investigated in piglets. We also found the genetic changes in the S protein during passage in vitro and the gut microbiota may contribute to the pathogenesis of PDCoV, while their interaction molecular mechanism would need to be explored further.

Dynamics of extrachromosomal circular DNA in rice

The genome's dynamic nature, exemplified by elements like extrachromosomal circular DNA (eccDNA), is crucial for biodiversity and adaptation. Yet, the role of eccDNA in plants, particularly rice, remains underexplored. Here, we identify 25,598 eccDNAs, unveiling the widespread presence of eccDNA across six rice tissues and revealing its formation as a universal and random process. Interestingly, we discover that direct repeats play a pivotal role in eccDNA formation, pointing to a unique origin mechanism. Despite eccDNA's prevalence in coding sequences, its impact on gene expression is minimal, implying its roles beyond gene regulation. We also observe the association between eccDNA's formation and minor chromosomal deletions, providing insights of its possible function in regulating genome stability. Further, we discover eccDNA specifically accumulated in rice leaves, which may be associated with DNA damage caused by environmental stressors like intense light. In summary, our research advances understanding of eccDNA's role in the genomic architecture and offers valuable insights for rice cultivation and breeding.

Dominance reversals: the resolution of genetic conflict and maintenance of genetic variation

Beneficial reversals of dominance reduce the costs of genetic trade-offs and can enable selection to maintain genetic variation for fitness. Beneficial dominance reversals are characterized by the beneficial allele for a given context (e.g. habitat, developmental stage, trait or sex) being dominant in that context but recessive where deleterious. This context dependence at least partially mitigates the fitness consequence of heterozygotes carrying one non-beneficial allele for their context and can result in balancing selection that maintains alternative alleles. Dominance reversals are theoretically plausible and are supported by mounting empirical evidence. Here, we highlight the importance of beneficial dominance reversals as a mechanism for the mitigation of genetic conflict and review the theory and empirical evidence for them. We identify some areas in need of further research and development and outline three methods that could facilitate the identification of antagonistic genetic variation (dominance ordination, allele-specific expression and allele-specific ATAC-Seq (assay for transposase-accessible chromatin with sequencing)). There is ample scope for the development of new empirical methods as well as reanalysis of existing data through the lens of dominance reversals. A greater focus on this topic will expand our understanding of the mechanisms that resolve genetic conflict and whether they maintain genetic variation.

Gene expression plasticity followed by genetic change during colonization in a high-elevation environment

Phenotypic plasticity facilitates organismal invasion of novel environments, and the resultant phenotypic change may later be modified by genetic change, so called 'plasticity first.' Herein, we quantify gene expression plasticity and regulatory adaptation in a wild bird (Eurasian Tree Sparrow) from its original lowland (ancestral stage), experimentally implemented hypoxia acclimation (plastic stage), and colonized highland (colonized stage). Using a group of co-expressed genes from the cardiac and flight muscles, respectively, we demonstrate that gene expression plasticity to hypoxia tolerance is more often reversed than reinforced at the colonized stage. By correlating gene expression change with muscle phenotypes, we show that colonized tree sparrows reduce maladaptive plasticity that largely associated with decreased hypoxia tolerance. Conversely, adaptive plasticity that is congruent with increased hypoxia tolerance is often reinforced in the colonized tree sparrows. Genes displaying large levels of reinforcement or reversion plasticity (i.e. 200% of original level) show greater genetic divergence between ancestral and colonized populations. Overall, our work demonstrates that gene expression plasticity at the initial stage of high-elevation colonization can be reversed or reinforced through selection-driven adaptive modification.

A new hybrid incompatibility locus between Drosophila melanogaster and Drosophila sechellia

Despite the fundamental importance of hybrid incompatibilities to the process of speciation, there are few cases where the evolution and genetic architecture of hybrid incompatibilities are understood. One of the longest studied hybrid incompatibilities causes F1 hybrid male inviability in crosses between Drosophila melanogaster females and males from the Drosophila simulans clade of species-Drosophila simulans, Drosophila mauritiana, and Drosophila sechellia. Here, we discover dramatic differences in the manifestation of this lethal hybrid incompatibility among the D. simulans clade of species. In particular, F1 hybrid males between D. melanogaster and D. sechellia are resistant to hybrid rescue through RNAi knockdown of an essential hybrid incompatibility gene. To understand the genetic basis of this inter-species difference in hybrid rescue, we developed a triple-hybrid mapping method. Our results show that 2 discrete large effect loci and many dispersed small effect changes across the genome underlie D. sechellia aversion to hybrid rescue. The large effect loci encompass a known incompatibility gene Lethal hybrid rescue (Lhr) and previously unknown factor, Sechellia aversion to hybrid rescue (Satyr). These results show that the genetic architecture of F1 hybrid male inviability is overlapping but not identical in the 3 inter-species crosses. Our results raise questions about whether new hybrid incompatibility genes can integrate into an existing hybrid incompatibility thus increasing in complexity over time, or if the continued evolution of genes can gradually strengthen an existing hybrid incompatibility.

The prevalence and linkage disequilibrium of 21 genetic variations related to thrombophilia, folate cycle, and hypertension in reproductive age women of Rostov region (Russia)

Several maternal genetic variations are known to play an important role during pregnancy since they can affect mother health and/or fetal growth. The frequency of these variants is variable among different populations. This study aimed to investigate thrombophilia, folate metabolism and hypertension genetic variants in reproductive age women of Rostov region (Russia) and then assess their linkage disequilibrium (LD) and heterogeneity among populations. A total of 3108 reproductive age women were included (33.75 ± 5.13 years). Twenty-one genetic variants were detected with RT-PCR. LD was tested according to (D') coefficient and p value. The highest frequency of mutant allele in studied population was as follows: PAI-1 rs1799768, MTRR rs1801394, AGT rs699, and AGTR2 rs1403543. We showed a high possibility of coinheritance of MTHFR rs1801133 with rs1801131 and AGT rs699 with rs4762 (D'=0.992 and 0.999, respectively). In addition, comparative analysis showed F7 rs6046, FGB rs1800790, MTR rs1805087, and AGT rs699 significantly more frequent among Rostov females by 1.3-1.5 times than European. MTHFR rs1801133, ADD1 rs4961, AGTR2 rs1403543, NOS3 rs2070744, and rs1799983 were with higher frequencies in Europeans than those in the studied group. Our data could be used as a reference for further associative studies of targeted genetic variations in different pregnancy complications specifically in this population.

A genomic toolkit for winged bean Psophocarpus tetragonolobus

A sustainable supply of plant protein is critical for future generations and needs to be achieved while reducing green house gas emissions from agriculture and increasing agricultural resilience in the face of climate volatility. Agricultural diversification with more nutrient-rich and stress tolerant crops could provide the solution. However, this is often hampered by the limited availability of genomic resources and the lack of understanding of the genetic structure of breeding germplasm and the inheritance of important traits. One such crop with potential is winged bean (Psophocarpus tetragonolobus), a high seed protein tropical legume which has been termed 'the soybean for the tropics'. Here, we present a chromosome level winged bean genome assembly, an investigation of the genetic diversity of 130 worldwide accessions, together with two linked genetic maps and a trait QTL analysis (and expression studies) for regions of the genome with desirable ideotype traits for breeding, namely architecture, protein content and phytonutrients.

Effects of genetic variation on the structure of RNA and protein

Changes in the structure of RNA and protein, have an important impact on biological functions and are even important determinants of disease pathogenesis and treatment. Some genetic variations, including copy number variation, single nucleotide variation, and so on, can lead to changes in biological function and increased susceptibility to certain diseases by changing the structure of RNA or protein. With the development of structural biology and sequencing technology, a large amount of RNA and protein structure data and genetic variation data resources has emerged to be used to explain biological processes. Here, we reviewed the effects of genetic variation on the structure of RNAs and proteins, and investigated their impact on several diseases. An online resource (http://www.onethird-lab.com/gems/) to support convenient retrieval of common tools is also built. Finally, the challenges and future development of the effects of genetic variation on RNA and protein were discussed.

Are Genetics the Predicting Factor for the Success of Migraine Surgery? A Report on Identical Twins

Migraine affects more than 1 billion people globally, with distinct genetic variations influencing susceptibility. Thereby, genetic variations play a key role in the probability of developing migraine. However, personalized genetic analysis-based treatment options in migraine treatments are limited. Notably, surgical deactivation of extracranial trigger has shown efficacy in the treatment of migraine patients with identifiable trigger points in specific anatomical locations in the head and neck region. We present the first case of monozygotic twin sisters, both experiencing occipital and temporal-triggered migraine headaches with identical history and characteristics and without response to conservative migraine treatments. After surgical intervention, targeting the greater and lesser occipital nerves as well as auriculotemporal nerves, both twin sisters exhibited an over 99% reduction in symptoms without postoperative complications. This case suggests a potential correlation between genetic background, irrespective of environmental factors, and the effectiveness of surgical deactivation of trigger points in migraine management.

Genomic data in the All of Us Research Program

Comprehensively mapping the genetic basis of human disease across diverse individuals is a long-standing goal for the field of human genetics1-4. The All of Us Research Program is a longitudinal cohort study aiming to enrol a diverse group of at least one million individuals across the USA to accelerate biomedical research and improve human health5,6. Here we describe the programme's genomics data release of 245,388 clinical-grade genome sequences. This resource is unique in its diversity as 77% of participants are from communities that are historically under-represented in biomedical research and 46% are individuals from under-represented racial and ethnic minorities. All of Us identified more than 1 billion genetic variants, including more than 275 million previously unreported genetic variants, more than 3.9 million of which had coding consequences. Leveraging linkage between genomic data and the longitudinal electronic health record, we evaluated 3,724 genetic variants associated with 117 diseases and found high replication rates across both participants of European ancestry and participants of African ancestry. Summary-level data are publicly available, and individual-level data can be accessed by researchers through the All of Us Researcher Workbench using a unique data passport model with a median time from initial researcher registration to data access of 29 hours. We anticipate that this diverse dataset will advance the promise of genomic medicine for all.

Genetic variation in TAS2R38 bitterness receptor is associated with body composition in Korean females

Bitterness-receptor gene TAS2R38 is associated with taste sensitivity and dietary behaviour, and is known to play a critical role in adiposity. However, evidence regarding body composition from a large cohort is lacking. This study aimed to ascertain whether TAS2R38 rs10246939 C > T bitterness genetic variation is associated with body composition in Korean individuals. The TAS2R38 rs10246939 genotypes, anthropometric measurements, and body composition of 1,843 males and 1,801 females from the Korean Genome and Epidemiology Study were analysed. Findings suggested that there was a significant difference in body fat components by TAS2R38 genotype. Furthermore, the bitterness genotype exhibited a positive association with adiposity markers in females. The TT genotype showed greater body mass index, body fat percentage, and degree of obesity than those with the C allele. However, such an association was not observed in males. In conclusion, this study suggests that TAS2R38 rs10246939 is associated with fat tissue markers in Korean females.

Genetic constraints in genes exhibiting splicing plasticity in facultative diapause

Phenotypic plasticity is produced and maintained by processes regulating the transcriptome. While differential gene expression is among the most important of these processes, relatively little is known about other sources of transcriptional variation. Previous work suggests that alternative splicing plays an extensive and functionally unique role in transcriptional plasticity, though plastically spliced genes may be more constrained than the remainder of expressed genes. In this study, we explore the relationship between expression and splicing plasticity, along with the genetic diversity in those genes, in an ecologically consequential polyphenism: facultative diapause. Using 96 samples spread over two tissues and 10 timepoints, we compare the extent of differential splicing and expression between diapausing and direct developing pupae of the butterfly Pieris napi. Splicing differs strongly between diapausing and direct developing trajectories but alters a smaller and functionally unique set of genes compared to differential expression. We further test the hypothesis that among these expressed loci, plastically spliced genes are likely to experience the strongest purifying selection to maintain seasonally plastic phenotypes. Genes with unique transcriptional changes through diapause consistently had the lowest nucleotide diversity, and this effect was consistently stronger among genes that were differentially spliced compared to those with just differential expression through diapause. Further, the strength of negative selection was higher in the population expressing diapause every generation. Our results suggest that maintenance of the molecular mechanisms involved in diapause progression, including post-transcriptional modifications, are highly conserved and likely to experience genetic constraints, especially in northern populations of P. napi.

Gene expression and alternative splicing contribute to adaptive divergence of ecotypes

Regulation of gene expression is a critical link between genotype and phenotype explaining substantial heritable variation within species. However, we are only beginning to understand the ways that specific gene regulatory mechanisms contribute to adaptive divergence of populations. In plants, the post-transcriptional regulatory mechanism of alternative splicing (AS) plays an important role in both development and abiotic stress response, making it a compelling potential target of natural selection. AS allows organisms to generate multiple different transcripts/proteins from a single gene and thus may provide a source of evolutionary novelty. Here, we examine whether variation in alternative splicing and gene expression levels might contribute to adaptation and incipient speciation of dune-adapted prairie sunflowers in Great Sand Dunes National Park, Colorado, USA. We conducted a common garden experiment to assess transcriptomic variation among ecotypes and analyzed differential expression, differential splicing, and gene coexpression. We show that individual genes are strongly differentiated for both transcript level and alternative isoform proportions, even when grown in a common environment, and that gene coexpression networks are disrupted between ecotypes. Furthermore, we examined how genome-wide patterns of sequence divergence correspond to divergence in transcript levels and isoform proportions and find evidence for both cis and trans-regulation. Together, our results emphasize that alternative splicing has been an underappreciated mechanism providing source material for natural selection at short evolutionary time scales.

Genome-Enabled Analysis of Population Dynamics and Virulence-Associated Loci in the Oat Crown Rust Fungus Puccinia coronata f. sp. avenae

Puccinia coronata f. sp. avenae (Pca) is an important fungal pathogen causing crown rust that impacts oat production worldwide. Genetic resistance for crop protection against Pca is often overcome by the rapid virulence evolution of the pathogen. This study investigated the factors shaping adaptive evolution of Pca using pathogen populations from distinct geographic regions within the United States and South Africa. Phenotypic and genome-wide sequencing data of these diverse Pca collections, including 217 isolates, uncovered phylogenetic relationships and established distinct genetic composition between populations from northern and southern regions from the United States and South Africa. The population dynamics of Pca involve a bidirectional movement of inoculum between northern and southern regions of the United States and contributions from clonality and sexuality. The population from South Africa is solely clonal. A genome-wide association study (GWAS) employing a haplotype-resolved Pca reference genome was used to define 11 virulence-associated loci corresponding to 25 oat differential lines. These regions were screened to determine candidate Avr effector genes. Overall, the GWAS results allowed us to identify the underlying genetic factors controlling pathogen recognition in an oat differential set used in the United States to assign pathogen races (pathotypes). Key GWAS findings support complex genetic interactions in several oat lines, suggesting allelism among resistance genes or redundancy of genes included in the differential set, multiple resistance genes recognizing genetically linked Avr effector genes, or potentially epistatic relationships. A careful evaluation of the composition of the oat differential set accompanied by the development or implementation of molecular markers is recommended. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Evaluating the clinical validity of genes related to hemostasis and thrombosis using the Clinical Genome Resource gene curation framework

BACKGROUND

Inherited bleeding, thrombotic, and platelet disorders (BTPDs) are a heterogeneous set of diseases, many of which are very rare globally. Over the past 5 decades, the genetic basis of some of these disorders has been identified, and recently, high-throughput sequencing has become the primary means of identifying disease-causing genetic variants.

OBJECTIVES

Knowledge of the clinical validity of a gene-disease relationship is essential to provide an accurate diagnosis based on results of diagnostic gene panel tests and inform the construction of such panels. The Scientific and Standardization Committee for Genetics in Thrombosis and Hemostasis undertook a curation process for selecting 96 TIER1 genes for BTPDs. The purpose of the process was to evaluate the evidence supporting each gene-disease relationship and provide an expert-reviewed classification for the clinical validity of genes associated with BTPDs.

METHODS

The Clinical Genome Resource (ClinGen) Hemostasis/Thrombosis Gene Curation Expert Panel assessed the strength of evidence for TIER1 genes using the semiquantitative ClinGen gene-disease clinical validity framework. ClinGen Lumping and Splitting guidelines were used to determine the appropriate disease entity or entities for each gene, and 101 gene-disease relationships were identified for curation.

RESULTS

The final outcome included 68 Definitive (67%), 26 Moderate (26%), and 7 Limited (7%) classifications. The summary of each curation is available on the ClinGen website.

CONCLUSION

Expert-reviewed assignment of gene-disease relationships by the ClinGen Hemostasis/Thrombosis Gene Curation Expert Panel facilitates accurate molecular diagnoses of BTPDs by clinicians and diagnostic laboratories. These curation efforts can allow genetic testing to focus on genes with a validated role in disease.

Genetic Analyses of Flower Main Traits from Two Pitayas and Their Progenies: A Cactus Plant

Elucidation of the genetic foundation governing crucial traits in pitaya flowers is imperative for enhancing both the ornamental and economic values. In this study, the dynamic variation in flower genetics, segregation variation patterns, and a mixed inheritance model of the major and multigene flower traits of 'Dahong' and 'Honghuaqinglong' pitayas and their progenies were explored. The results showed that the main traits of flowers exhibited varying degrees of variation among the reciprocal F1 hybrids, with the data exhibiting the characteristics of quantitative traits. The betalain content, petal number, and stigma number exhibited values below the median values of the parents, suggesting a genetic inclination towards lower values. Perianth width, calyx tube width, petal number, and stigma number had the same genetic effects and significant correlation. Stigma-related traits had a clear maternal inheritance tendency. The heritability of flower length, stigma relative to anther distance, and petal betalain content was governed by two pairs of additive-dominant major genes. Perianth width, calyx tube width, petal number, and stigma number all conformed to the model of two pairs of equal-additive-dominant major genes. This study provides valuable information for parental selection, cross-breeding, and the enhancement of pitaya varieties to meet market preferences and environmental conditions.

Protein-folding chaperones predict structure-function relationships and cancer risk in BRCA1 mutation carriers

Predicting the risk of cancer mutations is critical for early detection and prevention, but differences in allelic severity of human carriers confound risk predictions. Here, we elucidate protein folding as a cellular mechanism driving differences in mutation severity of tumor suppressor BRCA1. Using a high-throughput protein-protein interaction assay, we show that protein-folding chaperone binding patterns predict the pathogenicity of variants in the BRCA1 C-terminal (BRCT) domain. HSP70 selectively binds 94% of pathogenic BRCA1-BRCT variants, most of which engage HSP70 more than HSP90. Remarkably, the magnitude of HSP70 binding linearly correlates with loss of folding and function. We identify a prevalent class of human hypomorphic BRCA1 variants that bind moderately to chaperones and retain partial folding and function. Furthermore, chaperone binding signifies greater mutation penetrance and earlier cancer onset in the clinic. Our findings demonstrate the utility of chaperones as quantitative cellular biosensors of variant folding, phenotypic severity, and cancer risk.

One-pot DTECT enables rapid and efficient capture of genetic signatures for precision genome editing and clinical diagnostics

The detection of genomic sequences and their alterations is crucial for basic research and clinical diagnostics. However, current methodologies are costly and time-consuming and require outsourcing sample preparation, processing, and analysis to genomic companies. Here, we establish One-pot DTECT, a platform that expedites the detection of genetic signatures, only requiring a short incubation of a PCR product in an optimized one-pot mixture. One-pot DTECT enables qualitative, quantitative, and visual detection of biologically relevant variants, such as cancer mutations, and nucleotide changes introduced by prime editing and base editing into cancer cells and human primary T cells. Notably, One-pot DTECT achieves quantification accuracy for targeted genetic signatures comparable with Sanger and next-generation sequencing. Furthermore, its effectiveness as a diagnostic platform is demonstrated by successfully detecting sickle cell variants in blood and saliva samples. Altogether, One-pot DTECT offers an efficient, versatile, adaptable, and cost-effective alternative to traditional methods for detecting genomic signatures.

Evidence review and considerations for use of first line genome sequencing to diagnose rare genetic disorders

Early use of genome sequencing (GS) in the diagnostic odyssey can reduce suffering and improve care, but questions remain about which patient populations are most amenable to GS as a first-line diagnostic test. To address this, the Medical Genome Initiative conducted a literature review to identify appropriate clinical indications for GS. Studies published from January 2011 to August 2022 that reported on the diagnostic yield (DY) or clinical utility of GS were included. An exploratory meta-analysis using a random effects model evaluated DY based on cohort size and diagnosed cases per cohort. Seventy-one studies met inclusion criteria, comprising over 13,000 patients who received GS in one of the following settings: hospitalized pediatric patients, pediatric outpatients, adult outpatients, or mixed. GS was the first-line test in 38% (27/71). The unweighted mean DY of first-line GS was 45% (12-73%), 33% (6-86%) in cohorts with prior genetic testing, and 33% (9-60%) in exome-negative cohorts. Clinical utility was reported in 81% of first-line GS studies in hospitalized pediatric patients. Changes in management varied by cohort and underlying molecular diagnosis (24-100%). To develop evidence-informed points to consider, the quality of all 71 studies was assessed using modified American College of Radiology (ACR) criteria, with five core points to consider developed, including recommendations for use of GS in the N/PICU, in lieu of sequential testing and when disorders with substantial allelic heterogeneity are suspected. Future large and controlled studies in the pediatric and adult populations may support further refinement of these recommendations.

Exploring the Genetic Landscape of Mild Behavioral Impairment as an Early Marker of Cognitive Decline: An Updated Review Focusing on Alzheimer's Disease

The clinical features and pathophysiology of neuropsychiatric symptoms (NPSs) in dementia have been extensively studied. However, the genetic architecture and underlying neurobiological mechanisms of NPSs at preclinical stages of cognitive decline and Alzheimer's disease (AD) remain largely unknown. Mild behavioral impairment (MBI) represents an at-risk state for incident cognitive impairment and is defined by the emergence of persistent NPSs among non-demented individuals in later life. These NPSs include affective dysregulation, decreased motivation, impulse dyscontrol, abnormal perception and thought content, and social inappropriateness. Accumulating evidence has recently begun to shed more light on the genetic background of MBI, focusing on its potential association with genetic factors related to AD. The Apolipoprotein E (APOE) genotype and the MS4A locus have been associated with affective dysregulation, ZCWPW1 with social inappropriateness and psychosis, BIN1 and EPHA1 with psychosis, and NME8 with apathy. The association between MBI and polygenic risk scores (PRSs) in terms of AD dementia has been also explored. Potential implicated mechanisms include neuroinflammation, synaptic dysfunction, epigenetic modifications, oxidative stress responses, proteosomal impairment, and abnormal immune responses. In this review, we summarize and critically discuss the available evidence on the genetic background of MBI with an emphasis on AD, aiming to gain insights into the potential underlying neurobiological mechanisms, which till now remain largely unexplored. In addition, we propose future areas of research in this emerging field, with the aim to better understand the molecular pathophysiology of MBI and its genetic links with cognitive decline.

The HSP90AA1 Gene Is Involved in Heat Stress Responses and Its Functional Genetic Polymorphism Are Associated with Heat Tolerance in Holstein Cows

As the stress-inducible isoform of the Heat Shock Protein 90 (HSP90), the HSP90AA1 gene encodes HSP90α and plays an important role in heat stress (HS) response. Therefore, this study aimed to investigate the role of the HSP90AA1 gene in cellular responses during HS and to identify functional single nucleotide polymorphisms (SNP) associated with thermotolerance in Holstein cattle. For the in vitro validation experiment of acute HS, cells from the Madin-Darby bovine kidney (MDBK) cell line were exposed to 42°C for 1 h, and various parameters were assessed, including cell apoptosis, cell autophagy, and the cellular functions of HSP90α by using its inhibitor 17-allylamino-17-demethoxygeldanamycin (17-AAG). Furthermore, the polymorphisms identified in the HSP90AA1 gene and their functions related to HS were in vitro validated. Acute HS exposure induced cell apoptosis, cell autophagy, and upregulated expression of the HSP90AA1 gene. Inhibition of HSP90α by 17-AAG treatment had a significant effect on the expression of the HSP90α protein (P < 0.05) and increased cell apoptosis. However, autophagy decreased in comparison to the control treatment when cells were exposed to 42°C for 1 h. Five SNPs identified in the HSP90AA1 gene were significantly associated with rectal temperature (RT; P < 0.05) and respiration score (RS; P < 0.05) in Holstein cows, in which the rs109256957 SNP is located in the 3' untranslated region (3' UTR). Furthermore, we demonstrated that the 3' UTR of HSP90AA1 is a direct target of bta-miR-1224 by cell transfection with exogenous miRNA mimic and inhibitor. The luciferase assays revealed that the SNP rs109256957 affects the regulation of bta-miR-1224 binding activity and alters the expression of the HSP90AA1 gene. Heat stress-induced HSP90AA1 expression maintains cell survival by inhibiting cell apoptosis and increasing cell autophagy. The rs109256957 SNP located in the 3' UTR region is a functional variation and it affects the HSP90AA1 expression by altering its binding activity with bta-miR-1224, thereby associating with the physiological parameters of Holstein cows.

Variant reclassification and clinical implications

Genomic technologies have transformed clinical genetic testing, underlining the importance of accurate molecular genetic diagnoses. Variant classification, ranging from benign to pathogenic, is fundamental to these tests. However, variant reclassification, the process of reassigning the pathogenicity of variants over time, poses challenges to diagnostic legitimacy. This review explores the medical and scientific literature available on variant reclassification, focusing on its clinical implications.Variant reclassification is driven by accruing evidence from diverse sources, leading to variant reclassification frequency ranging from 3.6% to 58.8%. Recent studies have shown that significant changes can occur when reviewing variant classifications within 1 year after initial classification, illustrating the importance of early, accurate variant assignation for clinical care.Variants of uncertain significance (VUS) are particularly problematic. They lack clear categorisation but have influenced patient treatment despite recommendations against it. Addressing VUS reclassification is essential to enhance the credibility of genetic testing and the clinical impact. Factors affecting reclassification include standardised guidelines, clinical phenotype-genotype correlations through deep phenotyping and ancestry studies, large-scale databases and bioinformatics tools. As genomic databases grow and knowledge advances, reclassification rates are expected to change, reducing discordance in future classifications.Variant reclassification affects patient diagnosis, precision therapy and family screening. The exact patient impact is yet unknown. Understanding influencing factors and adopting standardised guidelines are vital for precise molecular genetic diagnoses, ensuring optimal patient care and minimising clinical risk.

A genetic variant in the immune-related gene ERAP1 affects colorectal cancer prognosis

BACKGROUND

Findings on the association of genetic factors and colorectal cancer (CRC) survival are limited and inconsistent, and revealing the mechanism underlying their prognostic roles is of great importance. This study aimed to explore the relationship between functional genetic variations and the prognosis of CRC and further reveal the possible mechanism.

METHODS

We first systematically performed expression quantitative trait locus (eQTL) analysis using The Cancer Genome Atlas (TCGA) dataset. Then, the Kaplan-Meier analysis was used to filter out the survival-related eQTL target genes of CRC patients in two public datasets (TCGA and GSE39582 dataset from the Gene Expression Omnibus database). The seven most potentially functional eQTL single nucleotide polymorphisms (SNPs) associated with six survival-related eQTL target genes were genotyped in 907 Chinese CRC patients with clinical prognosis data. The regulatory mechanism of the survival-related SNP was further confirmed by functional experiments.

RESULTS

The rs71630754 regulating the expression of endoplasmic reticulum aminopeptidase 1 ( ERAP1 ) was significantly associated with the prognosis of CRC (additive model, hazard ratio [HR]: 1.43, 95% confidence interval [CI]: 1.08-1.88, P = 0.012). The results of dual-luciferase reporter assay and electrophoretic mobility shift assay showed that the A allele of the rs71630754 could increase the binding of transcription factor 3 (TCF3) and subsequently reduce the expression of ERAP1 . The results of bioinformatic analysis showed that lower expression of ERAP1 could affect the tumor immune microenvironment and was significantly associated with severe survival outcomes.

CONCLUSION

The rs71630754 could influence the prognosis of CRC patients by regulating the expression of the immune-related gene ERAP1 .

TRIAL REGISTRATION

No. NCT00454519 ( https://clinicaltrials.gov/ ).

Cepharanthine analogs mining and genomes of Stephania accelerate anti-coronavirus drug discovery

Cepharanthine is a secondary metabolite isolated from Stephania. It has been reported that it has anti-conronaviruses activities including severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Here, we assemble three Stephania genomes (S. japonica, S. yunnanensis, and S. cepharantha), propose the cepharanthine biosynthetic pathway, and assess the antiviral potential of compounds involved in the pathway. Among the three genomes, S. japonica has a near telomere-to-telomere assembly with one remaining gap, and S. cepharantha and S. yunnanensis have chromosome-level assemblies. Following by biosynthetic gene mining and metabolomics analysis, we identify seven cepharanthine analogs that have broad-spectrum anti-coronavirus activities, including SARS-CoV-2, Guangxi pangolin-CoV (GX_P2V), swine acute diarrhoea syndrome coronavirus (SADS-CoV), and porcine epidemic diarrhea virus (PEDV). We also show that two other genera, Nelumbo and Thalictrum, can produce cepharanthine analogs, and thus have the potential for antiviral compound discovery. Results generated from this study could accelerate broad-spectrum anti-coronavirus drug discovery.

Natural selection and genetic diversity maintenance in a parasitic wasp during continuous biological control application

Aphidius gifuensis is a parasitoid wasp and primary endoparasitoid enemy of the peach potato aphid, Myzus persicae. Artificially reared, captive wasps of this species have been extensively and effectively used to control populations of aphids and limit crop loss. However, the consequences of large-scale releasing of captive A. gifuensis, such as genetic erosion and reduced fitness in wild populations of this species, remains unclear. Here, we sequence the genomes of 542 A. gifuensis individuals collected across China, including 265 wild and 277 human-intervened samples. Population genetic analyses on wild individuals recovered Yunnan populations as the ancestral group with the most complex genetic structure. We also find genetic signature of environmental adaptation during the dispersal of wild populations from Yunnan to other regions. While comparative genomic analyses of captive wasps revealed a decrease in genetic diversity during long-term rearing, population genomic analyses revealed signatures of natural selection by several biotic (host plants) or abiotic (climate) factors, which support maintenance of the gene pool of wild populations in spite of the introduction of captive wasps. Therefore, the impact of large-scale release is reduced. Our study suggests that A. gifuensis is a good system for exploring the genetic and evolutionary effects of mass rearing and release on species commonly used as biocontrol agents.

Integration of risk factor polygenic risk score with disease polygenic risk score for disease prediction

Polygenic risk score (PRS) is useful for capturing an individual's genetic susceptibility. However, previous studies have not fully exploited the potential of the risk factor PRS (RFPRS) for disease prediction. We explored the potential of integrating disease-related RFPRSs with disease PRS to enhance disease prediction performance. We constructed 112 RFPRSs and analyzed the association of RFPRSs with diseases to identify disease-related RFPRSs in 700 diseases, using the UK Biobank dataset. We uncovered 6157 statistically significant associations between 247 diseases and 109 RFPRSs. We estimated the disease PRSs of 70 diseases that exhibited statistically significant heritability, to generate RFDiseasemetaPRS-a combined PRS integrating RFPRSs and disease PRS-and compare the prediction performance metrics between RFDiseasemetaPRS and disease PRS. RFDiseasemetaPRS showed better performance for Nagelkerke's pseudo-R2, odds ratio (OR) per 1 SD, net reclassification improvement (NRI) values and difference of R2 considered by variance of R2 in 31 out of 70 diseases. Additionally, we assessed risk classification between two models by examining OR between the top 10% and remaining 90% individuals for the 31 diseases; RFDiseasemetaPRS exhibited better R2, NRI and OR than disease PRS. These findings highlight the importance of utilizing RFDiseasemetaPRS, which can provide personalized healthcare and tailored prevention strategies.

The Allen Ancient DNA Resource (AADR) a curated compendium of ancient human genomes

More than two hundred papers have reported genome-wide data from ancient humans. While the raw data for the vast majority are fully publicly available testifying to the commitment of the paleogenomics community to open data, formats for both raw data and meta-data differ. There is thus a need for uniform curation and a centralized, version-controlled compendium that researchers can download, analyze, and reference. Since 2019, we have been maintaining the Allen Ancient DNA Resource (AADR), which aims to provide an up-to-date, curated version of the world's published ancient human DNA data, represented at more than a million single nucleotide polymorphisms (SNPs) at which almost all ancient individuals have been assayed. The AADR has gone through six public releases at the time of writing and review of this manuscript, and crossed the threshold of >10,000 individuals with published genome-wide ancient DNA data at the end of 2022. This note is intended as a citable descriptor of the AADR.

Regional environmental differences significantly affect the genetic structure and genetic differentiation of Carpinus tientaiensis Cheng, an endemic and extremely endangered species from China

Differences in topography and environment greatly affect the genetic structure and genetic differentiation of species, and endemic or endangered species with limited geographic ranges seem to be more sensitive to changes in climate and other environmental factors. The complex topography of eastern China is likely to affect genetic differentiation of plants there. Carpinus tientaiensis Cheng is a native and endangered plants from China, and exploring its genetic diversity has profound significance for protection and the collection of germplasm resources. Based on AFLP markers, this study found that C. tientaiensis has low genetic diversity, which mainly came from within populations, while Shangshantou and Tiantai Mountain populations have relatively high genetic diversity. The Nei genetic distance was closely related to geographical distance, and temperature and precipitation notablely affected the genetic variation and genetic differentiation of C. tientaiensis. Based on cpDNA, this study indicated that C. tientaiensis exhibits a moderate level of genetic diversity, and which mainly came from among populations, while Tiantai Mountain population have the highest genetic diversity. It demonstrated that there was genetic differentiation between populations, which can be divided into two independent geographical groups, but there was no significant phylogeographic structure between them. The MaxEnt model showed that climate change significantly affects its distribution, and the suitable distribution areas in Zhejiang were primarily divided into two regions, eastern Zhejiang and southern Zhejiang, and there was niche differentiation in its suitable distribution areas. Therefore, this study speculated that the climate and the terrain of mountains and hills in East China jointly shape the genetic structure of C. tientaiensis, which gived rise to an obvious north-south differentiation trend of these species, and the populations located in the hilly areas of eastern Zhejiang and the mountainous areas of southern Zhejiang formed two genetic branches respectively.

Heritable epigenetic variation facilitates long-term maintenance of epigenetic and genetic variation

How genetic and phenotypic variation are maintained has long been one of the fundamental questions in population and quantitative genetics. A variety of factors have been implicated to explain the maintenance of genetic variation in some contexts (e.g. balancing selection), but the potential role of epigenetic regulation to influence population dynamics has been understudied. It is well recognized that epigenetic regulation, including histone methylation, small RNA expression, and DNA methylation, helps to define differences between cell types and facilitate phenotypic plasticity. In recent years, empirical studies have shown the potential for epigenetic regulation to also be heritable for at least a few generations without selection, raising the possibility that differences in epigenetic regulation can act alongside genetic variation to shape evolutionary trajectories. Heritable differences in epigenetic regulation that arise spontaneously are termed "epimutations." Epimutations differ from genetic mutations in 2 key ways-they occur at a higher rate and the loci at which they occur often revert back to their original state within a few generations. Here, we present an extension of the standard population genetic model with selection to incorporate epigenetic variation arising via epimutation. Our model assumes a diploid, sexually reproducing population with random mating. In addition to spontaneous genetic mutation, we included parameters for spontaneous epimutation and back-epimutation, allowing for 4 potential epialleles at a single locus (2 genetic alleles, each with 2 epigenetic states), each of which affect fitness. We then analyzed the conditions under which stable epialleles were maintained. Our results show that highly reversible epialleles can be maintained in long-term equilibrium under neutral conditions in a manner that depends on the epimutation and back-epimutation rates, which we term epimutation-back-epimutation equilibrium. On the other hand, epialleles that compensate for deleterious mutations cause deviations from the expectations of mutation-selection balance by a simple factor that depends on the epimutation and back-epimutation rates. We also numerically analyze several sets of fitness parameters for which large deviations from mutation-selection balance occur. Together, these results demonstrate that transient epigenetic regulation may be an important factor in the maintenance of both epigenetic and genetic variation in populations.

Comparative Analysis of Chloroplast Pan-Genomes and Transcriptomics Reveals Cold Adaptation in Medicago sativa

Alfalfa (Medicago sativa) is a perennial forage legume that is widely distributed all over the world; therefore, it has an extremely complex genetic background. Though population structure and phylogenetic studies have been conducted on a large group of alfalfa nuclear genomes, information about the chloroplast genomes is still lacking. Chloroplast genomes are generally considered to be conservative and play an important role in population diversity analysis and species adaptation in plants. Here, 231 complete alfalfa chloroplast genomes were successfully assembled from 359 alfalfa resequencing data, on the basis of which the alfalfa chloroplast pan-genome was constructed. We investigated the genetic variations of the alfalfa chloroplast genome through comparative genomic, genetic diversity, phylogenetic, population genetic structure, and haplotype analysis. Meanwhile, the expression of alfalfa chloroplast genes under cold stress was explored through transcriptome analysis. As a result, chloroplast genomes of 231 alfalfa lack an IR region, and the size of the chloroplast genome ranges from 125,192 bp to 126,105 bp. Using population structure, haplotypes, and construction of a phylogenetic tree, it was found that alfalfa populations could be divided into four groups, and multiple highly variable regions were found in the alfalfa chloroplast genome. Transcriptome analysis showed that tRNA genes were significantly up-regulated in the cold-sensitive varieties, while rps7, rpl32, and ndhB were down-regulated, and the editing efficiency of ycf1, ycf2, and ndhF was decreased in the cold-tolerant varieties, which may be due to the fact that chloroplasts store nutrients through photosynthesis to resist cold. The huge number of genetic variants in this study provide powerful resources for molecular markers.

Variation in phenotypic plasticity across age-at-maturity genotypes in wild Atlantic salmon

Evolution of phenotypic plasticity requires genotype-environment interaction. The discovery of two large-effect loci in the vgll3 and six6 genomic regions associated with the number of years the Atlantic salmon spend feeding at sea before maturation (sea age), provides a unique opportunity to study evolutionary potential of phenotypic plasticity. Using data on 1246 Atlantic salmon caught in the River Surna in Norway, we show that variation in mean sea age among years (smolt cohorts 2013-2018) is influenced by genotype frequencies as well as interaction effects between genotype and year. Genotype-year interactions suggest that genotypes may differ in their response to environmental variation across years, implying genetic variation in phenotypic plasticity. Our results also imply that plasticity in sea age will evolve as an indirect response to selection on mean sea age due to a shared genetic basis. Furthermore, we demonstrate differences between years in the additive and dominance functional genetic effects of vgll3 and six6 on sea age, suggesting that evolutionary responses will vary across environments. Considering the importance of age at maturity for survival and reproduction, genotype-environment interactions likely play an important role in local adaptation and population demography in Atlantic salmon.

Whole genome sequencing of 4,787 individuals identifies gene-based rare variants in age-related macular degeneration

Genome-wide association studies have contributed extensively to the discovery of disease-associated common variants. However, the genetic contribution to complex traits is still largely difficult to interpret. We report a genome-wide association study of 2394 cases and 2393 controls for age-related macular degeneration (AMD) via whole-genome sequencing, with 46.9 million genetic variants. Our study reveals significant single-variant association signals at four loci and independent gene-based signals in CFH, C2, C3, and NRTN. Using data from the Exome Aggregation Consortium (ExAC) for a gene-based test, we demonstrate an enrichment of predicted rare loss-of-function variants in CFH, CFI, and an as-yet unreported gene in AMD, ORMDL2. Our method of using a large variant list without individual-level genotypes as an external reference provides a flexible and convenient approach to leverage the publicly available variant datasets to augment the search for rare variant associations, which can explain additional disease risk in AMD.

Disentangling the contemporary and historical effects of landscape on the population genomic variation of two bird species restricted to the highland forest enclaves of northeastern Brazil

Investigating the impact of landscape features on patterns of genetic variation is crucial to understand spatially dependent evolutionary processes. Here, we assess the population genomic variation of two bird species (Conopophaga cearae and Sclerurus cearensis) through the Caatinga moist forest enclaves in northeastern Brazil. To infer the evolutionary dynamics of bird populations through the Late Quaternary, we used genome-wide polymorphism data obtained from double-digestion restriction-site-associated DNA sequencing (ddRADseq), and integrated population structure analyses, historical demography models, paleodistribution modeling, and landscape genetics analyses. We found the population differentiation among enclaves to be significantly related to the geographic distance and historical resistance across the rugged landscape. The climate changes at the end of the Pleistocene to the Holocene likely triggered synchronic population decline in all enclaves for both species. Our findings revealed that both geographic distance and historical connectivity through highlands are important factors that can explain the current patterns of genetic variation. Our results further suggest that levels of population differentiation and connectivity cannot be explained purely on the basis of contemporary environmental conditions. By combining historical demographic analyses and niche modeling predictions in a historical framework, we provide strong evidence that climate fluctuations of the Quaternary promoted population differentiation and a high degree of temporal synchrony among population size changes in both species.

Monitoring of species' genetic diversity in Europe varies greatly and overlooks potential climate change impacts

Genetic monitoring of populations currently attracts interest in the context of the Convention on Biological Diversity but needs long-term planning and investments. However, genetic diversity has been largely neglected in biodiversity monitoring, and when addressed, it is treated separately, detached from other conservation issues, such as habitat alteration due to climate change. We report an accounting of efforts to monitor population genetic diversity in Europe (genetic monitoring effort, GME), the evaluation of which can help guide future capacity building and collaboration towards areas most in need of expanded monitoring. Overlaying GME with areas where the ranges of selected species of conservation interest approach current and future climate niche limits helps identify whether GME coincides with anticipated climate change effects on biodiversity. Our analysis suggests that country area, financial resources and conservation policy influence GME, high values of which only partially match species' joint patterns of limits to suitable climatic conditions. Populations at trailing climatic niche margins probably hold genetic diversity that is important for adaptation to changing climate. Our results illuminate the need in Europe for expanded investment in genetic monitoring across climate gradients occupied by focal species, a need arguably greatest in southeastern European countries. This need could be met in part by expanding the European Union's Birds and Habitats Directives to fully address the conservation and monitoring of genetic diversity.

Evolution of chemosensory genes in Colorado potato beetle, Leptinotarsa decemlineata

Associating with plant hosts is thought to have elevated the diversification of insect herbivores, which comprise the majority of global species diversity. In particular, there is considerable interest in understanding the genetic changes that allow host-plant shifts to occur in pest insects and in determining what aspects of functional genomic diversity impact host-plant breadth. Insect chemoreceptors play a central role in mediating insect-plant interactions, as they directly influence plant detection and sensory stimuli during feeding. Although chemosensory genes evolve rapidly, it is unclear how they evolve in response to host shifts and host specialization. We investigate whether selection at chemosensory genes is linked to host-plant expansion from the buffalo burr, Solanum rostratum, to potato, Solanum tuberosum, in the super-pest Colorado potato beetle (CPB), Leptinotarsa decemlineata (Coleoptera: Chrysomelidae). First, to refine our knowledge of CPB chemosensory genes, we developed novel gene expression data for the antennae and maxillary-labial palps. We then examine patterns of selection at these loci within CPB, as well as compare whether rates of selection vary with respect to 9 closely related, non-pest Leptinotarsa species that vary in diet breadth. We find that rates of positive selection on olfactory receptors are higher in host-plant generalists, and this signal is particularly strong in CPB. These results provide strong candidates for further research on the genetic basis of variation in insect chemosensory performance and novel targets for pest control of a notorious super-pest.

A chromosome-scale assembly reveals chromosomal aberrations and exchanges generating genetic diversity in Coffea arabica germplasm

In order to better understand the mechanisms generating genetic diversity in the recent allotetraploid species Coffea arabica, here we present a chromosome-level assembly obtained with long read technology. Two genomic compartments with different structural and functional properties are identified in the two homoeologous genomes. The resequencing data from a large set of accessions reveals low intraspecific diversity in the center of origin of the species. Across a limited number of genomic regions, diversity increases in some cultivated genotypes to levels similar to those observed within one of the progenitor species, Coffea canephora, presumably as a consequence of introgressions deriving from the so-called Timor hybrid. It also reveals that, in addition to few, early-occurring exchanges between homoeologous chromosomes, there are numerous recent chromosomal aberrations including aneuploidies, deletions, duplications and exchanges. These events are still polymorphic in the germplasm and could represent a fundamental source of genetic variation in such a lowly variable species.