BCFtools/csq: haplotype-aware variant consequences

AI-guided identification of risk variants for adrenocortical tumours in TP53 p.R337H carrier children: a genetic association study

Background

Adrenocortical tumours (ACT) in children are part of the Li-Fraumeni cancer spectrum and are frequently associated with a germline TP53 pathogenic variant. TP53 p.R337H is highly prevalent in the south and southeast of Brazil and predisposes to ACT with low penetrance. Thus, we aimed to investigate whether genetic variants exist which are associated with an increased risk of developing ACT in TP53 p.R337H carrier children.

Methods

A genetic association study was conducted in trios of children (14 girls, 7 boys) from southern Brazil carriers of TP53 p.R337H with (n = 18) or without (n = 3) ACT and their parents, one of whom also carries this pathogenic variant (discovery cohort). Results were confirmed in a validation cohort of TP53 p.R337H carriers with (n = 90; 68 girls, 22 boys) or without ACT (n = 302; 165 women, 137 men).

Findings

We analysed genomic data from whole exome sequencing of blood DNA from the trios. Using deep learning algorithms, according to a model where the affected child inherits from the non-carrier parent variant(s) increasing the risk of developing ACT, we found a significantly enriched representation of non-coding variants in genes involved in the cyclic AMP (cAMP) pathway known to be involved in adrenocortical tumorigenesis. One among those variants (rs2278986 in the SCARB1 gene) was confirmed to be significantly enriched in the validation cohort of TP53 p.R337H carriers with ACT compared to carriers without ACT (OR 1.858; 95% CI 1.146, 3.042, p = 0.01).

Interpretation

Profiling of the variant rs2278986 is a candidate for future confirmation and possible use as a tool for ACT risk stratification in TP53 p.R337H carriers.

Funding

Centre National de la Recherche Scientifique (CNRS), Behring Foundation, Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq).

A reference genome for the Harpy Eagle reveals steady demographic decline and chromosomal rearrangements in the origin of Accipitriformes

The Harpy Eagle (Harpia harpyja) is an iconic species that inhabits forested landscapes in Neotropical regions, with decreasing population trends mainly due to habitat loss, and currently classified as vulnerable. Here, we report on a chromosome-scale genome assembly for a female individual combining long reads, optical mapping, and chromatin conformation capture reads. The final assembly spans 1.35 Gb, with N50scaffold equal to 58.1 Mb and BUSCO completeness of 99.7%. We built the first extensive transposable element (TE) library for the Accipitridae to date and identified 7,228 intact TEs. We found a burst of an unknown TE ~ 13-22 million years ago (MYA), coincident with the split of the Harpy Eagle from other Harpiinae eagles. We also report a burst of solo-LTRs and CR1 retrotransposons ~ 31-33 MYA, overlapping with the split of the ancestor to all Harpiinae from other Accipitridae subfamilies. Comparative genomics with other Accipitridae, the closely related Cathartidae and Galloanserae revealed major chromosome-level rearrangements at the basal Accipitriformes genome, in contrast to a conserved ancient genome architecture for the latter two groups. A historical demography reconstruction showed a rapid decline in effective population size over the last 20,000 years. This reference genome serves as a crucial resource for future conservation efforts towards the Harpy Eagle.

Distinct Genes with Similar Functions Underlie Convergent Evolution in Myotis Bat Ecomorphs

Convergence offers an opportunity to explore to what extent evolution can be predictable when genomic composition and environmental triggers are similar. Here, we present an emergent model system to study convergent evolution in nature in a mammalian group, the bat genus Myotis. Three foraging strategies-gleaning, trawling, and aerial hawking, each characterized by different sets of phenotypic features-have evolved independently multiple times in different biogeographic regions in isolation for millions of years. To investigate the genomic basis of convergence and explore the functional genomic changes linked to ecomorphological convergence, we sequenced and annotated 17 new genomes and screened 16,426 genes for positive selection and associations between relative evolutionary rates and foraging strategies across 30 bat species representing all Myotis ecomorphs across geographic regions as well as among sister groups. We identify genomic changes that describe both phylogenetic and ecomorphological trends. We infer that colonization of new environments may have first required changes in genes linked to hearing sensory perception, followed by changes linked to fecundity and development, metabolism of carbohydrates, and heme degradation. These changes may be linked to prey acquisition and digestion and match phylogenetic trends. Our findings also suggest that the repeated evolution of ecomorphs does not always involve changes in the same genes but rather in genes with the same molecular functions such as developmental and cellular processes.

Resistance of zebu cattle (Bos indicus) to colonization by major ruminant hoof pathogens

Zebu cattle (Bos indicus) is reported to be more resistant towards harmful environmental factors than taurine cattle (Bos taurus). A few hundred zebu cattle are kept in Switzerland and in contrast to the Swiss indigenous breeds, infectious hoof disease in zebu is not observed. Therefore, we compared the prevalence of three ruminant hoof pathogens in zebu and taurine cattle. These included Treponema spp., Fusobacterium necrophorum and Dichelobacter nodosus which are associated with bovine digital dermatitis (BDD), different bovine hoof diseases and ovine footrot, respectively. Interdigital swabs and punch biopsies from hind feet of slaughter animals were tested for the three pathogens by PCR. Sixty zebu from eight farms were compared to a convenience sample of 20 taurine cattle from 17 farms. Treponema spp. associated with BDD were not detected in zebu while 23 % of animals and 50 % of farms were positive for benign D. nodosus, with results indicating environmental contamination rather than colonization. Taurine cattle showed 35 % of animals and 41 % of farms positive for T. phagedenis while 90 % of animals and 94 % of farms were colonized by D. nodosus as indicated by a 500-fold higher bacterial load than in zebu. The difference in prevalence of the two pathogens between zebu and taurine cattle was highly significant. F. necrophorum was as well only detected in taurine cattle with values of 15 % of animals and 17.7 % of farms, being significantly different at the animal level. Furthermore, genetic analysis of Swiss zebu indicates high genomic diversity and clear separation from taurine cattle. This is the first evidence that zebu show resistance towards colonization by bacterial hoof pathogens in contrast to taurine cattle.

HASCH - A high-throughput amplicon-based SNP-platform for medicinal cannabis and industrial hemp genotyping applications

BACKGROUND

Cannabis sativa is seeing a global resurgence as a food, fiber and medicinal crop for industrial hemp and medicinal Cannabis industries respectively. However, a widespread moratorium on the use and research of C. sativa throughout most of the 20th century has seen the development of improved cultivars for specific end uses lag behind that of conventional crops. While C. sativa research and development has seen significant investments in the recent past, resulting in a suite of publicly available genomic resources and tools, a versatile and cost-effective mid-density genotyping platform for applied purposes in breeding and pre-breeding is lacking. Here we report on a first mid-density fixed-target SNP platform for C. sativa.

RESULTS

The High-throughput Amplicon-based SNP-platform for medicinal Cannabis and industrial Hemp (HASCH) was designed using a combination of filtering and Integer Linear Programming on publicly available whole-genome sequencing and RNA sequencing data, supplemented with in-house generated genotyping-by-sequencing (GBS) data. HASCH contains 1,504 genome-wide targets of high call rate (97% mean) and even distribution across the genome, designed to be highly informative (> 0.3 minor allele frequency) across both medicinal cannabis and industrial hemp gene pools. Average numbers of mismatch SNP between any two accessions were 251 for medicinal cannabis (N = 116) and 272 for industrial hemp (N = 87). Comparing HASCH data with corresponding GBS data on a collection of diverse C. sativa accessions demonstrated high concordance and resulted in comparable phylogenies and genetic distance matrices. Using HASCH on a segregating F2 population derived from a cross between a tetrahydrocannabinol (THC)-dominant and a cannabidiol (CBD)-dominant accession resulted in a genetic map consisting of 310 markers, comprising 10 linkage groups and a total size of 582.7 cM. Quantitative Trait Locus (QTL) mapping identified a major QTL for CBD content on chromosome 7, consistent with previous findings.

CONCLUSION

HASCH constitutes a versatile, easy to use and cost-effective genotyping solution for the rapidly growing Cannabis research community. It provides consistent genetic fingerprints of 1504 SNPs with wide applicability genetic resource management, quantitative genetics and breeding.

Detection of distant relatedness in biobanks to identify undiagnosed cases of Mendelian disease as applied to Long QT syndrome

Rare genetic diseases are typically studied in referral populations, resulting in underdiagnosis and biased assessment of penetrance and phenotype. To address this, we develop a generalizable method of genotype inference based on distant relatedness and deploy this to identify undiagnosed Type 5 Long QT Syndrome (LQT5) rare variant carriers in a non-referral population. We identify 9 LQT5 families referred to a single specialty clinic, each carrying p.Asp76Asn, the most common LQT5 variant. We uncover recent common ancestry and a single shared haplotype among probands. Application to a non-referral population of 69,819 BioVU biobank subjects identifies 22 additional subjects sharing this haplotype, which we confirm to carry p.Asp76Asn. Referral and non-referral carriers have prolonged QT interval corrected for heart rate (QTc) compared to controls, and, among carriers, the QTc polygenic score is independently associated with QTc prolongation. Thus, our innovative analysis of shared chromosomal segments identifies undiagnosed cases of genetic disease and refines the understanding of LQT5 penetrance and phenotype.

Complete mitochondrial genome assembly of Zizania latifolia and comparative genome analysis

Zizania latifolia (Griseb.) Turcz. ex Stapf has been cultivated as a popular aquatic vegetable in China due to its important nutritional, medicinal, ecological, and economic values. The complete mitochondrial genome (mitogenome) of Z. latifolia has not been previously studied and reported, which has hindered its molecular systematics and understanding of evolutionary processes. Here, we assembled the complete mitogenome of Z. latifolia and performed a comprehensive analysis including genome organization, repetitive sequences, RNA editing event, intercellular gene transfer, phylogenetic analysis, and comparative mitogenome analysis. The mitogenome of Z. latifolia was estimated to have a circular molecule of 392,219 bp and 58 genes consisting of three rRNA genes, 20 tRNA genes, and 35 protein-coding genes (PCGs). There were 46 and 20 simple sequence repeats (SSRs) with different motifs identified from the mitogenome and chloroplast genome of Z. latifolia, respectively. Furthermore, 49 homologous fragments were observed to transfer from the chloroplast genome to the mitogenome of Z. latifolia, accounting for 47,500 bp, presenting 12.1% of the whole mitogenome. In addition, there were 11 gene-containing homologous regions between the mitogenome and chloroplast genome of Z. latifolia. Also, approximately 85% of fragments from the mitogenome were duplicated in the Z. latifolia nuclear genome. Selection pressure analysis revealed that most of the mitochondrial genes were highly conserved except for ccmFc, ccmFn, matR, rps1, and rps3. A total of 93 RNA editing sites were found in the PCGs of the mitogenome. Z. latifolia and Oryza minuta are the most closely related, as shown by collinear analysis and the phylogenetic analysis. We found that repeat sequences and foreign sequences in the mitogenomes of Oryzoideae plants were associated with genome rearrangements. In general, the availability of the Z. latifolia mitogenome will contribute valuable information to our understanding of the molecular and genomic aspects of Zizania.

Molecular Mechanisms of Temperature Tolerance Plasticity in an Arthropod

How species thrive in a wide range of environments is a major focus of evolutionary biology. For many species, limited genetic diversity or gene flow among habitats means that phenotypic plasticity must play an important role in their capacity to tolerate environmental heterogeneity and to colonize new habitats. However, we have a limited understanding of the molecular components that govern plasticity in ecologically relevant phenotypes. We examined this hypothesis in a spider species (Stegodyphus dumicola) with extremely low species-wide genetic diversity that nevertheless occupies a broad range of thermal environments. We determined phenotypic responses to temperature stress in individuals from four climatic zones using common garden acclimation experiments to disentangle phenotypic plasticity from genetic adaptations. Simultaneously, we created data sets on multiple molecular modalities: the genome, the transcriptome, the methylome, the metabolome, and the bacterial microbiome to determine associations with phenotypic responses. Analyses of phenotypic and molecular associations reveal that acclimation responses in the transcriptome and metabolome correlate with patterns of phenotypic plasticity in temperature tolerance. Surprisingly, genes whose expression seemed to be involved in plasticity in temperature tolerance were generally highly methylated contradicting the idea that DNA methylation stabilizes gene expression. This suggests that the function of DNA methylation in invertebrates varies not only among species but also among genes. The bacterial microbiome was stable across the acclimation period; combined with our previous demonstrations that the microbiome is temporally stable in wild populations, this is convincing evidence that the microbiome does not facilitate plasticity in temperature tolerance. Our results suggest that population-specific variation in temperature tolerance among acclimation temperatures appears to result from the evolution of plasticity in mainly gene expression.

Plasmodium falciparum population dynamics in East Africa and genomic surveillance along the Kenya-Uganda border

East African countries accounted for ~ 10% of all malaria prevalence worldwide in 2022, with an estimated 23.8 million cases and > 53,000 deaths. Despite recent increases in malaria incidence, high-resolution genome-wide analyses of Plasmodium parasite populations are sparse in Kenya, Tanzania, and Uganda. The Kenyan-Ugandan border region is a particular concern, with Uganda confirming the emergence and spread of artemisinin resistant P. falciparum parasites. To establish genomic surveillance along the Kenyan-Ugandan border and analyse P. falciparum population dynamics within East Africa, we generated whole-genome sequencing (WGS) data for 38 parasites from Bungoma, Western Kenya. These sequences were integrated into a genomic analysis of available East African isolate data (n = 599) and revealed parasite subpopulations with distinct genetic structure and diverse ancestral origins. Ancestral admixture analysis of these subpopulations alongside isolates from across Africa (n = 365) suggested potential independent ancestral populations from other major African populations. Within isolates from Western Kenya, the prevalence of biomarkers associated with chloroquine resistance (e.g. Pfcrt K76T) were significantly reduced compared to wider East African populations and a single isolate contained the PfK13 V568I variant, potentially linked to reduced susceptibility to artemisinin. Overall, our work provides baseline WGS data and analysis for future malaria genomic surveillance in the region.

Stem lodging Resistance-1 controls stem strength by positively regulating the biosynthesis of cell wall components in Capsicum annuum L

Lodging presents a significant challenge in cultivating high-yield crops with extensive above-ground biomass, yet the molecular mechanisms underlying this phenomenon in the Solanaceae family remain largely unexplored. In this study, we identified a gene, CaSLR1 (Capsicum annuum Stem Lodging Resistance 1), which encodes a MYELOBLASTOSIS (MYB) family transcription factor, from a lodging-affected C. annuum EMS mutant. The suppression of CaSLR1 expression in pepper led to notable stem lodging, reduced thickness of the secondary cell wall, and decreased stem strength. A similar phenotype was observed in tomato with the knockdown of SlMYB61, the orthologous gene to CaSLR1. Further investigations demonstrated that CaNAC6, a gene involved in secondary cell wall (SCW) formation, is co-expressed with CaSLR1 and acts as a positive regulator of its expression, as confirmed through yeast one-hybrid, dual-luciferase reporter assays, and electrophoretic mobility shift assays. These findings elucidate the CaNAC6-CaSLR1 module that contributes to lodging resistance, emphasizing the critical role of CaSLR1 in the lodging resistance regulatory network.

Comparative Analysis Highlights Uniconazole's Efficacy in Enhancing the Cold Stress Tolerance of Mung Beans by Targeting Photosynthetic Pathways

Soybean (Glycine max) and mung bean (Vigna radiata) are key legumes with global importance, but their mechanisms for coping with cold stress-a major challenge in agriculture-have not been thoroughly investigated, especially in a comparative study. This research aimed to fill this gap by examining how these two major legumes respond differently to cold stress and exploring the role of uniconazole, a potential stress mitigator. Our comprehensive approach involved transcriptomic and metabolomic analyses, revealing distinct responses between soybean and mung bean under cold stress conditions. Notably, uniconazole was found to significantly enhance cold tolerance in mung bean by upregulating genes associated with photosynthesis, while its impact on soybean was either negligible or adverse. To further understand the molecular interactions, we utilized advanced machine learning algorithms for protein structure prediction, focusing on photosynthetic pathways. This enabled us to identify LOC106780309 as a direct binding target for uniconazole, confirmed through isothermal titration calorimetry. This research establishes a new comparative approach to explore how soybean and mung bean adapt to cold stress, offers key insights to improve the hardiness of legumes against environmental challenges, and contributes to sustainable agricultural practices and food security.

Diverse Gene Regulatory Mechanisms Alter Rattlesnake Venom Gene Expression at Fine Evolutionary Scales

Understanding and predicting the relationships between genotype and phenotype is often challenging, largely due to the complex nature of eukaryotic gene regulation. A step towards this goal is to map how phenotypic diversity evolves through genomic changes that modify gene regulatory interactions. Using the Prairie Rattlesnake (Crotalus viridis) and related species, we integrate mRNA-seq, proteomic, ATAC-seq and whole-genome resequencing data to understand how specific evolutionary modifications to gene regulatory network components produce differences in venom gene expression. Through comparisons within and between species, we find a remarkably high degree of gene expression and regulatory network variation across even a shallow level of evolutionary divergence. We use these data to test hypotheses about the roles of specific trans-factors and cis-regulatory elements, how these roles may vary across venom genes and gene families, and how variation in regulatory systems drive diversity in venom phenotypes. Our results illustrate that differences in chromatin and genotype at regulatory elements play major roles in modulating expression. However, we also find that enhancer deletions, differences in transcription factor expression, and variation in activity of the insulator protein CTCF also likely impact venom phenotypes. Our findings provide insight into the diversity and gene-specificity of gene regulatory features and highlight the value of comparative studies to link gene regulatory network variation to phenotypic variation.

The whole genome sequence of Polish vaccine strain Mycobacterium bovis BCG Moreau

Currently, tuberculosis immunoprophylaxis is based solely on Bacillus Calmette-Guérin (BCG) vaccination, and some of the new potential tuberculosis vaccines are based on the BCG genome. Therefore, it is reasonable to analyze the genomes of individual BCG substrains. The aim of this study was the genetic characterization of the BCG-Moreau Polish (PL) strain used for the production of the BCG vaccine in Poland since 1955. Sequencing of different BCG lots showed that the strain was stable over a period of 59 years. As a result of comparison, BCG-Moreau PL with BCG-Moreau Rio de Janeiro (RDJ) 143 single nucleotide polymorphisms (SNPs) and 32 insertion/deletion mutations (INDELs) were identified. However, the verification of these mutations showed that the most significant were accumulated in the BCG-Moreau RDJ genome. The mutations unique to the Polish strain genome are 1 SNP and 2 INDEL. The strategy of combining short-read sequencing with long-read sequencing is currently the most optimal approach for sequencing bacterial genomes. With this approach, the only available genomic sequence of BCG-Moreau PL was obtained. This sequence will primarily be a reference point in the genetic control of the stability of the vaccine strain in the future. The results enrich knowledge about the microevolution and attenuation of the BCG vaccine substrains.

IMPORTANCE

The whole genome sequence obtained is the only genomic sequence of the strain that has been used for vaccine production in Poland since 1955. Sequencing of different BCG lots showed that the strain was stable over a period of 59 years. The comprehensive genomic analysis performed not only enriches knowledge about the microevolution and attenuation of the BCG vaccine substrains but also enables the utilization of identified markers as a reference point in the genetic control and identity tests of the stability of the vaccine strain in the future.

Genetic diversity and population structure in five Inner Mongolia cashmere goat populations using whole-genome genotyping

OBJECTIVE

As a charismatic species, cashmere goats have rich genetic resources. In the Inner Mongolia Autonomous Region, there are three cashmere goat varieties named and approved by the state. These goats are renowned for their high cashmere production and superior cashmere quality. Therefore, it is vitally important to protect their genetic resources as they will serve as breeding material for developing new varieties in the future.

METHODS

Three breeds including Inner Mongolia cashmere goats (IMCG), Hanshan White cashmere goats (HS), and Ujimqin white cashmere goats (WZMQ) were studied. IMCG were of three types: Aerbas (AEBS), Erlangshan (ELS), and Alashan (ALS). Nine DNA samples were collected for each population, and they were genomically re-sequenced to obtain high-depth data. The genetic diversity parameters of each population were estimated to determine selection intensity. Principal component analysis, phylogenetic tree construction and genetic differentiation parameter estimation were performed to determine genetic relationships among populations.

RESULTS

Samples from the 45 individuals from the five goat populations were sequenced, and 30,601,671 raw single nucleotide polymorphisms (SNPs) obtained. Then, variant calling was conducted using the reference genome, and 17,214,526 SNPs were retained after quality control. Individual sequencing depth of individuals ranged from 21.13× to 46.18×, with an average of 28.5×. In the AEBS, locus polymorphism (79.28) and expected heterozygosity (0.2554) proportions were the lowest, and the homologous consistency ratio (0.1021) and average inbreeding coefficient (0.1348) were the highest, indicating that this population had strong selection intensity. Conversely, ALS and WZMQ selection intensity was relatively low. Genetic distance between HS and the other four populations was relatively high, and genetic exchange existed among the other four populations.

CONCLUSION

The Inner Mongolia cashmere goat (AEBS type) population has a relatively high selection intensity and a low genetic diversity. The IMCG (ALS type) and WZMQ populations had relatively low selection intensity and high genetic diversity. The genetic distance between HS and the other four populations was relatively high, with a moderate degree of differentiation. Overall, these genetic variations provide a solid foundation for resource identification of Inner Mongolia Autonomous Region cashmere goats in the future.

GraphSlimmer: Preserving Read Mappability with the Minimum Number of Variants

Modern genomic datasets, like those generated under the 1000 Genome Project, contain millions of variants belonging to known haplotypes. Although these datasets are more representative than a single reference sequence and can alleviate issues like reference bias, they are significantly more computationally burdensome to work with, often involving large-indexed genome graph data structures for tasks such as read mapping. The construction, preprocessing, and mapping algorithms can require substantial computational resources depending on the size of these variant sets. Moreover, the accuracy of mapping algorithms has been shown to decrease when working with complete variant sets. Therefore, a drastically reduced set of variants that preserves important properties of the original set is desirable. This work provides a technique for finding a minimal subset of variants S such that for given parameters α and δ, all substrings up to length α in the haplotypes are guaranteed to be still alignable to the appropriate locations with either Hamming or edit distance at most δ, using only S . Our contributions include showing the NP-hardness and inapproximability of these optimization problems and providing Integer Linear Programming (ILP) formulations. Our edit distance ILP formulation carefully decomposes the problem according to variant locations, which allows it to scale to support all of chromosome 22's variants from the 1000 Genome Project. Our experiments also demonstrate a significant reduction in the number of variants. For example, for moderately long reads, e.g., α = 1000, over 75% of the variants can be removed while preserving read mappability with edit distance at most one.

The de novo, chromosome-level genome assembly of the sweet chestnut (Castanea sativa Mill.) Cv. Marrone Di Chiusa Pesio

OBJECTIVES

The sweet chestnut Castanea sativa Mill. is the only native Castanea species in Europe, and it is a tree of high economic value that provides appreciated fruits and valuable wood. In this study, we assembled a high-quality nuclear genome of the ancient Italian chestnut variety 'Marrone di Chiusa Pesio' using a combination of Oxford Nanopore Technologies long reads, whole-genome and Omni-C Illumina short reads.

DATA DESCRIPTION

The genome was assembled into 238 scaffolds with an N50 size of 21.8 Mb and an N80 size of 7.1 Mb for a total assembled sequence of 750 Mb. The BUSCO assessment revealed that 98.6% of the genome matched the embryophyte dataset, highlighting good completeness of the genetic space. After chromosome-level scaffolding, 12 chromosomes with a total length of 715.8 and 713.0 Mb were constructed for haplotype 1 and haplotype 2, respectively. The repetitive elements represented 37.3% and 37.4% of the total assembled genome in haplotype 1 and haplotype 2, respectively. A total of 57,653 and 58,146 genes were predicted in the two haplotypes, and approximately 73% of the genes were functionally annotated using the EggNOG-mapper. The assembled genome will be a valuable resource and reference for future chestnut breeding and genetic improvement.

De novo TLK1 and MDM1 mutations in a patient with a neurodevelopmental disorder and immunodeficiency

The Tousled-like kinases 1 and 2 (TLK1/TLK2) regulate DNA replication, repair and chromatin maintenance. TLK2 variants underlie the neurodevelopmental disorder (NDD) 'Intellectual Disability, Autosomal Dominant 57' (MRD57), characterized by intellectual disability and microcephaly. Several TLK1 variants have been reported in NDDs but their functional significance is unknown. A male patient presenting with ID, seizures, global developmental delay, hypothyroidism, and primary immunodeficiency was determined to have a heterozygous TLK1 variant (c.1435C>G, p.Q479E), as well as a mutation in MDM1 (c.1197dupT, p.K400∗). Cells expressing TLK1 p.Q479E exhibited reduced cytokine responses and elevated DNA damage, but not increased radiation sensitivity or DNA repair defects. The TLK1 p.Q479E variant impaired kinase activity but not proximal protein interactions. Our study provides the first functional characterization of NDD-associated TLK1 variants and suggests that, such as TLK2, TLK1 variants may impact development in multiple tissues and should be considered in the diagnosis of rare NDDs.

High heteroplasmy is associated with low mitochondrial copy number and selection against non-synonymous mutations in the snail Cepaea nemoralis

Molluscan mitochondrial genomes are unusual because they show wide variation in size, radical genome rearrangements and frequently show high variation (> 10%) within species. As progress in understanding this variation has been limited, we used whole genome sequencing of a six-generation matriline of the terrestrial snail Cepaea nemoralis, as well as whole genome sequences from wild-collected C. nemoralis, the sister species C. hortensis, and multiple other snail species to explore the origins of mitochondrial DNA (mtDNA) variation. The main finding is that a high rate of SNP heteroplasmy in somatic tissue was negatively correlated with mtDNA copy number in both Cepaea species. In individuals with under ten mtDNA copies per nuclear genome, more than 10% of all positions were heteroplasmic, with evidence for transmission of this heteroplasmy through the germline. Further analyses showed evidence for purifying selection acting on non-synonymous mutations, even at low frequency of the rare allele, especially in cytochrome oxidase subunit 1 and cytochrome b. The mtDNA of some individuals of Cepaea nemoralis contained a length heteroplasmy, including up to 12 direct repeat copies of tRNA-Val, with 24 copies in another snail, Candidula rugosiuscula, and repeats of tRNA-Thr in C. hortensis. These repeats likely arise due to error prone replication but are not correlated with mitochondrial copy number in C. nemoralis. Overall, the findings provide key insights into mechanisms of replication, mutation and evolution in molluscan mtDNA, and so will inform wider studies on the biology and evolution of mtDNA across animal phyla.

Detection of rare variants among nuclei populating the arbuscular mycorrhizal fungal model species Rhizophagus irregularis DAOM197198

Identifying genuine polymorphic variants is a significant challenge in sequence data analysis, although detecting low-frequency variants in sequence data is essential for estimating demographic parameters and investigating genetic processes, such as selection, within populations. Arbuscular mycorrhizal (AM) fungi are multinucleate organisms, in which individual nuclei collectively operate as a population, and the extent of genetic variation across nuclei has long been an area of scientific interest. In this study, we investigated the patterns of polymorphism discovery and the alternate allele frequency distribution by comparing polymorphism discovery in 2 distinct genomic sequence datasets of the AM fungus model species, Rhizophagus irregularis strain DAOM197198. The 2 datasets used in this study are publicly available and were generated either from pooled spores and hyphae or amplified single nuclei from a single spore. We also estimated the intraorganismal variation within the DAOM197198 strain. Our results showed that the 2 datasets exhibited different frequency patterns for discovered variants. The whole-organism dataset showed a distribution spanning low-, intermediate-, and high-frequency variants, whereas the single-nucleus dataset predominantly featured low-frequency variants with smaller proportions in intermediate and high frequencies. Furthermore, single nucleotide polymorphism density estimates within both the whole organism and individual nuclei confirmed the low intraorganismal variation of the DAOM197198 strain and that most variants are rare. Our study highlights the methodological challenges associated with detecting low-frequency variants in AM fungal whole-genome sequence data and demonstrates that alternate alleles can be reliably identified in single nuclei of AM fungi.

Comparative genomics analysis reveals sequence characteristics potentially related to host preference in Cryptosporidium xiaoi

Cryptosporidium spp. are important diarrhea-associated pathogens in humans and livestock. Among the known species, Cryptosporidium xiaoi, which causes cryptosporidiosis in sheep and goats, was previously recognized as a genotype of the bovine-specific Cryptosporidium bovis based on their high sequence identity in the ssrRNA gene. However, the lack of genomic data has limited characterization of the genetic differences between the two closely related species. In this study, we sequenced the genomes of two C. xiaoi isolates and performed comparative genomic analysis to identify the sequence uniqueness of this ovine-adapted species compared with other Cryptosporidium spp. Our results showed that C. xiaoi is genetically related to C. bovis as shown by their 95.8% genomic identity and similar gene content. Consistent with this, both C. xiaoi and C. bovis appear to have fewer genes encoding mitochondrial metabolic enzymes and invasion-related protein families. However, they appear to possess several species-specific genes. Further analysis indicates that the sequence differences between these two Cryptosporidium spp. are mainly in 24 highly polymorphic genes, half of which are located in the subtelomeric regions. Some of these subtelomeric genes encode secretory proteins that have undergone positive selection. In addition, the genomes of two C. xiaoi isolates, identified as subtypes XXIIIf and XXIIIh, share 99.9% nucleotide sequence identity, with six highly divergent genes encoding putative secretory proteins. Therefore, these species-specific genes and sequence polymorphism in subtelomeric genes probably contribute to the different host preference of C. xiaoi and C. bovis.

Integrating Multiple Database Resources to Elucidate the Gene Flow in Southeast Asian Pig Populations

The domestic pig (Sus scrofa) and its subfamilies have experienced long-term and extensive gene flow, particularly in Southeast Asia. Here, we analyzed 236 pigs, focusing on Yunnan indigenous, European commercial, East Asian, and Southeast Asian breeds, using the Pig Genomics Reference Panel (PGRP v1) of Pig Genotype-Tissue Expression (PigGTEx) to investigate gene flow and associated complex traits by integrating multiple database resources. In this study, we discovered evidence of admixtures from European pigs into the genome of Yunnan indigenous pigs. Additionally, we hypothesized that a potential conceptual gene flow route that may have contributed to the genetic composition of the Diannan small-ear pig is a gene exchange from the Vietnamese pig. Based on the most stringent gene introgression scan using the fd statistic, we identified three specific loci on chromosome 8, ranging from 51.65 to 52.45 Mb, which exhibited strong signatures of selection and harbored the NAF1, NPY1R, and NPY5R genes. These genes are associated with complex traits, such as fat mass, immunity, and litter weight, in pigs, as supported by multiple bio-functionalization databases. We utilized multiple databases to explore the potential dynamics of genetic exchange in Southeast Asian pig populations and elucidated specific gene functionalities.

The Landscape of Presence/Absence Variations during the Improvement of Rice

Rice is one of the most important staple crops in the world; therefore, the improvement of rice holds great significance for enhancing agricultural production and addressing food security challenges. Although there have been numerous studies on the role of single-nucleotide polymorphisms (SNPs) in rice improvement with the development of next-generation sequencing technologies, research on the role of presence/absence variations (PAVs) in the improvement of rice is limited. In particular, there is a scarcity of studies exploring the traits and genes that may be affected by PAVs in rice. Here, we extracted PAVs utilizing resequencing data from 148 improved rice varieties distributed in Asia. We detected a total of 33,220 PAVs and found that the number of variations decreased gradually as the length of the PAVs increased. The number of PAVs was the highest on chromosome 1. Furthermore, we identified a 6 Mb hotspot region on chromosome 11 containing 1091 PAVs in which there were 29 genes related to defense responses. By conducting a genome-wide association study (GWAS) using PAV variation data and phenotypic data for five traits (flowering time, plant height, flag leaf length, flag leaf width, and panicle number) across all materials, we identified 186 significantly associated PAVs involving 20 cloned genes. A haplotype analysis and expression analysis of candidate genes revealed that important genes might be affected by PAVs, such as the flowering time gene OsSFL1 and the flag leaf width gene NAL1. Our work investigated the pattern in PAVs and explored important PAV key functional genes associated with agronomic traits. Consequently, these results provide potential and exploitable genetic resources for rice breeding.

Genome-wide identification analysis in wild-type Solanum pinnatisectum reveals some genes defending against Phytophthora infestans

Solanum pinnatisectum exhibits strong resistance to late blight caused by Phytophthora infestans but only an incomplete genome assembly based on short Illumina reads has been published. In this study, we generated the first chromosome-level draft genome for the wild-type potato species S. pinnatisectum in China using Oxford Nanopore technology sequencing and Hi-C technology. The high-quality assembled genome size is 664 Mb with a scaffold N50 value of 49.17 Mb, of which 65.87% was occupied by repetitive sequences, and predominant long terminal repeats (42.51% of the entire genome). The genome of S. pinnatisectum was predicted to contain 34,245 genes, of which 99.34% were functionally annotated. Moreover, 303 NBS-coding disease resistance (R) genes were predicted in the S. pinnatisectum genome to investigate the potential mechanisms of resistance to late blight disease. The high-quality chromosome-level reference genome of S. pinnatisectum is expected to provide potential valuable resources for intensively and effectively investigating molecular breeding and genetic research in the future.

Genetic Causes and Genomic Consequences of Breakdown of Distyly in Linum trigynum

Distyly is an iconic floral polymorphism governed by a supergene, which promotes efficient pollen transfer and outcrossing through reciprocal differences in the position of sexual organs in flowers, often coupled with heteromorphic self-incompatibility. Distyly has evolved convergently in multiple flowering plant lineages, but has also broken down repeatedly, often resulting in homostylous, self-compatible populations with elevated rates of self-fertilization. Here, we aimed to study the genetic causes and genomic consequences of the shift to homostyly in Linum trigynum, which is closely related to distylous Linum tenue. Building on a high-quality genome assembly, we show that L. trigynum harbors a genomic region homologous to the dominant haplotype of the distyly supergene conferring long stamens and short styles in L. tenue, suggesting that loss of distyly first occurred in a short-styled individual. In contrast to homostylous Primula and Fagopyrum, L. trigynum harbors no fixed loss-of-function mutations in coding sequences of S-linked distyly candidate genes. Instead, floral gene expression analyses and controlled crosses suggest that mutations downregulating the S-linked LtWDR-44 candidate gene for male self-incompatibility and/or anther height could underlie homostyly and self-compatibility in L. trigynum. Population genomic analyses of 224 whole-genome sequences further demonstrate that L. trigynum is highly self-fertilizing, exhibits significantly lower genetic diversity genome-wide, and is experiencing relaxed purifying selection and less frequent positive selection on nonsynonymous mutations relative to L. tenue. Our analyses shed light on the loss of distyly in L. trigynum, and advance our understanding of a common evolutionary transition in flowering plants.

Biallelic NAA60 variants with impaired n-terminal acetylation capacity cause autosomal recessive primary familial brain calcifications

Primary familial brain calcification (PFBC) is characterized by calcium deposition in the brain, causing progressive movement disorders, psychiatric symptoms, and cognitive decline. PFBC is a heterogeneous disorder currently linked to variants in six different genes, but most patients remain genetically undiagnosed. Here, we identify biallelic NAA60 variants in ten individuals from seven families with autosomal recessive PFBC. The NAA60 variants lead to loss-of-function with lack of protein N-terminal (Nt)-acetylation activity. We show that the phosphate importer SLC20A2 is a substrate of NAA60 in vitro. In cells, loss of NAA60 caused reduced surface levels of SLC20A2 and a reduction in extracellular phosphate uptake. This study establishes NAA60 as a causal gene for PFBC, provides a possible biochemical explanation of its disease-causing mechanisms and underscores NAA60-mediated Nt-acetylation of transmembrane proteins as a fundamental process for healthy neurobiological functioning.

Chromosome-level genome assembly and annotation of the Yunling cattle with PacBio and Hi-C sequencing data

Yunling cattle is a new breed of beef cattle bred in Yunnan Province, China. It is bred by crossing the Brahman, the Murray Grey and the Yunnan Yellow cattle. Yunling cattle can adapt to the tropical and subtropical climate environment, and has good reproductive ability and growth speed under high temperature and high humidity conditions, it also has strong resistance to internal and external parasites and with good beef performance. In this study, we generated a high-quality chromosome-level genome assembly of a male Yunling cattle using a combination of short reads sequencing, PacBio HiFi sequencing and Hi-C scaffolding technologies. The genome assembly(3.09 Gb) is anchored to 31 chromosomes(29 autosomes plus one X and Y), with a contig N50 of 35.97 Mb and a scaffold N50 of 112.01 Mb. It contains 1.62 Gb of repetitive sequences and 20,660 protein-coding genes. This first construction of the Yunling cattle genome provides a valuable genetic resource that will facilitate further study of the genetic diversity of bovine species and accelerate Yunling cattle breeding efforts.

Population genomics of Agrotis segetum provide insights into the local adaptive evolution of agricultural pests

BACKGROUND

The adaptive mechanisms of agricultural pests are the key to understanding the evolution of the pests and to developing new control strategies. However, there are few studies on the genetic basis of adaptations of agricultural pests. The turnip moth, Agrotis segetum (Lepidoptera: Noctuidae) is an important underground pest that affects a wide range of host plants and has a strong capacity to adapt to new environments. It is thus a good model for studying the adaptive evolution of pest species.

RESULTS

We assembled a high-quality reference genome of A. segetum using PacBio reads. Then, we constructed a variation map of A. segetum by resequencing 98 individuals collected from six natural populations in China. The analysis of the population structure showed that all individuals were divided into four well-differentiated populations, corresponding to their geographical distribution. Selective sweep analysis and environmental association studies showed that candidate genes associated with local adaptation were functionally correlated with detoxification metabolism and glucose metabolism.

CONCLUSIONS

Our study of A. segetum has provided insights into the genetic mechanisms of local adaptation and evolution; it has also produced genetic resources for developing new pest management strategies.

Karyotypic stasis and swarming influenced the evolution of viral tolerance in a species-rich bat radiation

The emergence of COVID-19 and severe acute respiratory syndrome (SARS) has prioritized understanding bats' viral tolerance. Myotis bats are exceptionally species rich and have evolved viral tolerance. They also exhibit swarming, a cryptic behavior where large, multi-species assemblages gather for mating, which has been hypothesized to promote interspecific hybridization. To resolve the coevolution of genome architecture and their unusual antiviral tolerance, we undertook a phylogenomic analysis of 60 Old World Myotis genomes. We demonstrate an extensive history of introgressive hybridization that has replaced the species phylogeny across 17%-93% of the genome except for pericentromeric regions of macrochromosomes. Introgression tracts were enriched on microchromosome regions containing key antiviral pathway genes overexpressed during viral challenge experiments. Together, these results suggest that the unusual Myotis karyotype may have evolved to selectively position immune-related genes in high recombining genomic regions prone to introgression of divergent alleles, including a diversity of interleukin loci responsible for the release of pro-inflammatory cytokines.

Oxidative stress resistance prompts pyrroloquinoline quinone biosynthesis in Hyphomicrobium denitrificans H4-45

Pyrroloquinoline quinone (PQQ) is a natural antioxidant with diverse applications in food and pharmaceutical industries. A lot of effort has been devoted toward the discovery of PQQ high-producing microbial species and characterization of biosynthesis, but it is still challenging to achieve a high PQQ yield. In this study, a combined strategy of random mutagenesis and adaptive laboratory evolution (ALE) with fermentation optimization was applied to improve PQQ production in Hyphomicrobium denitrificans H4-45. A mutant strain AE-9 was obtained after nearly 400 generations of UV-LiCl mutagenesis, followed by an ALE process, which was conducted with a consecutive increase of oxidative stress generated by kanamycin, sodium sulfide, and potassium tellurite. In the flask culture condition, the PQQ production in mutant strain AE-9 had an 80.4% increase, and the cell density increased by 14.9% when compared with that of the initial strain H4-45. Moreover, batch and fed-batch fermentation processes were optimized to further improve PQQ production by pH control strategy, methanol and H2O2 feed flow, and segmented fermentation process. Finally, the highest PQQ production and productivity of the mutant strain AE-9 reached 307 mg/L and 4.26 mg/L/h in a 3.7-L bioreactor, respectively. Whole genome sequencing analysis showed that genetic mutations in the ftfL gene and thiC gene might contribute to improving PQQ production by enhancing methanol consumption and cell growth in the AE-9 strain. Our study provided a systematic strategy to obtain a PQQ high-producing mutant strain and achieve high production of PQQ in fermentation. These practical methods could be applicable to improve the production of other antioxidant compounds with uncleared regulation mechanisms. KEY POINTS: • Improvement of PQQ production by UV-LiCl mutagenesis combined with adaptive laboratory evolution (ALE) and fermentation optimization. • A consecutive increase of oxidative stress could be used as the antagonistic factor for ALE to enhance PQQ production. • Mutations in the ftfL gene and thiC gene indicated that PQQ production might be increased by enhancing methanol consumption and cell growth.

Rapid diversification of grey mangroves (Avicennia marina) driven by geographic isolation and extreme environmental conditions in the Arabian Peninsula

Biological systems occurring in ecologically heterogeneous and spatially discontinuous habitats provide an ideal opportunity to investigate the relative roles of neutral and selective factors in driving lineage diversification. The grey mangroves (Avicennia marina) of Arabia occur at the northern edge of the species' range and are subject to variable, often extreme, environmental conditions, as well as historic large fluctuations in habitat availability and connectivity resulting from Quaternary glacial cycles. Here, we analyse fully sequenced genomes sampled from 19 locations across the Red Sea, the Arabian Sea and the Persian/Arabian Gulf (PAG) to reconstruct the evolutionary history of the species in the region and to identify adaptive mechanisms of lineage diversification. Population structure and phylogenetic analyses revealed marked genetic structure correlating with geographic distance and highly supported clades among and within the seas surrounding the Arabian Peninsula. Demographic modelling showed times of divergence consistent with recent periods of geographic isolation and low marine connectivity during glaciations, suggesting the presence of (cryptic) glacial refugia in the Red Sea and the PAG. Significant migration was detected within the Red Sea and the PAG, and across the Strait of Hormuz to the Arabian Sea, suggesting gene flow upon secondary contact among populations. Genetic-environment association analyses revealed high levels of adaptive divergence and detected signs of multi-loci local adaptation driven by temperature extremes and hypersalinity. These results support a process of rapid diversification resulting from the combined effects of historical factors and ecological selection and reveal mangrove peripheral environments as relevant drivers of lineage diversity.

The Addis Ababa Lions: Whole-Genome Sequencing of a Rare and Precious Population

Lions are widely known as charismatic predators that once roamed across the globe, but their populations have been greatly affected by environmental factors and human activities over the last 150 yr. Of particular interest is the Addis Ababa lion population, which has been maintained in captivity at around 20 individuals for over 75 yr, while many wild African lion populations have become extinct. In order to understand the molecular features of this unique population, we conducted a whole-genome sequencing study on 15 Addis Ababa lions and detected 4.5 million distinct genomic variants compared with the reference African lion genome. Using functional annotation, we identified several genes with mutations that potentially impact various traits such as mane color, body size, reproduction, gastrointestinal functions, cardiovascular processes, and sensory perception. These findings offer valuable insights into the genetics of this threatened lion population.

HiPhase: jointly phasing small, structural, and tandem repeat variants from HiFi sequencing

MOTIVATION

In diploid organisms, phasing is the problem of assigning the alleles at heterozygous variants to one of two haplotypes. Reads from PacBio HiFi sequencing provide long, accurate observations that can be used as the basis for both calling and phasing variants. HiFi reads also excel at calling larger classes of variation, such as structural or tandem repeat variants. However, current phasing tools typically only phase small variants, leaving larger variants unphased.

RESULTS

We developed HiPhase, a tool that jointly phases SNVs, indels, structural, and tandem repeat variants. The main benefits of HiPhase are (i) dual mode allele assignment for detecting large variants, (ii) a novel application of the A*-algorithm to phasing, and (iii) logic allowing phase blocks to span breaks caused by alignment issues around reference gaps and homozygous deletions. In our assessment, HiPhase produced an average phase block NG50 of 480 kb with 929 switchflip errors and fully phased 93.8% of genes, improving over the current state of the art. Additionally, HiPhase jointly phases SNVs, indels, structural, and tandem repeat variants and includes innate multi-threading, statistics gathering, and concurrent phased alignment output generation.

AVAILABILITY AND IMPLEMENTATION

HiPhase is available as source code and a pre-compiled Linux binary with a user guide at https://github.com/PacificBiosciences/HiPhase.

Chromosome-level genome assembly and annotation of the yellow grouper, Epinephelus awoara

Epinephelus awoara, as known as yellow grouper, is a significant economic marine fish that has been bred artificially in China. However, the genetic structure and evolutionary history of yellow grouper remains largely unknown. Here, this work presents the high-quality chromosome-level genome assembly of yellow grouper using PacBio single molecule sequencing technique (SMRT) and High-through chromosome conformation capture (Hi-C) technologies. The 984.48 Mb chromosome-level genome of yellow grouper was assembled, with a contig N50 length of 39.77 Mb and scaffold N50 length of 41.39 Mb. Approximately 99.76% of assembled sequences were anchored into 24 pseudo-chromosomes with the assistance of Hi-C reads. Furthermore, approximately 41.17% of the genome was composed of repetitive elements. In total, 24,541 protein-coding genes were predicted, of which 22,509 (91.72%) genes were functionally annotated. The highly accurate, chromosome-level reference genome assembly and annotation are crucial to the understanding of population genetic structure, adaptive evolution and speciation of the yellow grouper.

Monoallelically expressed noncoding RNAs form nucleolar territories on NOR-containing chromosomes and regulate rRNA expression

Out of the several hundred copies of rRNA genes arranged in the nucleolar organizing regions (NOR) of the five human acrocentric chromosomes, ~50% remain transcriptionally inactive. NOR-associated sequences and epigenetic modifications contribute to the differential expression of rRNAs. However, the mechanism(s) controlling the dosage of active versus inactive rRNA genes within each NOR in mammals is yet to be determined. We have discovered a family of ncRNAs, SNULs (Single NUcleolus Localized RNA), which form constrained sub-nucleolar territories on individual NORs and influence rRNA expression. Individual members of the SNULs monoallelically associate with specific NOR-containing chromosomes. SNULs share sequence similarity to pre-rRNA and localize in the sub-nucleolar compartment with pre-rRNA. Finally, SNULs control rRNA expression by influencing pre-rRNA sorting to the DFC compartment and pre-rRNA processing. Our study discovered a novel class of ncRNAs influencing rRNA expression by forming constrained nucleolar territories on individual NORs.

Assessing genome-wide adaptations associated with range expansion in the pink rice borer, Sesamia inferens

Understanding the genetic basis of adaptive evolution following habitat expansion can have important implications for pest management. The pink rice borer (PRB), Sesamia inferens (Walker), is a destructive pest of rice that was historically restricted to regions south of 34° N latitude in China. However, with changes in global climate and farming practices, the distribution of this moth has progressively expanded, encompassing most regions in North China. Here, 3 highly differentiated subpopulations were discovered using high-quality single-nucleotide polymorphism and structural variant datasets across China, corresponding to northern, southern China regions, and the Yunnan-Guizhou Plateau, with significant patterns of isolation by geographic and environmental distances. Our estimates of evolutionary history indicate asymmetric migration with varying population sizes across the 3 subpopulations. Selective sweep analyses estimated strong selection at insect cuticle glycine-rich cuticular protein genes which are associated with enhanced desiccation adaptability in the northern group, and at the histone-lysine-N-methyltransferase gene associated with range expansion and local adaptation in the Shandong population. Our findings have significant implications for the development of effective strategies to control this pest.

CaeNDR, the Caenorhabditis Natural Diversity Resource

Studies of model organisms have provided important insights into how natural genetic differences shape trait variation. These discoveries are driven by the growing availability of genomes and the expansive experimental toolkits afforded to researchers using these species. For example, Caenorhabditis elegans is increasingly being used to identify and measure the effects of natural genetic variants on traits using quantitative genetics. Since 2016, the C. elegans Natural Diversity Resource (CeNDR) has facilitated many of these studies by providing an archive of wild strains, genome-wide sequence and variant data for each strain, and a genome-wide association (GWA) mapping portal for the C. elegans community. Here, we present an updated platform, the Caenorhabditis Natural Diversity Resource (CaeNDR), that enables quantitative genetics and genomics studies across the three Caenorhabditis species: C. elegans, C. briggsae and C. tropicalis. The CaeNDR platform hosts several databases that are continually updated by the addition of new strains, whole-genome sequence data and annotated variants. Additionally, CaeNDR provides new interactive tools to explore natural variation and enable GWA mappings. All CaeNDR data and tools are accessible through a freely available web portal located at caendr.org.

BCFtools/liftover: an accurate and comprehensive tool to convert genetic variants across genome assemblies

MOTIVATION

Many genetics studies report results tied to genomic coordinates of a legacy genome assembly. However, as assemblies are updated and improved, researchers are faced with either realigning raw sequence data using the updated coordinate system or converting legacy datasets to the updated coordinate system to be able to combine results with newer datasets. Currently available tools to perform the conversion of genetic variants have numerous shortcomings, including poor support for indels and multi-allelic variants, that lead to a higher rate of variants being dropped or incorrectly converted. As a result, many researchers continue to work with and publish using legacy genomic coordinates.

RESULTS

Here we present BCFtools/liftover, a tool to convert genomic coordinates across genome assemblies for variants encoded in the variant call format with improved support for indels represented by different reference alleles across genome assemblies and full support for multi-allelic variants. It further supports variant annotation fields updates whenever the reference allele changes across genome assemblies. The tool has the lowest rate of variants being dropped with an order of magnitude less indels dropped or incorrectly converted and is an order of magnitude faster than other tools typically used for the same task. It is particularly suited for converting variant callsets from large cohorts to novel telomere-to-telomere assemblies as well as summary statistics from genome-wide association studies tied to legacy genome assemblies.

AVAILABILITY AND IMPLEMENTATION

The tool is written in C and freely available under the MIT open source license as a BCFtools plugin available at http://github.com/freeseek/score.

A splicing site change between exon 5 and 6 of the nuclear-encoded chloroplast-localized HvYGL8 gene results in reduced chlorophyll content and plant height in barley

The chloroplast is an important cellular organelle and metabolic hub, which is not only responsible for plant photosynthesis but is also involved in the de novo biosynthesis of pigments, fatty acids, and hormone metabolisms. Several genes that are responsible for rice leaf color variations have been reported to be directly or indirectly involved in chlorophyll biosynthesis and chloroplast development, whereas a few genes have been functionally confirmed to be responsible for leaf color changes in barley at the molecular level. In this study, we obtained a yellow leaf and dwarf ygl8 mutant from the progeny of Morex (a variety of barley) seeds treated with EMS. We performed bulked-segregant analysis (BSA) and RNA-seq analysis and targeted a UMP kinase encoding gene, YGL8, which generated a splicing site change between exon 5 and 6 of YGL8 due to a G to A single-nucleotide transition in the 5th exon/intron junction in the ygl8 mutant. The splicing site change between exon 5 and 6 of YGL8 had no effects on chloroplast subcellular localization but resulted in an additional loop in the UMP kinase domain, which might disturb the access of the substrates. On one hand, the splicing site change between exon 5 and 6 of YGL8 downregulated the transcriptional expression of chloroplast-encoded genes and chlorophyll-biosynthesis-related genes in a temperature-dependent manner in the ygl8 mutant. On the other hand, the downregulation of bioactive GA-biosynthesis-related GA20ox genes and cell-wall-cellulose-biosynthesis-related CesA genes was also observed in the ygl8 mutant, which led to a reduction in plant height. Our study will facilitate the understanding of the regulation of leaf color and plant height in barley.

The loss of neoantigens is an important reason for immune escape in multiple myeloma patients with high intratumor heterogeneity

OBJECTIVES

Intratumor heterogeneity (ITH) is an important factor for clinical outcomes in patients with multiple myeloma (MM). High ITH has been proven to be a key reason for tumor immune escape and treatment resistance. Neoantigens are thought to be associated with ITH, but the specific correlation and functional basis for this remains unclear.

METHODS

We study this question through the whole-exome sequencing (WES) data from 43 high ITH newly diagnosed MM patients in our center. Mutant allele tumor heterogeneity (MATH) was conducted to quantify ITH. The cutoff value for high intratumor heterogeneity was determined by comparing MATH of different kinds of tumors. NeoPredPipe was performed to predict neoantigens and binding affinity.

RESULTS

Compared to other tumors, MM has a relatively low tumor mutation burden but a high ITH. Patients with high MATH had significantly shorter progression-free survival times than those with low MATH (p = 0.001). In high ITH samples, there is a decrease in strong-binding neoantigens (p = 0.019). The loss of strong-binding neoantigens is a key factor for insensitivity to therapy (p = 0.015). Loss of heterozygosity in HLA was not observed. In addition, patients with fewer neoantigens loss had higher rates of disease remission (p = 0.047). CD8 + T cells (p = 0.012) and NK cells (p = 0.011) decreased significantly in patients with high neoantigens loss rate. A prediction model based on neoantigens was built to evaluate the strength of immune escape.

CONCLUSION

The loss of strong-binding neoantigens explains why tumors with high ITH have a higher degree of immune escape and may be feasible for deciding the clinical treatment of MM.

EpiTyping: analysis of epigenetic aberrations in parental imprinting and X-chromosome inactivation using RNA-seq

Human pluripotent stem cells (hPSCs) hold a central role in studying human development, in disease modeling and in regenerative medicine. These cells not only acquire genetic modifications when kept in culture, but they may also harbor epigenetic aberrations, mainly involving parental imprinting and X-chromosome inactivation. Here we present a detailed bioinformatic protocol for detecting such aberrations using RNA sequencing data. We provide a pipeline designed to process and analyze RNA sequencing data for the identification of abnormal biallelic expression of imprinted genes, and thus detect loss of imprinting. Furthermore, we show how to differentiate among X-chromosome inactivation, full activation and aberrant erosion of X chromosome in female hPSCs. In addition to providing bioinformatic tools, we discuss the impact of such epigenetic variations in hPSCs on their utility for various purposes. This pipeline can be used by any user with basic understanding of the Linux command line. It is available on GitHub as a software container ( https://github.com/Gal-Keshet/EpiTyping ) and produces reliable results in 1-4 d.

Genetic and environmental drivers of migratory behavior in western burrowing owls and implications for conservation and management

Migration is driven by a combination of environmental and genetic factors, but many questions remain about those drivers. Potential interactions between genetic and environmental variants associated with different migratory phenotypes are rarely the focus of study. We pair low coverage whole genome resequencing with a de novo genome assembly to examine population structure, inbreeding, and the environmental factors associated with genetic differentiation between migratory and resident breeding phenotypes in a species of conservation concern, the western burrowing owl (Athene cunicularia hypugaea). Our analyses reveal a dichotomy in gene flow depending on whether the population is resident or migratory, with the former being genetically structured and the latter exhibiting no signs of structure. Among resident populations, we observed significantly higher genetic differentiation, significant isolation-by-distance, and significantly elevated inbreeding. Among migratory breeding groups, on the other hand, we observed lower genetic differentiation, no isolation-by-distance, and substantially lower inbreeding. Using genotype-environment association analysis, we find significant evidence for relationships between migratory phenotypes (i.e., migrant versus resident) and environmental variation associated with cold temperatures during the winter and barren, open habitats. In the regions of the genome most differentiated between migrants and residents, we find significant enrichment for genes associated with the metabolism of fats. This may be linked to the increased pressure on migrants to process and store fats more efficiently in preparation for and during migration. Our results provide a significant contribution toward understanding the evolution of migratory behavior and vital insight into ongoing conservation and management efforts for the western burrowing owl.

COVID-19 Genomic Surveillance in Bangui (Central African Republic) Reveals a Landscape of Circulating Variants Linked to Validated Antiviral Targets of SARS-CoV-2 Proteome

Since its outbreak, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) spread rapidly, causing the Coronavirus Disease 19 (COVID-19) pandemic. Even with the vaccines' administration, the virus continued to circulate due to inequal access to prevention and therapeutic measures in African countries. Information about COVID-19 in Africa has been limited and contradictory, and thus regional studies are important. On this premise, we conducted a genomic surveillance study about COVID-19 lineages circulating in Bangui, Central African Republic (CAR). We collected 2687 nasopharyngeal samples at four checkpoints in Bangui from 2 to 22 July 2021. Fifty-three samples tested positive for SARS-CoV-2, and viral genomes were sequenced to look for the presence of different viral strains. We performed phylogenetic analysis and described the lineage landscape of SARS-CoV-2 circulating in the CAR along 15 months of pandemics and in Africa during the study period, finding the Delta variant as the predominant Variant of Concern (VoC). The deduced aminoacidic sequences of structural and non-structural genes were determined and compared to reference and reported isolates from Africa. Despite the limited number of positive samples obtained, this study provides valuable information about COVID-19 evolution at the regional level and allows for a better understanding of SARS-CoV-2 circulation in the CAR.

The complete mitochondrial genome of the lowland paca (Cuniculus paca) and its phylogenetic relationship with other New World hystricognath rodents

The lowland paca (Cuniculus paca) is a nocturnal, widespread, and solitary large-sized rodent in the family Cuniculidae, and one of the most frequently hunted mammals in the Neotropical forests of Latin America. We assembled the first complete mitochondrial genome of lowland paca using three closely related hystricognath species as reference sequences. The mitochondrial genome is 16,770 basepairs (bp) in length, with similar characteristics of vertebrate mitochondrial genomes. We performed phylogenetic analyses using 26 mitochondrial genome of hystricognath species based on thirteen protein-coding genes. The result confirms the taxonomical placement among the New World hystricognath rodents with high support. The placement is consistent with previous phylogenetic studies based on individual mitochondrial and nuclear genes. The current study improves the phylogenic resolution of hystricognath rodents.

Comprehensive detection of structural variation and transposable element differences between wild type laboratory lineages of C. elegans

Genomic structural variations (SVs) and transposable elements (TEs) can be significant contributors to genome evolution, altered gene expression, and risk of genetic diseases. Recent advancements in long-read sequencing have greatly improved the quality of de novo genome assemblies and enhanced the detection of sequence variants at the scale of hundreds or thousands of bases. Comparisons between two diverged wild isolates of Caenorhabditis elegans, the Bristol and Hawaiian strains, have been widely utilized in the analysis of small genetic variations. Genetic drift, including SVs and rearrangements of repeated sequences such as TEs, can occur over time from long-term maintenance of wild type isolates within the laboratory. To comprehensively detect both large and small structural variations as well as TEs due to genetic drift, we generated de novo genome assemblies and annotations for each strain from our lab collection using both long- and short-read sequencing and compared our assemblies and annotations with that of other lab wild type strains. Within our lab assemblies, we annotate over 3.1Mb of sequence divergence between the Bristol and Hawaiian isolates: 337,584 SNPs, 94,503 small insertion-deletions (<50bp), and 4,334 structural variations (>50bp). Further, we define the location and movement of specific DNA TEs between N2 Bristol and CB4856 Hawaiian wild type isolates. Specifically, we find the N2 Bristol genome has 20.6% more TEs from the Tc1/mariner family than the CB4856 Hawaiian genome. Moreover, we identified Zator elements as the most abundant and mobile TE family in the genome. Using specific TE sequences with unique SNPs, we also identify 38 TEs that moved intrachromosomally and 9 TEs that moved interchromosomally between the N2 Bristol and CB4856 Hawaiian genomes. By comparing the de novo genome assembly of our lab collection Bristol isolate to the VC2010 Bristol assembly, we also reveal that lab lineages display over 2 Mb of total variation: 1,162 SNPs, 1,528 indels, and 897 SVs with 95% of the variation due to SVs. Overall, our work demonstrates the unique contribution of SVs and TEs to variation and genetic drift between wild type laboratory strains assumed to be isogenic despite growing evidence of genetic drift and phenotypic variation.

Production of leishmanin skin test antigen from Leishmania donovani for future reintroduction in the field

The leishmanin skin test was used for almost a century to detect exposure and immunity to Leishmania, the causative agent of leishmaniasis, a major neglected tropical disease. Due to a lack of antigen used for the intradermal injection, the leishmanin skin test is no longer available. As leishmaniasis control programs are advancing and new vaccines are entering clinical trials, it is essential to re-introduce the leishmanin skin test. Here we establish a Leishmania donovani strain and describe the production, under Good Laboratory Practice conditions, of leishmanin soluble antigen used to induce the leishmanin skin test in animal models of infection and vaccination. Using a mouse model of cutaneous leishmaniasis and a hamster model of visceral leishmaniasis, soluble antigen induces a leishmanin skin test response following infection and vaccination with live attenuated Leishmania major (LmCen-/-). Both the CD4+ and CD8+ T-cells are necessary for the leishmanin skin test response. This study demonstrates the feasibility of large-scale production of leishmanin antigen addressing a major bottleneck for performing the leishmanin skin test in future surveillance and vaccine clinical trials.

Investigation of phylogenetic relationships within Saxifraga diversifolia complex (Saxifragaceae) based on restriction-site associated DNA sequence markers

Subsect. Hirculoideae Engl. & Irmsch., belonging to Saxifraga sect. Ciliatae Haw., has high species richness. It can be divided into S. diversifolia, S. pseudohirculus, and S. sinomontana complexes based on morphological characteristics. The species with prominent leaf veins on the posterior leaf edge were placed in the S. diversifolia complex, which is mainly distributed on the eastern and southern margins of the Qinghai-Tibetan Plateau. In this study, 53 samples, representing 15 of the 33 described species in the S. diversifolia complex, were sequenced using the Restriction-site Associated DNA Sequence (RAD-seq) technique. A total of 111,938 high-quality SNP loci were screened to investigate the phylogenetic relationships within the S. diversifolia complex. The result of the neighbor-joining (NJ) tree shows that the S. diversifolia complex is a paraphyletic group. Despite of some inconsistencies as revealed by genetic structural analysis, clustering results of representative species reconstructed by both NJ and principal component analysis analyses support previous biogeographic and morphological evidences. In addition, long-distance gene flow events for 11 taxa were detected in the S. diversifolia complex, respectively from S. implicans 1 to S. implicans 2, S. diversifolia and S. maxionggouensis, and from S. maxionggouensis to S. nigroglandulifera. These findings may improve our comprehension of the phylogeny, classification, and evolution of the S. diversifolia complex.

Whole genome sequencing and annotation of Scleroderma yunnanense, the only edible Scleroderma species

Scleroderma yunnanense, an ectomycorrhizal fungus, is a popular edible mushroom within the Yunnan Province of Southwest China that holds great ecological and economic implications. However, despite its significance, there remains limited information about this species. Therefore, we sequenced S. yunnanense genome to identify the functional genes of S. yunnanense involved in secondary metabolite and carbohydrate production pathways. First, we present the 40.43 Mb high-quality reference genome for S. yunnanense, distributed across 35 contigs; moreover, the N50 contig size was found to reach 3.31 Mb and contained 8877 functional genes. Finally, genome annotation was conducted to compare the functional genes of S. yunnanense with protein sequences from different publicly available databases. Taken together, we identified 12 biosynthetic gene clusters across 10 contigs; among these were 13 key mevalonate (MVA) pathway enzymes, a key tyrosinase enzyme in the 3,4-dihydroxyphenylalanine (DOPA) pathway that is responsible for producing DOPA melanins, and 16 enzymes involved in uridine diphosphate glucose biosynthesis. Overall, this study presents the first genome assembly and annotation of S. yunnanense; ultimately, this information will be important in the elucidation of the biological activities and artificial domestication of this fungus.

Allele-specific expression reveals genetic drivers of tissue regeneration in mice

In adult mammals, skin wounds typically heal by scarring rather than through regeneration. In contrast, "super-healer" Murphy Roths Large (MRL) mice have the unusual ability to regenerate ear punch wounds; however, the molecular basis for this regeneration remains elusive. Here, in hybrid crosses between MRL and non-regenerating mice, we used allele-specific gene expression to identify cis-regulatory variation associated with ear regeneration. Analyzing three major cell populations (immune, fibroblast, and endothelial), we found that genes with cis-regulatory differences specifically in fibroblasts were associated with wound-healing pathways and also co-localized with quantitative trait loci for ear wound-healing. Ectopic treatment with one of these proteins, complement factor H (CFH), accelerated wound repair and induced regeneration in typically fibrotic wounds. Through single-cell RNA sequencing (RNA-seq), we observed that CFH treatment dramatically reduced immune cell recruitment to wounds, suggesting a potential mechanism for CFH's effect. Overall, our results provide insights into the molecular drivers of regeneration with potential clinical implications.

Self-sampling monkeypox virus testing in high-risk populations, asymptomatic or with unrecognized Mpox, in Spain

The recent monkeypox virus (MPXV) outbreak was of global concern and has mainly affected gay, bisexual and other men who have sex with men (GBMSM). Here we assess prevalence of MPXV in high-risk populations of GBMSM, trans women (TW) and non-binary people without symptoms or with unrecognized monkeypox (Mpox) symptoms, using a self-sampling strategy. Anal and pharyngeal swabs are tested by MPXV real-time PCR and positive samples are tested for cytopathic effect (CPE) in cell culture. 113 individuals participated in the study, 89 (78.76%) were cis men, 17 (15.04%) were TW. The median age was 35.0 years (IQR: 30.0-43.0), 96 (85.02%) individuals were gay or bisexual and 72 (63.72%) were migrants. Seven participants were MPXV positive (6.19% (95% CI: 1.75%-10.64%)). Five tested positive in pharyngeal swabs, one in anal swab and one in both. Six did not present symptoms recognized as MPXV infection. Three samples were positive for CPE, and showed anti-vaccinia pAb staining by FACS and confocal microscopy. This suggests that unrecognized Mpox cases can shed infectious virus. Restricting testing to individuals reporting Mpox symptoms may not be sufficient to contain outbreaks.

Evolution of the viral genomes of SARS-CoV-2 in association with the changes in local condition: a genomic epidemiological study of a suburban city of Japan

Understanding the factors driving the spread and evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) at the local, regional, national, and international levels is important in protecting against future pandemics. By exploring their viral genomes, we attempted to analyse the spread of SARS-CoV-2 and its evolutionary convergence in Kashiwa City, as an example of a representative commuter town in Japan. From September 2020 to January 2023, a total of 47,134 nasopharyngeal swab and saliva specimens were collected from patients in 47 local clinics and hospitals, covering the vast majority of healthcare facilities. All SARS-CoV-2-positive samples were subjected to whole genome sequencing. Based on the analysis of 5,536 identified genomes, all major strains were represented. Unique regional mutations were occasionally identified in each strain. Inspection of these mutations revealed that the overall base substitution rate increased with progressive waves of the pandemic, at an overall rate of 2.56 bases/year. Interestingly, the spread and evolutionary patterns appeared to be distinct between regions and between individual clinics. Further analysis of the synonymous base substitution rate showed that the speed of viral evolution accelerated coincident with the beginning of public vaccination. Comprehensive genomic epidemiological studies, as presented here, should be useful in precisely understanding the pandemic and preparing for possible future pandemics.