The variation and evolution of complete human centromeres

Human centromeres have been traditionally very difficult to sequence and assemble owing to their repetitive nature and large size1. As a result, patterns of human centromeric variation and models for their evolution and function remain incomplete, despite centromeres being among the most rapidly mutating regions2,3. Here, using long-read sequencing, we completely sequenced and assembled all centromeres from a second human genome and compared it to the finished reference genome4,5. We find that the two sets of centromeres show at least a 4.1-fold increase in single-nucleotide variation when compared with their unique flanks and vary up to 3-fold in size. Moreover, we find that 45.8% of centromeric sequence cannot be reliably aligned using standard methods owing to the emergence of new α-satellite higher-order repeats (HORs). DNA methylation and CENP-A chromatin immunoprecipitation experiments show that 26% of the centromeres differ in their kinetochore position by >500 kb. To understand evolutionary change, we selected six chromosomes and sequenced and assembled 31 orthologous centromeres from the common chimpanzee, orangutan and macaque genomes. Comparative analyses reveal a nearly complete turnover of α-satellite HORs, with characteristic idiosyncratic changes in α-satellite HORs for each species. Phylogenetic reconstruction of human haplotypes supports limited to no recombination between the short (p) and long (q) arms across centromeres and reveals that novel α-satellite HORs share a monophyletic origin, providing a strategy to estimate the rate of saltatory amplification and mutation of human centromeric DNA.

Demographic fluctuations and selection during host-parasite co-evolution interactively increase genetic diversity

Host-parasite interactions can cause strong demographic fluctuations accompanied by selective sweeps of resistance/infectivity alleles. Both demographic bottlenecks and frequent sweeps are expected to reduce the amount of segregating genetic variation and therefore might constrain adaptation during co-evolution. Recent studies, however, suggest that the interaction of demographic and selective processes is a key component of co-evolutionary dynamics and may rather positively affect levels of genetic diversity available for adaptation. Here, we provide direct experimental testing of this hypothesis by disentangling the effects of demography, selection and their interaction in an experimental host-parasite system. We grew 12 populations of a unicellular, asexually reproducing algae (Chlorella variabilis) that experienced either growth followed by constant population sizes (three populations), demographic fluctuations (three populations), selection induced by exposure to a virus (three populations), or demographic fluctuations together with virus-induced selection (three populations). After 50 days (~50 generations), we conducted whole-genome sequencing of each algal host population. We observed more genetic diversity in populations that jointly experienced selection and demographic fluctuations than in populations where these processes were experimentally separated. In addition, in those three populations that jointly experienced selection and demographic fluctuations, experimentally measured diversity exceeds expected values of diversity that account for the cultures' population sizes. Our results suggest that eco-evolutionary feedbacks can positively affect genetic diversity and provide the necessary empirical measures to guide further improvements of theoretical models of adaptation during host-parasite co-evolution.

New insight into genetic diversity of zoonotic-potential Ancylostoma ceylanicum in stray cats living in Bangkok, Thailand, based on deep amplicon sequencing

AIMS

This study aimed to characterize feline hookworms from stray cats living in Bangkok.

METHODS AND RESULTS

A total of 56 hookworm-positive faecal samples were identified for hookworm species by using PCR targeting the ITS1, 5.8S, and ITS2 fragment and qPCR targeting ITS2. Of 56 samples, 96.4% (54/56) were identified as Ancylostoma ceylanicum and 1.8% (1/56) as Ancylostoma caninum. With qPCR, 89.3% (50/56) were identified as single A. ceylanicum infection and 5.4% (3/56) as coinfection of A. ceylanicum and A. caninum. For genetic characterization of A. ceylanicum, 10 samples were pooled, and the partial COI gene was amplified, followed by deep amplicon sequencing. Five pooled samples were analysed, and 99.73% were identified with A. ceylanicum sequences, which were allocated into 19 haplotypes (AC01-AC19). Genetic diversity findings for A. ceylanicum in Asia revealed that three of eight haplotypes considered of zoonotic significance occurred in humans, dogs, and cats, including haplotypes H01, H20, and H21. The predominant haplotype in this study, AC01, was clustered with H01-a zoonotic haplotype.

CONCLUSIONS

The diversity obtained by deep amplicon sequencing supported that the A. ceylanicum community had high genetic variation. Deep amplicon sequencing was a useful method to determine source, zoonotic potential, and host-parasite relationship.

Evidence for the additivity of rare and common variant burden throughout the spectrum of intellectual disability

Intellectual disability (ID) is a common disorder, yet there is a wide spectrum of impairment from mild to profoundly affected individuals. Mild ID is seen as the low extreme of the general distribution of intelligence, while severe ID is often seen as a monogenic disorder caused by rare, pathogenic, highly penetrant variants. To investigate the genetic factors influencing mild and severe ID, we evaluated rare and common variation in the Northern Finland Intellectual Disability cohort (n = 1096 ID patients), a cohort with a high percentage of mild ID (n = 550) and from a population bottleneck enriched in rare, damaging variation. Despite this enrichment, we found only a small percentage of ID was due to recessive Finnish-enriched variants (0.5%). A larger proportion was linked to dominant variation, with a significant burden of rare, damaging variation in both mild and severe ID. This rare variant burden was enriched in more severe ID (p = 2.4e-4), patients without a relative with ID (p = 4.76e-4), and in those with features associated with monogenic disorders. We also found a significant burden of common variants associated with decreased cognitive function, with no difference between mild and more severe ID. When we included common and rare variants in a joint model, the rare and common variants had additive effects in both mild and severe ID. A multimodel inference approach also found that common and rare variants together best explained ID status (ΔAIC = 16.8, ΔBIC = 10.2). Overall, we report evidence for the additivity of rare and common variant burden throughout the spectrum of intellectual disability.

Genetic characteristics of Blastocystis sp. in cattle from Hebei Province, China

Blastocystis sp. is a protozoan parasite that infects the intestines of humans and animals, causing chronic diseases such as skin rashes, abdominal pain, and irritable bowel syndrome. A survey was conducted to determine the prevalence and genetic diversity of Blastocystis sp. infection in cattle, in Hebei Province, China. 2746 cattle fecal samples were collected from 11 cities in Hebei Province and analyzed using polymerase chain reaction targeting the Blastocystis sp. barcoding gene. MEGA, PhyloSuite, and PopART were used to analyze the subtype, sequence signature, pairwise genetic distance, and genetic diversity indices. The results showed that the Blastocystis sp. detection rate was 12.60% (346/2746). The infection rate in different herds was affected by region, age, breeding mode, and variety; that is, the infection rates in areas of southern Hebei, cattle under one year old, intensive raising, and dairy cattle were higher than the infection rates in northern Hebei, cattle over one year old, scatter feeding, and beef cattle. Seven Blastocystis subtypes were identified, namely, ST1, ST2, ST5, ST10, ST14, ST21, and ST26; ST10 was the dominant subtype, and ST14 was the second most common subtype. A total of 374 polymorphic and conserved sites were obtained, including 273 invariable (monomorphic) sites and 101 variable (polymorphic) sites, accounting for 27.01% of all nucleotides. The nucleotide diversity index (Pi) was 0.07749, and the haplotype (gene) diversity index (Hd) was 0.946. This study provides the first comprehensive information on the epidemiological situation of Blastocystis sp. infection in cattle from Hebei Province, China, and revealed rich genetic diversity of Blastocystis sp.

Subtypes of Blastocystis in Tibetan Antelope (Pantholops hodgsonii)

Blastocystis is a protist that is distributed in the gut tract of humans and animals. However, the reports about Blastocystis infection in Tibetan antelope are scarce. We collected 173 Tibetan antelope feces samples from Xinjiang, Qinghai and Xizang, and amplified the SSU rRNA gene of 600 bp region of Blastocystis in our research. Fifty-one samples in total were positive for Blastocystis, with all subtypes being ST31. The lowest prevalence of Blastocystis was observed in Xizang (2/20, 9.1%), followed by Qinghai (18/92, 16.4%), Xinjiang (31/61, 33.7%). The highest prevalence of Blastocystis in Tibetan antelope was detected during the summer was (19/30, 38.8%). This is the first research work regarding the Blastocystis subtypes ST31 in Tibetan antelope. Our research provides information for future researches on the distribution of this Blastocystis subtype and the control of Blastocystis infection.

Genetically transitional disease: conceptual understanding and applicability to rheumatic disease

In genomic medicine, the concept of genetically transitional disease (GTD) refers to cases in which gene mutation is necessary but not sufficient to cause disease. In this Perspective, we apply this novel concept to rheumatic diseases, which have been linked to hundreds of genetic variants via association studies. These variants are in the 'grey zone' between monogenic variants with large effect sizes and common susceptibility alleles with small effect sizes. Among genes associated with rare autoinflammatory diseases, many low-frequency and/or low-penetrance variants are known to increase susceptibility to systemic inflammation. In autoimmune diseases, hundreds of HLA and non-HLA genetic variants have been revealed to be modest- to moderate-risk alleles. These diseases can be reclassified as GTDs. The same concept could apply to many other human diseases. GTD could improve the reporting of genetic testing results, diagnostic yields, genetic counselling and selection of therapy, as well as facilitating research using a novel approach to human genetic diseases.

Chromosome-level genome assembly and population genomic analysis reveal evolution and local adaptation in common hairfin anchovy (Setipinna tenuifilis)

Understanding the genetic structure and the factors associated with adaptive diversity has significant implications for the effective management of wild populations under threat from overfishing and climate change. The common hairfin anchovy (Setipinna tenuifilis) is an economically and ecologically important pelagic fish species, spanning a broad latitudinal gradient along marginal seas of the Northwest Pacific. In this study, we constructed the first reference genome of S. tenuifilis using PacBio long reads and high-resolution chromosome conformation capture (Hi-C) technology. The assembled genome was 798.38 Mb with a contig N50 of 1.43 Mb and a scaffold N50 of 32.42 Mb, which were anchored onto 24 pseudochromosomes. A total of 22,019 genes were functionally annotated, which accounted for 95.27% of the predicted protein-coding genes. Chromosomal collinearity analysis revealed chromosome fusion or fission events in Clupeiformes species. Three genetic groups of S. tenuifilis were revealed along the Chinese coast using restriction site-associated DNA sequencing (RADseq). We investigated the influence of four bioclimatic variables as potential drivers of adaptive divergence in S. tenuifilis, suggesting that these environmental variables, especially sea surface temperature, may play important roles as drivers of spatially varying selection for S. tenuifilis. We also identified candidate functional genes underlying adaptive mechanisms and ecological tradeoffs using redundancy analysis (RDA) and BayeScan analysis. In summary, this study sheds light on the evolution and spatial patterns of genetic variation of S. tenuifilis, providing a valuable genomic resource for further biological and genetic studies on this species and other closely related Clupeiformes.

Genetic variation in parasite avoidance, yet no evidence for constitutive fitness costs

Behavioral avoidance of parasites is a widespread strategy among animal hosts and in human public health. Avoidance has repercussions for both individual and population-level infection risk. Although most cases of parasite avoidance are viewed as adaptive, there is little evidence that the basic assumptions of evolution by natural selection are met. This study addresses this gap by testing whether there is a heritable variation in parasite avoidance behavior. We quantified behavioral avoidance of the bacterial parasite Serratia marcescens for 12 strains of the nematode host Caenorhabditis elegans. We found that these strains varied in their magnitude of avoidance, and we estimated the broad-sense heritability of this behavior to be in the range of 11%-26%. We then asked whether avoidance carries a constitutive fitness cost. We did not find evidence of one. Rather, strains with higher avoidance had higher fitness, measured as population growth rate. Together, these results direct future theoretical and empirical work to identify the forces maintaining genetic variation in parasite avoidance.

Molnupiravir increases SARS-CoV-2 genome diversity and complexity: A case-control cohort study

Molnupiravir, an oral direct-acting antiviral effective in vitro against SARS-CoV-2, has been largely employed during the COVID-19 pandemic, since December 2021. After marketing and widespread usage, a progressive increase in SARS-CoV-2 lineages characterized by a higher transition/transversion ratio, a characteristic signature of molnupiravir action, appeared in the Global Initiative on Sharing All Influenza Data (GISAID) and International Nucleotide Sequence Database Collaboration (INSDC) databases. Here, we assessed the drug effects by SARS-CoV-2 whole-genome sequencing on 38 molnupiravir-treated persistently positive COVID-19 outpatients tested before and after treatment. Seventeen tixagevimab/cilgavimab-treated outpatients served as controls. Mutational analyses confirmed that SARS-CoV-2 exhibits an increased transition/transversion ratio seven days after initiation of molnupiravir. Moreover we observed an increased G->A ratio compared to controls, which was not related to apolipoprotein B mRNAediting enzyme, catalytic polypeptide-like (APOBEC) activity. In addition, we demonstrated for the first time an increased diversity and complexity of the viral quasispecies.

Rapid adaptation of recombining populations on tunable fitness landscapes

How does standing genetic variation affect polygenic adaptation in recombining populations? Despite a large body of work in quantitative genetics, epistatic and weak additive fitness effects among simultaneously segregating genetic variants are difficult to capture experimentally or to predict theoretically. In this study, we simulated adaptation on fitness landscapes with tunable ruggedness driven by standing genetic variation in recombining populations. We confirmed that recombination hinders the movement of a population through a rugged fitness landscape. When surveying the effect of epistasis on the fixation of alleles, we found that the combined effects of high ruggedness and high recombination probabilities lead to preferential fixation of alleles that had a high initial frequency. This indicates that positive epistatic alleles escape from being broken down by recombination when they start at high frequency. We further extract direct selection coefficients and pairwise epistasis along the adaptive path. When taking the final fixed genotype as the reference genetic background, we observe that, along the adaptive path, beneficial direct selection appears stronger and pairwise epistasis weaker than in the underlying fitness landscape. Quantitatively, the ratio of epistasis and direct selection is smaller along the adaptive path (≈ 1 ) than expected. Thus, adaptation on a rugged fitness landscape may lead to spurious signals of direct selection generated through epistasis. Our study highlights how the interplay of epistasis and recombination constrains the adaptation of a diverse population to a new environment.

Low-frequency and rare genetic variants associated with rheumatoid arthritis risk

Rheumatoid arthritis (RA) has an estimated heritability of nearly 50%, which is particularly high in seropositive RA. HLA alleles account for a large proportion of this heritability, in addition to many common single-nucleotide polymorphisms with smaller individual effects. Low-frequency and rare variants, such as those captured by next-generation sequencing, can also have a large role in heritability in some individuals. Rare variant discovery has informed the development of drugs such as inhibitors of PCSK9 and Janus kinases. Some 34 low-frequency and rare variants are currently associated with RA risk. One variant (19:10352442G>C in TYK2) was identified in five separate studies, and might therefore represent a promising therapeutic target. Following a set of best practices in future studies, including studying diverse populations, using large sample sizes, validating RA and serostatus, replicating findings, adjusting for other variants and performing functional assessment, could help to ensure the relevance of identified variants. Exciting opportunities are now on the horizon for genetics in RA, including larger datasets and consortia, whole-genome sequencing and direct applications of findings in the management, and especially treatment, of RA.

Investigation of Chlamydia pecorum in livestock from Switzerland reveals a high degree of diversity in bovine strains

Chlamydia pecorum is a widespread veterinary chlamydial species causing endemic infections in livestock, such as ruminants and pigs, globally. However, there is limited contemporary knowledge on infecting strain diversity in various hosts. This study aimed to evaluate the genetic diversity of C. pecorum strains infecting Swiss livestock through C. pecorum genotyping and phylogenetic analyses in comparison to the global population, while also assessing chlamydial strains for plasmid carriage. A total of 263 C. pecorum positive samples from clinically healthy ruminant and pig herds (Bovines = 216, sheep = 25, pigs = 14) as well as placentae from eight C. pecorum positive ruminant abortion cases from other Swiss herds were investigated. The ompA and Multi-Locus sequence typing revealed novel C. pecorum genotypes, and bovine strains exhibited considerable genetic diversity, contrasting with lower diversity in sheep and pig strains. C. pecorum plasmid was detected in 100.0% of sheep (41/41) and pig (255/255) samples, and in 69.4% of bovine samples (150/216). In contrast, no plasmid was detected in the eight C. pecorum-positive ruminant abortion cases either representing plasmid-less strains or possibly escaping PCR detection due to autolysis of the placenta. This study supports the genetic diversity of C. pecorum strains, particularly in bovines, and identifies novel sequence types in Swiss livestock.

Toward breeding pigmented rice balancing nutrition and yield

A recent study by Sedeek et al. provides multiomic resources that illustrate the genetic diversity, metabolites, elemental composition, and the possibility of agronomic trait improvement, through genome-editing technology, for nutrient-rich pigmented rice. This will guide future rice breeding programs for balancing optimal agronomic traits and excellent nutritional quality.

Genomic signatures of strawberry domestication and diversification

Cultivated strawberry (Fragaria × ananassa) has a brief history of less than 300 yr, beginning with the hybridization of octoploids Fragaria chiloensis and Fragaria virginiana. Here we explored the genomic signatures of early domestication and subsequent diversification for different climates using whole-genome sequences of 289 wild, heirloom, and modern varieties from two major breeding programs in the United States. Four nonadmixed wild octoploid populations were identified, with recurrent introgression among the sympatric populations. The proportion of F. virginiana ancestry increased by 20% in modern varieties over initial hybrids, and the proportion of F. chiloensis subsp. pacifica rose from 0% to 3.4%. Effective population size rapidly declined during early breeding. Meanwhile, divergent selection for distinct environments reshaped wild allelic origins in 21 out of 28 chromosomes. Overlapping divergent selective sweeps in natural and domesticated populations revealed 16 convergent genomic signatures that may be important for climatic adaptation. Despite 20 breeding cycles since initial hybridization, more than half of loci underlying yield and fruit size are still not under artificial selection. These insights add clarity to the domestication and breeding history of what is now the most widely cultivated fruit in the world.

Genomic evolution of human respiratory syncytial virus during a decade (2013-2023): bridging the path to monoclonal antibody surveillance

OBJECTIVES

This study investigated the prevalence, genetic diversity, and evolution of human respiratory syncytial virus (HRSV) in Barcelona from 2013 to 2023.

METHODS

Respiratory specimens from patients with RTI suspicion at Hospital Universitari Vall d'Hebron were collected from October 2013 to May 2023 for laboratory-confirmation of respiratory viruses. Next-generation sequencing was performed in randomly-selected samples with Illumina technology. Phylogenetic analyses of whole genome sequences were performed with BEAST v1.10.4. Signals of selection and evolutionary pressures were inferred by population dynamics and evolutionary analyses. Mutations in major surface proteins were genetic and structurally characterised, emphasizing those within antigenic epitopes.

RESULTS

Analyzing 139,625 samples, 5.3% were HRSV-positive (3008 HRSV-A, 3882 HRSV-B, 56 HRSV-A and -B, and 495 unsubtyped HRSV), with a higher prevalence observed in the paediatric population. Pandemic-related shifts in seasonal patterns returned to normal in 2022-2023. A total of 198 whole-genome sequences were obtained for HRSV-A (6.6% of the HRSV-A positive samples) belonging to GA2.3.5 lineage. For HRSV-B, 167 samples were sequenced (4.3% of the HRSV-B positive samples), belonging to GB5.0.2, GB5.0.4a and GB5.0.5a. HRSV-B exhibited a higher evolution rate. Post-SARS-CoV-2 pandemic, both subtypes showed increased evolutionary rates and decreased effective population size initially, followed by a sharp increase. Analyses indicated negative selective pressure on HRSV. Mutations in antigenic epitopes, including S276N and M274I in palivizumab-targeted site II, and I206M, Q209R, and S211N in nirsevimab-targeted site Ø, were identified.

DISCUSSION

Particularly in the context of the large-scale use in 2023-2024 season of nirsevimab, continuous epidemiological and genomic surveillance is crucial.

Investigation of inherited noncoding genetic variation impacting the pharmacogenomics of childhood acute lymphoblastic leukemia treatment

Defining genetic factors impacting chemotherapy failure can help to better predict response and identify drug resistance mechanisms. However, there is limited understanding of the contribution of inherited noncoding genetic variation on inter-individual differences in chemotherapy response in childhood acute lymphoblastic leukemia (ALL). Here we map inherited noncoding variants associated with treatment outcome and/or chemotherapeutic drug resistance to ALL cis-regulatory elements and investigate their gene regulatory potential and target gene connectivity using massively parallel reporter assays and three-dimensional chromatin looping assays, respectively. We identify 54 variants with transcriptional effects and high-confidence gene connectivity. Additionally, functional interrogation of the top variant, rs1247117, reveals changes in chromatin accessibility, PU.1 binding affinity and gene expression, and deletion of the genomic interval containing rs1247117 sensitizes cells to vincristine. Together, these data demonstrate that noncoding regulatory variants associated with diverse pharmacological traits harbor significant effects on allele-specific transcriptional activity and impact sensitivity to antileukemic agents.

Plastic and genomic change of a newly established lizard population following a founder event

Understanding how phenotypic divergence arises among natural populations remains one of the major goals in evolutionary biology. As part of competitive exclusion experiment conducted in 1971, 10 individuals of Italian wall lizard (Podarcis siculus (Rafinesque-Schmaltz, 1810)) were transplanted from Pod Kopište Island to the nearby island of Pod Mrčaru (Adriatic Sea). Merely 35 years after the introduction, the newly established population on Pod Mrčaru Island had shifted their diet from predominantly insectivorous towards omnivorous and changed significantly in a range of morphological, behavioural, physiological and ecological characteristics. Here, we combine genomic and quantitative genetic approaches to determine the relative roles of genetic adaptation and phenotypic plasticity in driving this rapid phenotypic shift. Our results show genome-wide genetic differentiation between ancestral and transplanted population, with weak genetic erosion on Pod Mrčaru Island. Adaptive processes following the founder event are indicated by highly differentiated genomic loci associating with ecologically relevant phenotypic traits, and/or having a putatively adaptive role across multiple lizard populations. Diverged traits related to head size and shape or bite force showed moderate heritability in a crossing experiment, but between-population differences in these traits did not persist in a common garden environment. Our results confirm the existence of sufficient additive genetic variance for traits to evolve under selection while also demonstrating that phenotypic plasticity and/or genotype by environment interactions are the main drivers of population differentiation at this early evolutionary stage.

A comprehensive review of HVS-I mitochondrial DNA variation of 19 Iranian populations

Iran is located along the Central Asian corridor, a natural artery that has served as a cross-continental route since the first anatomically modern human populations migrated out of Africa. We compiled and reanalyzed the HVS-I (hypervariable segment-I) of 3840 mitochondrial DNA (mtDNA) sequences from 19 Iranian populations and from 26 groups from adjacent countries to give a comprehensive review of the maternal genetic variation and investigate the impact of historical events and cultural factors on the maternal genetic structure of modern Iranians. We conclude that Iranians have a high level of genetic diversity. Thirty-six haplogroups were observed in Iran's populations, and most of them belong to widespread West-Eurasian haplogroups, such as H, HV, J, N, T, and U. In contrast, the predominant haplogroups observed in most of the adjacent countries studied here are H, M, D, R, U, and C haplogroups. Using principal component analysis, clustering, and genetic distance-based calculations, we estimated moderate genetic relationships between Iranian and other Eurasian groups. Further, analyses of molecular variance and comparing geographic and genetic structures indicate that mtDNA HVS-I sequence diversity does not exhibit any sharp geographic structure in the country. Barring a few from some culturally distinct and naturally separated minorities, most Iranian populations have a homogenous maternal genetic structure.

Genetic control of root morphology in response to nitrogen across rapeseed diversity

Rapeseed (Brassica napus L.) is an oil-containing crop of great economic value but with considerable nitrogen requirement. Breeding root systems that efficiently absorb nitrogen from the soil could be a driver to ensure genetic gains for more sustainable rapeseed production. The aim of this study is to identify genomic regions that regulate root morphology in response to nitrate availability. The natural variability offered by 300 inbred lines was screened at two experimental locations. Seedlings grew hydroponically with low or elevated nitrate levels. Fifteen traits related to biomass production and root morphology were measured. On average across the panel, a low nitrate level increased the root-to-shoot biomass ratio and the lateral root length. A large phenotypic variation was observed, along with important heritability values and genotypic effects, but low genotype-by-nitrogen interactions. Genome-wide association study and bulk segregant analysis were used to identify loci regulating phenotypic traits. The first approach nominated 319 SNPs that were combined into 80 QTLs. Three QTLs identified on the A07 and C07 chromosomes were stable across nitrate levels and/or experimental locations. The second approach involved genotyping two groups of individuals from an experimental F2 population created by crossing two accessions with contrasting lateral root lengths. These individuals were found in the tails of the phenotypic distribution. Co-localized QTLs found in both mapping approaches covered a chromosomal region on the A06 chromosome. The QTL regions contained some genes putatively involved in root organogenesis and represent selection targets for redesigning the root morphology of rapeseed.

Genome organization and botanical diversity

The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes-both in content and in dynamics-has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.

The importance of variant reinterpretation in inherited cardiovascular diseases: Establishing the optimal timeframe

Inherited cardiovascular diseases are rare diseases that are difficult to diagnose by non-expert professionals. Genetic analyses play a key role in the diagnosis of these diseases, in which the identification of a pathogenic genetic variant is often a diagnostic criterion. Therefore, genetic variant classification and routine reinterpretation as data become available represent one of the main challenges associated with genetic analyses. Using the genetic variants identified in an inherited cardiovascular diseases unit during a 10-year period, the objectives of this study were: 1) to evaluate the impact of genetic variant reinterpretation, 2) to compare the reclassification rates between different cohorts of cardiac channelopathies and cardiomyopathies, and 3) to establish the most appropriate periodicity for genetic variant reinterpretation. All the evaluated cohorts (full cohort of inherited cardiovascular diseases, cardiomyopathies, cardiac channelopathies, hypertrophic cardiomyopathy, dilated cardiomyopathy, arrhythmogenic cardiomyopathy, Brugada syndrome, long QT syndrome and catecholaminergic polymorphic ventricular tachycardia) showed reclassification rates above 25%, showing even higher reclassification rates when there is definitive evidence of the association between the gene and the disease in the cardiac channelopathies. Evaluation of genetic variant reclassification rates based on the year of the initial classification showed that the most appropriate frequency for the reinterpretation would be 2 years, with the possibility of a more frequent reinterpretation if deemed convenient. To keep genetic variant classifications up to date, genetic counsellors play a critical role in the reinterpretation process, providing clinical evidence that genetic diagnostic laboratories often do not have at their disposal and communicating changes in classification and the potential implications of these reclassifications to patients and relatives.

Pedigree diversity and implications for genetic selection of Katahdin sheep

The Katahdin hair breed gained popularity in the United States as low input and prolific, with a propensity to exhibit parasite resistance. With the introduction of genomically enhanced estimated breeding values (GEBV) to the Katahdin genetic evaluation, defining the diversity present in the breed is pertinent. Utilizing pedigree records (n = 92,030) from 1984 to 2019 from the National Sheep Improvement Program, our objectives were to (i) estimate the completeness and quality of the pedigree, (ii) calculate diversity statistics for the whole pedigree and relevant reference subpopulations and (iii) assess the impact of current diversity on genomic selection. Reference 1 was Katahdins born from 2017 to 2019 (n = 23,494), while reference 2 was a subset with at least three generations of Katahdin ancestry (n = 9327). The completeness of the whole pedigree, and the pedigrees of reference 1 and reference 2, were above 50% through the fourth, fifth and seventh generation of ancestors, respectively. Effective population size (Ne) averaged 111 animals with a range from 42.2 to 451.0. The average generation interval was 2.9 years for the whole pedigree and reference 1, and 2.8 years for reference 2. The mean individual inbreeding and average relatedness coefficients were 1.62% and 0.91%, 1.74% and 0.90% and 2.94% and 1.46% for the whole pedigree, reference 1, and reference 2, respectively. There were over 300 effective founders in the whole pedigree and reference 1, with 169 in reference 2. Effective number of ancestors were over 150 for the whole pedigree and reference 1, while there were 67 for reference 2. Prediction accuracies increased as the reference population grew from 1k to 7.5k and plateaued at 15k animals. Given the large number of founders and ancestors contributing to the base genetic variation in the breed, the Ne is sufficient to maintain diversity while achieving progress with selection. Stable low rates of inbreeding and relatedness suggest that incorporating genetic conservation in breeding decisions is currently not of high priority. Current Ne suggests that with limited genotyping, high levels of accuracy for genomic prediction can be achieved. However, intense selection on GEBV may cause loss of genetic diversity long term.

Investigation of CaMV-host co-evolution through synonymous codon pattern

Cauliflower mosaic virus (CaMV) has a double-stranded DNA genome and is globally distributed. The phylogeny tree of 121 CaMV isolates was categorized into two primary groups, with Iranian isolates showing the greatest genetic variations. Nucleotide A demonstrated the highest percentage (36.95%) in the CaMV genome and the dinucleotide odds ratio analysis revealed that TC dinucleotide (1.34 ≥ 1.23) and CG dinucleotide (0.63 ≤ 0.78) are overrepresented and underrepresented, respectively. Relative synonymous codon usage (RSCU) analysis confirmed codon usage bias in CaMV and its hosts. Brassica oleracea and Brassica rapa, among the susceptible hosts of CaMV, showed a codon adaptation index (CAI) value above 0.8. Additionally, relative codon deoptimization index (RCDI) results exhibited the highest degree of deoptimization in Raphanus sativus. These findings suggest that the genes of CaMV underwent codon adaptation with its hosts. Among the CaMV open reading frames (ORFs), genes that produce reverse transcriptase and virus coat proteins showed the highest CAI value of 0.83. These genes are crucial for the creation of new virion particles. The results confirm that CaMV co-evolved with its host to ensure the optimal expression of its genes in the hosts, allowing for easy infection and effective spread. To detect the force behind codon usage bias, an effective number of codons (ENC)-plot and neutrality plot were conducted. The results indicated that natural selection is the primary factor influencing CaMV codon usage bias.

Changes in intrahost genetic diversity according to lesion severity in longitudinal HPV16 samples

Human papillomavirus type 16 (HPV16) is the most common cause of cervical cancer, but most infections are transient with lesions not progressing to cancer. There is a lack of specific biomarkers for early cancer risk stratification. This study aimed to explore the intrahost HPV16 genomic variation in longitudinal samples from HPV16-infected women with different cervical lesion severity (normal, low-grade, and high-grade). The TaME-seq deep sequencing protocol was used to generate whole genome HPV16 sequences of 102 samples collected over time from 40 individuals. Single nucleotide variants (SNVs) and intrahost SNVs (iSNVs) were identified in the viral genomes. A majority of individuals had a unique set of SNVs and these SNVs were stable over time. Overall, the number of iSNVs and APOBEC3-induced iSNVs were significantly lower in high-grade relative to normal and low-grade samples. A significant increase in the number of APOBEC3-induced iSNVs over time was observed for normal samples when compared to high-grade. Our results indicates that the lower incidence of iSNVs and APOBEC3-induced iSNVs in high-grade lesions may have implications for novel biomarkers discoveries, potentially aiding early stratification of HPV-induced cervical precancerous lesions.

Genetic variations of penicillin-binding protein 1A: insights into the current status of amoxicillin-based regimens for Helicobacter pylori eradication in Malaysia

Introduction. Resistance towards amoxicillin in Helicobacter pylori causes significant therapeutic impasse in healthcare settings worldwide. In Malaysia, the standard H. pylori treatment regimen includes a 14-day course of high-dose proton-pump inhibitor (rabeprazole, 20 mg) with amoxicillin (1000 mg) dual therapy.Hypothesis/Gap Statement. The high eradication rate with amoxicillin-based treatment could be attributed to the primary resistance rates of amoxicillin being relatively low at 0%, however, a low rate of secondary resistance has been documented in Malaysia recently.Aim. This study aims to investigate the amino acid mutations and related genetic variants in PBP1A of H. pylori, correlating with amoxicillin resistance in the Malaysian population.Methodology. The full-length pbp1A gene was amplified via PCR from 50 genomic DNA extracted from gastric biopsy samples of H. pylori-positive treatment-naïve Malaysian patients. The sequences were then compared with reference H. pylori strain ATCC 26695 for mutation and variant detection. A phylogenetic analysis of 50 sequences along with 43 additional sequences from the NCBI database was performed. These additional sequences included both amoxicillin-resistant strains (n=20) and amoxicillin-sensitive strains (n=23).Results. There was a total of 21 variants of amino acids, with three of them located in or near the PBP-motif (SKN402-404). The percentages of these three variants are as follows: K403X, 2%; S405I, 2% and E406K, 16%. Based on the genetic markers identified, the resistance rate for amoxicillin in our sample remained at 0%. The phylogenetic examination suggested that H. pylori might exhibit unique conserved pbp1A sequences within the Malaysian context.Conclusions. Overall, the molecular analysis of PBP1A supported the therapeutic superiority of amoxicillin-based regimens. Therefore, it is crucial to continue monitoring the amoxicillin resistance background of H. pylori with a larger sample size to ensure the sustained effectiveness of amoxicillin-based treatments in Malaysia.

Aspergillus flavus pangenome (AflaPan) uncovers novel aflatoxin and secondary metabolite associated gene clusters

BACKGROUND

Aspergillus flavus is an important agricultural and food safety threat due to its production of carcinogenic aflatoxins. It has high level of genetic diversity that is adapted to various environments. Recently, we reported two reference genomes of A. flavus isolates, AF13 (MAT1-2 and highly aflatoxigenic isolate) and NRRL3357 (MAT1-1 and moderate aflatoxin producer). Where, an insertion of 310 kb in AF13 included an aflatoxin producing gene bZIP transcription factor, named atfC. Observations of significant genomic variants between these isolates of contrasting phenotypes prompted an investigation into variation among other agricultural isolates of A. flavus with the goal of discovering novel genes potentially associated with aflatoxin production regulation. Present study was designed with three main objectives: (1) collection of large number of A. flavus isolates from diverse sources including maize plants and field soils; (2) whole genome sequencing of collected isolates and development of a pangenome; and (3) pangenome-wide association study (Pan-GWAS) to identify novel secondary metabolite cluster genes.

RESULTS

Pangenome analysis of 346 A. flavus isolates identified a total of 17,855 unique orthologous gene clusters, with mere 41% (7,315) core genes and 59% (10,540) accessory genes indicating accumulation of high genomic diversity during domestication. 5,994 orthologous gene clusters in accessory genome not annotated in either the A. flavus AF13 or NRRL3357 reference genomes. Pan-genome wide association analysis of the genomic variations identified 391 significant associated pan-genes associated with aflatoxin production. Interestingly, most of the significantly associated pan-genes (94%; 369 associations) belonged to accessory genome indicating that genome expansion has resulted in the incorporation of new genes associated with aflatoxin and other secondary metabolites.

CONCLUSION

In summary, this study provides complete pangenome framework for the species of Aspergillus flavus along with associated genes for pathogen survival and aflatoxin production. The large accessory genome indicated large genome diversity in the species A. flavus, however AflaPan is a closed pangenome represents optimum diversity of species A. flavus. Most importantly, the newly identified aflatoxin producing gene clusters will be a new source for seeking aflatoxin mitigation strategies and needs new attention in research.

High clonal diversity of Staphylococcus aureus isolates from children's playgrounds in Hungary

Staphylococcus aureus is one of the most important human pathogenic bacteria and environmental surfaces play an important role in the spread of the bacterium. Presence of S. aureus on children's playgrounds and on toys was described in international studies, however, little is known about the prevalence and characteristics of S. aureus at playgrounds in Europe. In this study, 355 samples were collected from playgrounds from 16 cities in Hungary. Antibiotic susceptibility of the isolates was tested for nine antibiotics. Presence of virulence factors was detected by PCR. Clonal diversity of the isolates was tested by PFGE and MLST. The overall prevalence of S. aureus was 2.81% (10/355) and no MRSA isolates were found. Presence of spa (10), fnbA (10), fnbB (5), icaA (8), cna (7), sea (2), hla (10), hlb (2) and hlg (6) virulence genes were detected. The isolates had diverse PFGE pulsotypes. With MLST, we have detected isolates belonging to ST8 (CC8), ST22 (CC22), ST944 and ST182 (CC182), ST398 (CC398), ST6609 (CC45), ST3029 and ST2816. We have identified a new sequence type, ST6609 of CC45. S. aureus isolates are present on Hungarian playgrounds, especially on plastic surfaces. The isolates were clonally diverse and showed resistance to commonly used antibiotics. These data reinforce the importance of the outdoor environment in the spread for S. aureus in the community.

Approximate Bayesian Computation applied to time series of population genetic data disentangles rapid genetic changes and demographic variations in a pathogen population

Adaptation can occur at remarkably short timescales in natural populations, leading to drastic changes in phenotypes and genotype frequencies over a few generations only. The inference of demographic parameters can allow understanding how evolutionary forces interact and shape the genetic trajectories of populations during rapid adaptation. Here we propose a new Approximate Bayesian Computation (ABC) framework that couples a forward and individual-based model with temporal genetic data to disentangle genetic changes and demographic variations in a case of rapid adaptation. We test the accuracy of our inferential framework and evaluate the benefit of considering a dense versus sparse sampling. Theoretical investigations demonstrate high accuracy in both model and parameter estimations, even if a strong thinning is applied to time series data. Then, we apply our ABC inferential framework to empirical data describing the population genetic changes of the poplar rust pathogen following a major event of resistance overcoming. We successfully estimate key demographic and genetic parameters, including the proportion of resistant hosts deployed in the landscape and the level of standing genetic variation from which selection occurred. Inferred values are in accordance with our empirical knowledge of this biological system. This new inferential framework, which contrasts with coalescent-based ABC analyses, is promising for a better understanding of evolutionary trajectories of populations subjected to rapid adaptation.

Vertically inherited microbiota and environment modifying behaviours conceal genetic variation in dung beetle life history

Diverse organisms actively manipulate their (sym)biotic and physical environment in ways that feed back on their own development. However, the degree to which these processes affect microevolution remains poorly understood. The gazelle dung beetle both physically modifies its ontogenetic environment and structures its biotic interactions through vertical symbiont transmission. By experimentally eliminating (i) physical environmental modifications and (ii) the vertical inheritance of microbes, we assess how environment modifying behaviour and microbiome transmission shape heritable variation and evolutionary potential. We found that depriving larvae of symbionts and environment modifying behaviours increased additive genetic variance and heritability for development time but not body size. This suggests that larvae's ability to manipulate their environment has the potential to modify heritable variation and to facilitate the accumulation of cryptic genetic variation. This cryptic variation may become released and selectable when organisms encounter environments that are less amenable to organismal manipulation or restructuring. Our findings also suggest that intact microbiomes, which are commonly thought to increase genetic variation of their hosts, may instead reduce and conceal heritable variation. More broadly, our findings highlight that the ability of organisms to actively manipulate their environment may affect the potential of populations to evolve when encountering novel, stressful conditions.

Boom-bust population dynamics drive rapid genetic change

Increasing environmental threats and more extreme environmental perturbations place species at risk of population declines, with associated loss of genetic diversity and evolutionary potential. While theory shows that rapid population declines can cause loss of genetic diversity, populations in some environments, like Australia's arid zone, are repeatedly subject to major population fluctuations yet persist and appear able to maintain genetic diversity. Here, we use repeated population sampling over 13 y and genotype-by-sequencing of 1903 individuals to investigate the genetic consequences of repeated population fluctuations in two small mammals in the Australian arid zone. The sandy inland mouse (Pseudomys hermannsburgensis) experiences marked boom-bust population dynamics in response to the highly variable desert environment. We show that heterozygosity levels declined, and population differentiation (FST) increased, during bust periods when populations became small and isolated, but that heterozygosity was rapidly restored during episodic population booms. In contrast, the lesser hairy-footed dunnart (Sminthopsis youngsoni), a desert marsupial that maintains relatively stable population sizes, showed no linear declines in heterozygosity. These results reveal two contrasting ways in which genetic diversity is maintained in highly variable environments. In one species, diversity is conserved through the maintenance of stable population sizes across time. In the other species, diversity is conserved through rapid genetic mixing during population booms that restores heterozygosity lost during population busts.

Multiple genetic variants at the SLC30A8 locus affect local super-enhancer activity and influence pancreatic β-cell survival and function

Variants at the SLC30A8 locus are associated with type 2 diabetes (T2D) risk. The lead variant, rs13266634, encodes an amino acid change, Arg325Trp (R325W), at the C-terminus of the secretory granule-enriched zinc transporter, ZnT8. Although this protein-coding variant was previously thought to be the sole driver of T2D risk at this locus, recent studies have provided evidence for lowered expression of SLC30A8 mRNA in protective allele carriers. In the present study, we examined multiple variants that influence SLC30A8 allele-specific expression. Epigenomic mapping has previously identified an islet-selective enhancer cluster at the SLC30A8 locus, hosting multiple T2D risk and cASE associations, which is spatially associated with the SLC30A8 promoter and additional neighboring genes. Here, we show that deletion of variant-bearing enhancer regions using CRISPR-Cas9 in human-derived EndoC-βH3 cells lowers the expression of SLC30A8 and several neighboring genes and improves glucose-stimulated insulin secretion. While downregulation of SLC30A8 had no effect on beta cell survival, loss of UTP23, RAD21, or MED30 markedly reduced cell viability. Although eQTL or cASE analyses in human islets did not support the association between these additional genes and diabetes risk, the transcriptional regulator JQ1 lowered the expression of multiple genes at the SLC30A8 locus and enhanced stimulated insulin secretion.

Unraveling the phylogenetics of genetically closely related species, Haemaphysalis japonica and Haemaphysalis megaspinosa, using entire tick mitogenomes and microbiomes

Ticks have a profound impact on public health. Haemaphysalis is one of the most widespread genera in Asia, including Japan. The taxonomy and genetic differentiation of Haemaphysalis spp. is challenging. For instance, previous studies struggled to distinguish Haemaphysalis japonica and Haemaphysalis megaspinosa due to the dearth of nucleotide sequence polymorphisms in widely used barcoding genes. The classification of H. japonica japonica and its related sub-species Haemaphysalis japonica douglasi or Haemaphysalis jezoensis is also confused due to their high morphological similarity and a lack of molecular data that support the current classification. We used mitogenomes and microbiomes of H. japonica and H. megaspinosa to gain deeper insights into the phylogenetic relationships and genetic divergence between two species. Phylogenetic analyses of concatenated nucleotide sequences of protein-coding genes and ribosomal DNA genes distinguished H. japonica and H. megaspinosa as monophyletic clades, with further subdivision within the H. japonica clade. The 16S rRNA and NAD5 genes were valuable markers for distinguishing H. japonica and H. megaspinosa. Population genetic structure analyses indicated that genetic variation within populations accounted for a large proportion of the total variation compared to variation between populations. Microbiome analyses revealed differences in alpha and beta diversity between H. japonica and H. megaspinosa: H. japonica had the higher diversity. Coxiella sp., a likely endosymbiont, was found in both Haemaphysalis species. The abundance profiles of likely endosymbionts, pathogens, and commensals differed between H. japonica and H. megaspinosa: H. megaspinosa was more diverse.

Copy number and selection of MHC genes in ruminants are related to habitat, average life span and diet

The ruminants, as the main group of livestock, have been extensively studied in terms of their physiology, endocrinology, biochemistry, genetics, and nutrition. Despite the wide geographic distribution and habitat diversity of animals in this group, their ecology and evolution remain poorly understood. In this study, we analyzed the gene copy number, selection, and ecological and evolutionary processes that have affected the evolution of major histocompatibility complex (MHC) genes across ruminant lineages based on available genomic data. The 51 species analyzed represented all six families of ruminants. Our finding indicated that the architecture of the MHC region is conserved in ruminants, but with variable copy numbers of MHC-I, MHC-IIA, and MHC-IIB genes. No lineage-specific gene duplication was observed in the MHC genes. The phylogenetic generalized least squares regression (PGLS) model revealed association between ecological and biological factors (habitat and lifespan) and gene duplication in DQA and DQB, but not in DRB. The selection pressure of DQA and DQB were related with lifespan, diet, and the ratio of genetic repeat elements. These results suggest that the MHC evolution in ruminants, including copy number and selection, has been influenced by genetic repeat elements, pathogen exposure risk, and intrinsic cost of possessing multiple MHC genes.

Harnessing the intragenomic variability of rRNA operons to improve differentiation of Vibrio species

Although the 16S rRNA gene is frequently used as a phylogenetic marker in analysis of environmental DNA, this marker often fails to distinguish closely related species, including those in the genus Vibrio. Here, we investigate whether inclusion and analysis of 23S rRNA sequence can help overcome the intrinsic weaknesses of 16S rRNA analyses for the differentiation of Vibrio species. We construct a maximum likelihood 16S rRNA gene tree to assess the use of this gene to identify clades of Vibrio species. Within the 16S rRNA tree, we identify the putative informative bases responsible for polyphyly, and demonstrate the association of these positions with tree topology. We demonstrate that concatenation of 16S and 23S rRNA genes increases the number of informative nucleotide positions, thereby overcoming ambiguities in 16S rRNA-based phylogenetic reconstructions. Finally, we experimentally demonstrate that this approach considerably improves the differentiation and identification of Vibrio species in environmental samples.

Finger millet (Eleusine coracana) enhancement through genomic resources and breeding methods: current implications and potential future interventions

Finger millet (Eleusine coracana) is an essential staple crop in many regions of Africa and Asia, valued for its nutritional content and resilience in challenging agro-ecological conditions. The enhancement of finger millet through genomic resources and breeding methods represents a promising avenue for addressing food and nutritional security. Current efforts in this field have harnessed genomic technologies to decipher the crop's genetic diversity and identify key traits related to yield, disease resistance, and nutritional content. These insights have facilitated the development of improved varieties through selective breeding, accelerating the crop's adaptation to changing environmental conditions. In the future, continued advancements in genomics and breeding methodologies hold the potential to further enhance finger millet's resilience, nutritional value, and productivity, ultimately benefiting both farmers and consumers. This review article synthesizes the current state of research and development in finger millet enhancement through the integration of genomic resources and innovative breeding methods. The utilization of these insights in selective breeding has already yielded promising results in developing improved finger millet varieties that meet the evolving needs of farmers and consumers. Moreover, this article discusses potential future interventions, including the continued advancement of genomics, precision breeding, and sustainable agricultural practices. These interventions hold the promise of further enhancing finger millet's adaptability to changing climates, its nutritional quality, and its overall productivity, thereby contributing to food security and improved livelihoods in finger millet-dependent regions.

Serial founder effects slow range expansion in an invasive social insect

Invasive populations often experience founder effects: a loss of genetic diversity relative to the source population, due to a small number of founders. Even where these founder effects do not impact colonization success, theory predicts they might affect the rate at which invasive populations expand. This is because secondary founder effects are generated at advancing population edges, further reducing local genetic diversity and elevating genetic load. We show that in an expanding invasive population of the Asian honey bee (Apis cerana), genetic diversity is indeed lowest at range edges, including at the complementary sex determiner, csd, a locus that is homozygous-lethal. Consistent with lower local csd diversity, range edge colonies had lower brood viability than colonies in the range centre. Further, simulations of a newly-founded and expanding honey bee population corroborate the spatial patterns in mean colony fitness observed in our empirical data and show that such genetic load at range edges will slow the rate of population expansion.

Stalk rot species diversity and molecular phylogeny associated with diseased maize in India

Stalk rot disease is a major constraint in maize production and till date reported to be caused by two to three species of phytopathogenic fungi but, in our present study, we disclose the first report of stalk rot is caused by complex species of phytopathogens, which belongs to five different genera. Therefore, to substantiate these findings, a total of 105 diseased samples of maize were collected from 21 different locations in six different geographical locations of India from which 48 isolates were used for the research study. Morphological features such as pigmentation, colony color, type of mycelium and pattern of mycelium was examined using macro and microscopic methods. A total of 11 different spp. of pathogens belonging to the five different genera: Fusarium verticillioides (56.25%), F. equiseti (14.5%), F. andiyazi (6.25%), F. solani (2.08%), F. proliferatum (2.08%), F. incarnatum (2.08%), Lasidioplodia theobrame (6.25%), Exserohilum rostrtum (4.16%), Nigrospora spp. (4.16%). and Schizophyllum commune (2.08%) were identified by different housekeeping genes (ITS, TEF-1α, RPB2 and Actin). Fusarium verticillioides, F. equiseti and F. andiyazi were major pathogens involved in stalk rot. This is the first report on F. proliferatum, F. solani, F. incarnatum, Lasidioplodia theobrame, Exserohilum rostrtum, Nigrospora spp. and Schizophyllum commune causing stalk rot of maize and their distribution in the different states of India. Studies on population dynamics of PFSR will enhance the understanding of pathogen behavior, virulence, or its association with different pathogens across India, which will facilitate the development of resistant maize genotypes against the PFSR.

Unprecedented variation pattern of plastid genomes and the potential role in adaptive evolution in Poales

BACKGROUND

The plastid is the photosynthetic organelle in plant cell, and the plastid genomes (plastomes) are generally conserved in evolution. As one of the most economically and ecologically important order of angiosperms, Poales was previously documented to exhibit great plastomic variation as an order of photoautotrophic plants.

RESULTS

We acquired 93 plastomes, representing all the 16 families and 5 major clades of Poales to reveal the extent of their variation and evolutionary pattern. Extensive variation including the largest one in monocots with 225,293 bp in size, heterogeneous GC content, and a wide variety of gene duplication and loss were revealed. Moreover, rare occurrences of three inverted repeat (IR) copies in angiosperms and one IR loss were observed, accompanied by short IR (sIR) and small direct repeat (DR). Widespread structural heteroplasmy, diversified inversions, and unusual genomic rearrangements all appeared in Poales, occasionally within a single species. Extensive repeats in the plastomes were found to be positively correlated with the observed inversions and rearrangements. The variation all showed a "small-large-moderate" trend along the evolution of Poales, as well as for the sequence substitution rate. Finally, we found some positively selected genes, mainly in C4 lineages, while the closely related lineages of those experiencing gene loss tended to have undergone more relaxed purifying selection.

CONCLUSIONS

The variation of plastomes in Poales may be related to its successful diversification into diverse habitats and multiple photosynthetic pathway transitions. Our order-scale analyses revealed unusual evolutionary scenarios for plastomes in the photoautotrophic order of Poales and provided new insights into the plastome evolution in angiosperms as a whole.

Using Portuguese BRCA pathogenic variation as a model to study the impact of human admixture on human health

BACKGROUND

Admixture occurs between different ethnic human populations. The global colonization in recent centuries by Europeans led to the most significant admixture in human history. While admixture may enhance genetic diversity for better fitness, it may also impact on human health by transmitting genetic variants for disease susceptibility in the admixture population. The admixture by Portuguese global exploration initiated in the 15th century has reached over 20 million of Portuguese-heritage population worldwide. It provides a valuable model to study the impact of admixture on human health. BRCA1 and BRCA2 (BRCA) are two of the important tumor suppressor genes. The pathogenic variation (PV) in BRCA is well determined to cause high risk of hereditary breast and ovarian cancer. Tracing the distribution of Portuguese BRCA PV in Portuguese-heritage population will help to understand the impact of admixture on cancer susceptibility in modern humans. In this study, we analyzed the distribution of the Portuguese-originated BRCA variation in Brazilian population, which has high degree Portuguese-heritage.

METHODS

By comprehensive data mining, standardization and annotation, we generated a Portuguese-derived BRCA variation dataset and a Brazilian-derived BRCA variation dataset. We compared the two BRCA variation datasets to identify the BRCA variants shared between the two populations.

RESULTS

The Portuguese-derived BRCA variation dataset consists of 220 BRCA variants including 78 PVs from 11,482 Portuguese cancer patients, 93 (42.2%) in BRCA1 and 127 (57.7%) in BRCA2. Of the 556 Portuguese BRCA PV carriers carrying the 78 PVs, 331 (59.5%) carried the three Portuguese-BRCA founder PVs of BRCA1 c.2037delinsCC, BRCA1 c.3331_3334del and BRCA2 c.156_157insAlu. The Brazilian-derived BRCA variation dataset consists of 255 BRCA PVs from 7,711 cancer patients, 136 (53.3%) in BRCA1 and 119 (46.6%) in BRCA2. We developed an open database named dbBRCA-Portuguese ( https://genemutation.fhs.um.edu.mo/dbbrca-portuguese/ ) and an open database named dbBRCA-Brazilian ( https://genemutation.fhs.um.edu.mo/dbbrca-brazilian ) to host the BRCA variation data from Portuguese and Brazilian populations. We compared the BRCA PV datasets between Portuguese and Brazilian populations, and identified 29 Portuguese-specific BRCA PVs shared between Portuguese and Brazilian populations, 14 in BRCA1 including the Portuguese founder BRCA1 c.3331_3334del and BRCA1 c.2037delinsCC, and 15 in BRCA2 including the Portuguese founder BRCA2 c.156_157insAlu. Searching the 78 Portuguese BRCA PVs in over 5,000 ancient human genomes identified evolution origin for only 8 PVs in Europeans dated between 37,470 and 3,818 years before present, confirming the Portuguese-specificity of Portuguese BRCA PVs; comparing the 78 Portuguese BRCA PVs Portuguese, 255 Brazilian BRCA PVs, and 134 African BRCA PVs showed little overlapping, ruling out the possibility that the BRCA PVs shared between Portuguese and Brazilian may also be contributed by African.

CONCLUSION

Our study provides evidence that the admixture in recent human history contributed to cancer susceptibility in modern humans.

Imputation of 3D genome structure by genetic-epigenetic interaction modeling in mice

Gene expression is known to be affected by interactions between local genetic variation and DNA accessibility, with the latter organized into three-dimensional chromatin structures. Analyses of these interactions have previously been limited, obscuring their regulatory context, and the extent to which they occur throughout the genome. Here, we undertake a genome-scale analysis of these interactions in a genetically diverse population to systematically identify global genetic-epigenetic interaction, and reveal constraints imposed by chromatin structure. We establish the extent and structure of genotype-by-epigenotype interaction using embryonic stem cells derived from Diversity Outbred mice. This mouse population segregates millions of variants from eight inbred founders, enabling precision genetic mapping with extensive genotypic and phenotypic diversity. With 176 samples profiled for genotype, gene expression, and open chromatin, we used regression modeling to infer genetic-epigenetic interactions on a genome-wide scale. Our results demonstrate that statistical interactions between genetic variants and chromatin accessibility are common throughout the genome. We found that these interactions occur within the local area of the affected gene, and that this locality corresponds to topologically associated domains (TADs). The likelihood of interaction was most strongly defined by the three-dimensional (3D) domain structure rather than linear DNA sequence. We show that stable 3D genome structure is an effective tool to guide searches for regulatory elements and, conversely, that regulatory elements in genetically diverse populations provide a means to infer 3D genome structure. We confirmed this finding with CTCF ChIP-seq that revealed strain-specific binding in the inbred founder mice. In stem cells, open chromatin participating in the most significant regression models demonstrated an enrichment for developmental genes and the TAD-forming CTCF-binding complex, providing an opportunity for statistical inference of shifting TAD boundaries operating during early development. These findings provide evidence that genetic and epigenetic factors operate within the context of 3D chromatin structure.

A hepatitis B virus (HBV) sequence variation graph improves alignment and sample-specific consensus sequence construction

Nearly 300 million individuals live with chronic hepatitis B virus (HBV) infection (CHB), for which no curative therapy is available. As viral diversity is associated with pathogenesis and immunological control of infection, improved methods to characterize this diversity could aid drug development efforts. Conventionally, viral sequencing data are mapped/aligned to a reference genome, and only the aligned sequences are retained for analysis. Thus, reference selection is critical, yet selecting the most representative reference a priori remains difficult. We investigate an alternative pangenome approach which can combine multiple reference sequences into a graph which can be used during alignment. Using simulated short-read sequencing data generated from publicly available HBV genomes and real sequencing data from an individual living with CHB, we demonstrate alignment to a phylogenetically representative 'genome graph' can improve alignment, avoid issues of reference ambiguity, and facilitate the construction of sample-specific consensus sequences more genetically similar to the individual's infection. Graph-based methods can, therefore, improve efforts to characterize the genetics of viral pathogens, including HBV, and have broader implications in host-pathogen research.

A systematic review on the biochemical threshold of mitochondrial genetic variants

Mitochondrial DNA (mtDNA) variants cause a range of diseases from severe pediatric syndromes to aging-related conditions. The percentage of mtDNA copies carrying a pathogenic variant, variant allele frequency (VAF), must reach a threshold before a biochemical defect occurs, termed the biochemical threshold. Whether the often-cited biochemical threshold of >60% VAF is similar across mtDNA variants and cell types is unclear. In our systematic review, we sought to identify the biochemical threshold of mtDNA variants in relation to VAF by human tissue/cell type. We used controlled vocabulary terms to identify articles measuring oxidative phosphorylation (OXPHOS) complex activities in relation to VAF. We identified 76 eligible publications, describing 69, 12, 16, and 49 cases for complexes I, III, IV, and V, respectively. Few studies evaluated OXPHOS activities in diverse tissue types, likely reflective of clinical access. A number of cases with similar VAFs for the same pathogenic variant had varying degrees of residual activity of the affected complex, alluding to the presence of modifying variants. Tissues and cells with VAFs <60% associated with low complex activities were described, suggesting the possibility of a biochemical threshold of <60%. Using Kendall rank correlation tests, the VAF of the m.8993T > G variant correlated with complex V activity in skeletal muscle (τ = -0.58, P = 0.01, n = 13); however, no correlation was observed in fibroblasts (P = 0.7, n = 9). Our systematic review highlights the need to investigate the biochemical threshold over a wider range of VAFs in disease-relevant cell types to better define the biochemical threshold for specific mtDNA variants.

Association of uncertain significance genetic variants with myocardial mechanics and morphometrics in patients with nonischemic dilated cardiomyopathy

BACKGROUND

Careful interpretation of the relation between phenotype changes of the heart and gene variants detected in dilated cardiomyopathy (DCM) is important for patient care and monitoring.

OBJECTIVE

We sought to assess the association between cardiac-related genes and whole-heart myocardial mechanics or morphometrics in nonischemic dilated cardiomyopathy (NIDCM).

METHODS

It was a prospective study consisting of patients with NIDCM. All patients were referred for genetic testing and a genetic analysis was performed using Illumina NextSeq 550 and a commercial gene capture panel of 233 genes (Systems Genomics, Cardiac-GeneSGKit®). It was analyzed whether there are significant differences in clinical, two-dimensional (2D) echocardiographic, and magnetic resonance imaging (MRI) parameters between patients with the genes variants and those without. 2D echocardiography and MRI were used to analyze myocardial mechanics and morphometrics.

RESULTS

The study group consisted of 95 patients with NIDCM and the average age was 49.7 ± 10.5. All echocardiographic and MRI parameters of myocardial mechanics (left ventricular ejection fraction 28.4 ± 8.7 and 30.7 ± 11.2, respectively) were reduced and all values of cardiac chambers were increased (left ventricular end-diastolic diameter 64.5 ± 5.9 mm and 69.5 ± 10.7 mm, respectively) in this group. It was noticed that most cases of whole-heart myocardial mechanics and morphometrics differences between patients with and without gene variants were in the genes GATAD1, LOX, RASA1, KRAS, and KRIT1. These genes have not been previously linked to DCM. It has emerged that KRAS and KRIT1 genes were associated with worse whole-heart mechanics and enlargement of all heart chambers. GATAD1, LOX, and RASA1 genes variants showed an association with better cardiac function and morphometrics parameters. It might be that these variants alone do not influence disease development enough to be selective in human evolution.

CONCLUSIONS

Combined variants in previously unreported genes related to DCM might play a significant role in affecting clinical, morphometrics, or myocardial mechanics parameters.

Genetic variant rs1205 is associated with COVID-19 outcomes: The Strong Heart Study and Strong Heart Family Study

BACKGROUND

Although COVID-19 infection has been associated with a number of clinical and environmental risk factors, host genetic variation has also been associated with the incidence and morbidity of infection. The CRP gene codes for a critical component of the innate immune system and CRP variants have been reported associated with infectious disease and vaccination outcomes. We investigated possible associations between COVID-19 outcome and a limited number of candidate gene variants including rs1205.

METHODOLOGY/PRINCIPAL FINDINGS

The Strong Heart and Strong Heart Family studies have accumulated detailed genetic, cardiovascular risk and event data in geographically dispersed American Indian communities since 1988. Genotypic data and 91 COVID-19 adjudicated deaths or hospitalizations from 2/1/20 through 3/1/23 were identified among 3,780 participants in two subsets. Among 21 candidate variants including genes in the interferon response pathway, APOE, TMPRSS2, TLR3, the HLA complex and the ABO blood group, only rs1205, a 3' untranslated region variant in the CRP gene, showed nominally significant association in T-dominant model analyses (odds ratio 1.859, 95%CI 1.001-3.453, p = 0.049) after adjustment for age, sex, center, body mass index, and a history of cardiovascular disease. Within the younger subset, association with the rs1205 T-Dom genotype was stronger, both in the same adjusted logistic model and in the SOLAR analysis also adjusting for other genetic relatedness.

CONCLUSION

A T-dominant genotype of rs1205 in the CRP gene is associated with COVID-19 death or hospitalization, even after adjustment for relevant clinical factors and potential participant relatedness. Additional study of other populations and genetic variants of this gene are warranted.

[Characterization of sequences, expression profiling, and natural allelic variation analysis of the MYC gene family in sorghum (Sorghum bicolor)]

Sorghum aphid (Melanaphis sacchari) and head smut fungi (Sporisorium reilianum) infesting sorghum cause delayed growth and development, and reduce yield and quality. This study use bioinformatics and molecular biological approaches to profile the gene expression pattern during sorghum development and under pest infestation, and analyzed the natural allelic DNA variation of sorghum MYC gene family. The findings provide insights for potential application in breeding the stress resistant and high productivity sorghum varieties. The results indicated that there are 28 MYC genes identified in sorghum genome, distributed on 10 chromosomes. The bHLH_MYC_N and HLH domains are the conserved domains of the MYC gene in sorghum. Gene expression analysis showed that SbbHLH35.7g exhibited high expression levels in leaves, SbAbaIn showed strong expression in early grains, and SbMYC2.1g showed high expression levels in mature pollen. In anti-aphid strains at the 5-leaf stage, SbAbaIn, SbLHW.4g and SbLHW.2g were significantly induced in leaves, while SbbHLH35.7g displayed the highest expression level in panicle tissue, which was significantly induced by the infection of head smut. Promoter cis-element analysis identified methyl jasmonate (MJ), abscisic acid (ABA), salicylic acid (SA) and MYB-binding sites related to drought-stress inducibility. Furthermore, genomic resequencing data analysis revealed natural allelic DNA variations such as single nucleotide polymorphism (SNP) and insertion-deletion (INDEL) for the key SbMYCs. Protein interaction network analysis using STRING indicated that SbAbaIn interacts with TIFYdomain protein, and SbbHLH35.7g interacts with MDR and imporin. SbMYCs exhibited temporal and spatial expression patterns and played vital roles during the sorghum development. Infestation by sugarcane aphids and head smut fungi induced the expression of SbAbaIn and SbbHLH35.7g, respectively. SbAbaIn modulated the jasmonic acid (JA) pathway to regulate the expression of defensive genes, conferring resistance to insects. On the other hand, SbbHLH35.7g participated in detoxification reactions to defend against pathogens.

Deep learning the cis-regulatory code for gene expression in selected model plants

Elucidating the relationship between non-coding regulatory element sequences and gene expression is crucial for understanding gene regulation and genetic variation. We explored this link with the training of interpretable deep learning models predicting gene expression profiles from gene flanking regions of the plant species Arabidopsis thaliana, Solanum lycopersicum, Sorghum bicolor, and Zea mays. With over 80% accuracy, our models enabled predictive feature selection, highlighting e.g. the significant role of UTR regions in determining gene expression levels. The models demonstrated remarkable cross-species performance, effectively identifying both conserved and species-specific regulatory sequence features and their predictive power for gene expression. We illustrated the application of our approach by revealing causal links between genetic variation and gene expression changes across fourteen tomato genomes. Lastly, our models efficiently predicted genotype-specific expression of key functional gene groups, exemplified by underscoring known phenotypic and metabolic differences between Solanum lycopersicum and its wild, drought-resistant relative, Solanum pennellii.