Genomic Regions Research Articles

Little evidence for homoeologous gene conversion and homoeologous exchange events in Gossypium allopolyploids.

A complicating factor in analyzing allopolyploid genomes is the possibility of physical interactions between homoeologous chromosomes during meiosis, resulting in either crossover (homoeologous exchanges) or non-crossover products (homoeologous gene conversion). Homoeologous gene conversion was first described in cotton by comparing SNP patterns in sequences from two diploid progenitors with those from the allopolyploid subgenomes. These analyses, however, did not explicitly consider other evolutionary scenarios that may give rise to similar SNP patterns as homoeologous gene conversion, creating uncertainties about the reality of the inferred gene conversion events. Here, we use an expanded phylogenetic sampling of high-quality genome assemblies from seven allopolyploid Gossypium species (all derived from the same polyploidy event), four diploid species (two closely related to each subgenome), and a diploid outgroup to derive a robust method for identifying potential genomic regions of gene conversion and homoeologous exchange. We found little evidence for homoeologous gene conversion in allopolyploid cottons, and that only two of the 40 best-supported events were shared by more than one species. We did, however, reveal a single, shared homoeologous exchange event at one end of chromosome 1, which occurred shortly after allopolyploidization but prior to divergence of the descendant species. Overall, our analyses demonstrated that homoeologous gene conversion and homoeologous exchanges are uncommon in Gossypium, affecting between zero and 24 genes per subgenome (0.0-0.065%) across the seven species. More generally, we highlighted the potential problems of using simple four-taxon tests to investigate patterns of homoeologous gene conversion in established allopolyploids.

Discovery of NFκB2-Coordinated Dual Regulation of Mitochondrial and Nuclear Genomes Leads to an Effective Therapy for Acute Myeloid Leukemia.

Acute myeloid leukemia (AML) has a poor survival rate for both pediatric and adult patients due to its frequent relapse. To elucidate the bioenergetic principle underlying AML relapse, we investigated the transcriptional regulation of mitochondrial-nuclear dual genomes responsible for metabolic plasticity in treatment-resistant blasts. Both the gain and loss of function results demonstrated that NFκB2, a noncanonical transcription factor (TF) of the NFκB (nuclear factor kappa-light-chain-enhancer of activated B cells) family, can control the expression of TFAM (mitochondrial transcription factor A), which is known to be essential for metabolic biogenesis. Furthermore, genetic tracking and promoter assays revealed that NFκB2 is in the mitochondria and can bind the specific "TTGGGGGGTG" region of the regulatory D-loop domain to activate the light-strand promoter (LSP) and heavy-strand promoter 1 (HSP1), promoters of the mitochondrial genome. Based on our discovery of NFκB2's novel function of regulating mitochondrial-nuclear dual genomes, we explored a novel triplet therapy including inhibitors of NFκB2, tyrosine kinase, and mitochondrial ATP synthase that effectively eliminated primary AML blasts with mutations of the FMS-related receptor tyrosine kinase 3 (FLT3) and displayed minimum toxicity to control cells ex vivo. As such, effective treatments for AML must include strong inhibitory actions on the dual genomes mediating metabolic plasticity to improve leukemia prognosis.

Draft Genome of Akame (Lates Japonicus) Reveals Possible Genetic Mechanisms for Long-Term Persistence and Adaptive Evolution with Low Genetic Diversity.

It is known that some endangered species have persisted for thousands of years despite their very small effective population sizes and low levels of genetic polymorphisms. To understand the genetic mechanisms of long-term persistence in threatened species, we determined the whole genome sequences of akame (Lates japonicus), which has survived for a long time with extremely low genetic variations. Genome-wide heterozygosity in akame was estimated to be 3.3 to 3.4 × 10-4/bp, one of the smallest values in teleost fishes. Analysis of demographic history revealed that the effective population size in akame was around 1,000 from 30,000 years ago to the recent past. The relatively high ratio of nonsynonymous to synonymous heterozygosity in akame indicated an increased genetic load. However, a detailed analysis of genetic diversity in the akame genome revealed that multiple genomic regions, including genes involved in immunity, synaptic development, and olfactory sensory systems, have retained relatively high nucleotide polymorphisms. This implies that the akame genome has preserved the functional genetic variations by balancing selection, to avoid a reduction in viability and loss of adaptive potential. Analysis of synonymous and nonsynonymous nucleotide substitution rates has detected signs of positive selection in many akame genes, suggesting adaptive evolution to temperate waters after the speciation of akame and its close relative, barramundi (Lates calcarifer). Our results indicate that the functional genetic diversity likely contributed to the long-term persistence of this species by avoiding the harmful effects of the population size reduction.

Leveraging Random Effects in Cistrome-Wide Association Studies for Decoding the Genetic Determinants of Prostate Cancer.

Cistrome-wide association studies (CWAS) are pivotal for identifying genetic determinants of diseases by correlating genetically regulated cistrome states with phenotypes. Traditional CWAS typically develops a model based on cistrome and genotype data to associate predicted cistrome states with phenotypes. The random effect cistrome-wide association study (RECWAS), reevaluates the necessity of cistrome state prediction in CWAS. RECWAS utilizes either a linear model or marginal effect for initial feature selection, followed by kernel-based feature aggregation for association testing is introduced. Through simulations and analysis of prostate cancer data, a thorough evaluation of CWAS and RECWAS is conducted. The results suggest that RECWAS offers improved power compared to traditional CWAS, identifying additional genomic regions associated with prostate cancer. CWAS identified 102 significant regions, while RECWAS found 50 additional significant regions compared to CWAS, many of which are validated. Validation encompassed a range of biological evidence, including risk signals from the GWAS catalog, susceptibility genes from the DisGeNET database, and enhancer-domain scores. RECWAS consistently demonstrated improved performance over traditional CWAS in identifying genomic regions associated with prostate cancer. These findings demonstrate the benefits of incorporating kernel methods into CWAS and provide new insights for genetic discovery in complex diseases.

Quantitative Trait Loci Mapping and Candidate Gene Analysis for Fiber Quality Traits in Upland Cotton

Superior fiber quality is one of the most important objectives in cotton breeding. To detect the genetic basis underlying fiber quality, an F2 population containing 413 plants was constructed by crossing Jifeng 914 and Jifeng 173, both of which have superior fiber quality, with Jifeng 173 being better. Five fiber quality traits were investigated in the F2, F2:3, F2:4, and F2:5 populations. Quantitative trait loci (QTL) mapping was conducted based on a high-density genetic map containing 11,488 single nucleotide polymorphisms (SNPs) and spanning 4202.12 cM in length. Transgressive segregation patterns and complex correlations in the five tested traits were observed. A total of 108 QTLs were found, including 13 major effect QTLs that contributed more than 10% toward phenotypic variation (PV) and 9 stable QTLs that could be repeatedly mapped in different generations. Chromosome A7 contained 12 QTL, ranking the first. No QTL was found on chromosomes D1 and D11. Two QTLs could be repeatedly detected in three populations, including qFL-D3-2 in F2, F2:4, and F2:5 with 9.18–21.45% of PV and qFS-A11-1 in F2:3, F2:4, and F2:5 with 6.05–10.41% of PV. Another seven stable QTLs could be detected in two populations, including four major effect QTLs: qFL-A12-3, qFS-D10-2, qMC-D6-2, and qMC-D8-1. Fourteen QTL-overlapping regions were found, which might explain the complex correlations among the five phenotypic traits. Four regions on chromosome A11, D3, D6, and D10 covered by both stable and major effect QTLs are promising for further fine mapping. The genomic regions of the two QTLs detected in three populations and the four major effect QTLs contain 810 genes. Gene functional analysis revealed that the annotated genes are mainly involved in protein binding and metabolic pathways. Fifteen candidate genes in the qFL-D3-2 region are highly expressed in fiber or ovules during fiber initiation, elongation, secondary cell wall thickening, or maturation stages. qRT-PCR revealed that Ghir_D03G005440.1 and Ghir_D03G011310.1 may play a role in promoting fiber initiation, while Ghir_D03G006470.1 may be beneficial for promoting fiber elongation. This study provides more information for revealing the molecular genetic basis underlying cotton fiber quality.

Long-term exposure to diesel exhaust particles induces concordant changes in DNA methylation and transcriptome in human adenocarcinoma alveolar basal epithelial cells

BackgroundDiesel exhaust particles (DEP), which contain hazardous compounds, are emitted during the combustion of diesel. As approximately one-third of the vehicles worldwide use diesel, there are growing concerns about the risks posed by DEP to human health. Long-term exposure to DEP is associated with airway hyperresponsiveness, pulmonary fibrosis, and inflammation; however, the molecular mechanisms behind the effects of DEP on the respiratory tract are poorly understood. Such mechanisms can be addressed by examining transcriptional and DNA methylation changes. Although several studies have focused on the effects of short-term DEP exposure on gene expression, research on the transcriptional effects and genome-wide DNA methylation changes caused by long-term DEP exposure is lacking. Hence, in this study, we investigated transcriptional and DNA methylation changes in human adenocarcinoma alveolar basal epithelial A549 cells caused by prolonged exposure to DEP and determined whether these changes are concordant.ResultsDNA methylation analysis using the Illumina Infinium MethylationEPIC BeadChips showed that the methylation levels of DEP-affected CpG sites in A549 cells changed in a dose-dependent manner; the extent of change increased with increasing dose reaching the statistical significance only in samples exposed to 30 µg/ml DEP. Four-week exposure to 30 µg/ml of DEP significantly induced DNA hypomethylation at 24,464 CpG sites, which were significantly enriched for DNase hypersensitive sites, genomic regions marked by H3K4me1 and H3K27ac, and several transcription factor binding sites. In contrast, 9,436 CpG sites with increased DNA methylation levels were significantly overrepresented in genomic regions marked by H3K27me3 as well as H3K4me1 and H3K27ac. In parallel, gene expression profiling by RNA sequencing demonstrated that long-term exposure to DEP altered the expression levels of 2,410 genes, enriching 16 gene sets including Xenobiotic metabolism, Inflammatory response, and Senescence. In silico analysis revealed that the expression levels of 854 genes correlated with the methylation levels of the DEP-affected cis-CpG sites.ConclusionsTo our knowledge, this is the first report of genome-wide transcriptional and DNA methylation changes and their associations in A549 cells following long-term exposure to DEP.

Genome-wide association studies from spoken phenotypic descriptions: a proof of concept from maize field studies.

We present a novel approach to genome-wide association studies (GWAS) by leveraging unstructured, spoken phenotypic descriptions to identify genomic regions associated with maize traits. Utilizing the Wisconsin Diversity panel, we collected spoken descriptions of Zea mays ssp. mays traits, converting these qualitative observations into quantitative data amenable to GWAS analysis. First, we determined that visually striking phenotypes could be detected from unstructured spoken phenotypic descriptions. Next, we developed two methods to process the same descriptions to derive the trait plant height, a well-characterized phenotypic feature in maize: (1) a semantic similarity metric that assigns a score based on the resemblance of each observation to the concept of 'tallness' and (2) a manual scoring system that categorizes and assigns values to phrases related to plant height. Our analysis successfully corroborated known genomic associations and uncovered novel candidate genes potentially linked to plant height. Some of these genes are associated with gene ontology terms that suggest a plausible involvement in determining plant stature. This proof-of-concept demonstrates the viability of spoken phenotypic descriptions in GWAS and introduces a scalable framework for incorporating unstructured language data into genetic association studies. This methodology has the potential not only to enrich the phenotypic data used in GWAS and to enhance the discovery of genetic elements linked to complex traits but also to expand the repertoire of phenotype data collection methods available for use in the field environment.

Digest: Evaluating sexual conflict in the genome.

Audet et al. (2024) investigate the genomic basis of sexual conflict in response to sexually discordant size selection in Drosophila melanogaster. They report interesting morphological changes in sexual dimorphism and multivariate allometry. Although they do not find any genetic variants that individually show a strong effect on these traits, as expected for a polygenic trait such as body size, they do find a region on chromosome 3L showing signs of sexually discordant selection (i.e., conflict between males and females). This study highlights potential genomic regions involved in sexual conflict, offering insights into sex-specific adaptations.

Into the Blue: Exploring genetic mechanisms behind the evolution of baleen whales

Marine ecosystems are ideal for studying evolutionary adaptations involved in lineage diversification due to few physical barriers and reduced opportunities for strict allopatry compared to terrestrial ecosystems. Cetaceans (whales, dolphins, and porpoises) are a diverse group of mammals that successfully adapted to various habitats within the aquatic environment around 50 million years ago. While the overall adaptive transition from terrestrial to fully aquatic species is relatively well understood, the radiation of modern whales is still unclear. Here high-quality genomes derived from previously published data were used to identify genomic regions that potentially underpinned the diversification of baleen whales (Balaenopteridae). A robust molecular phylogeny was reconstructed based on 10,159 single copy and complete genes for eight mysticetes, seven odontocetes and two cetacean outgroups. Analysis of positive selection across 3,150 genes revealed that balaenopterids have undergone numerous idiosyncratic and convergent genomic variations that may explain their diversification. Genes associated with aging, survival and homeostasis were enriched in all species. Additionally, positive selection on genes involved in the immune system were disclosed for the two largest species, blue and fin whales. Such genes can potentially be ascribed to their morphological evolution, allowing them to attain greater length and increased cell number. Further evidence is presented about gene regions that might have contributed to the extensive anatomical changes shown by cetaceans, including adaptation to distinct environments and diets. This study contributes to our understanding of the genomic basis of diversification in baleen whales and the molecular changes linked to their adaptive radiation, thereby enhancing our understanding of cetacean evolution.

Signatures of selection in Angus and Hanwoo beef cattle using imputed whole genome sequence data.

In this study, we detected signatures of selection in Hanwoo and Angus beef cattle using allele frequency and haplotype-based methods based on imputed whole genome sequence variants. Our dataset included 13,202 Angus animals with 10,057,633 imputed SNPs and 10,437 Hanwoo animals with 13,241,550 imputed SNPs. The dataset was subset down to 6,873,624 SNPs in common between the two populations to identify within population (runs of homozygosity, extended haplotype homozygosity) and between population signals of selection (allele fixation index, extended haplotype homozygosity). Assuming these selection signals were complementary to each other, they were combined into a decorrelated composite of multiple signals to identify regions under selection for each of the breeds. 27 genomic regions spanning 25.15Mb and harboring 360 genes were identified in Angus on chromosomes 1,3, 4, 5, 6, 7, 8, 12, 13, 14, 16, 20, 21 and 28. Similarly, in Hanwoo, 59 genes and 17 genomic regions spanning 5.21Mb on chromosomes 2, 4, 5, 6, 7, 8, 9, 10, 13, 17, 20 and 24 were identified. Apart from a small region on chromosome 13, there was no major overlap of selection signals between the two breeds reflecting their largely different selection histories, environmental challenges, breeding objectives and breed characteristics. Positional candidate genes identified in selected genomic regions in Angus have been previously associated with growth, immunity, reproductive development, feed efficiency and adaptation to environment while the candidate genes identified in Hanwoo included important genes regulating meat quality, fat deposition, cholesterol metabolism, lipid synthesis, neuronal development, and olfactory reception.

Blood Leukocyte DNA Methylation Markers of Periodontal Disease and Risk of Lung Cancer in a Case-Control Study Nested in the CLUE II Cohort.

Periodontal disease and DNA methylation markers have separately been associated with lung cancer risk. Examining methylation levels at genomic regions previously linked to periodontal disease may provide insights on the link between periodontal disease and lung cancer. In a nested case-control study drawn from the CLUE II cohort, we measured DNA methylation levels in 208 lung cancer cases and 208 controls. We examined the association between 37 DNA-methylated 5'-C-phosphate-G-3' (CpG) sites at three genomic regions, homeobox 4 (HOXA4), zinc finger protein (ZFP57), and a long noncoding RNA gene located in Chr10 (ENSG00000231601), and lung cancer risk. Statistically significant associations with lung cancer risk were observed for all 14 CpG sites from HOXA4 (OR ranging 1.41-1.62 for 1 SD increase in the DNA methylation level, especially within 15 years) and 1 CpG site on gene ENSG00000231601 (OR = 1.34 for 1 SD increase in the DNA methylation level). Although CpG sites on gene ZFP57 were not associated with lung cancer risk overall, statistically significant inverse associations were noted for six CpG sites when restricting follow-up to 15 years (OR = 0.73-0.77 for 1 SD increase in the DNA methylation level). Key methylation levels associated with periodontal disease are also associated with lung cancer risk. For both HOXA4 and ZFP57, the associations were stronger within 15 years of follow-up, which suggest that, if causal, the impact of methylation is acting late in the natural history of lung cancer. Identifying biological pathways that link periodontal disease and lung cancer could provide new opportunities for lung cancer detection and prevention.

Abstract Mo055: CHD4 Interacts With TBX5 to Maintain the Gene Regulatory Network of Postnatal Atrial Cardiomyocytes

Atrial fibrillation (AF) is the most common sustained arrhythmia in humans, leading to increased mortality due to increased incidence of stroke and heart failure. Over 200 genes are associated with AF by GWAS studies, suggesting a nuanced gene regulatory network (GRN) is required to maintain atrial cardiomyocyte (aCM) function and rhythm. Among these genes is the transcription factor TBX5, which we and others have shown is a master regulator of aCM identity that is essential for atrial rhythm homeostasis. We showed that TBX5 is required to maintain the chromatin looping, accessibility, and activity of atrial enhancers. Chromodomain helicase DNA binding protein 4 (CHD4), a chromatin remodeling protein that interacts with TBX5, is canonically thought to repress transcription as a component of the NuRD complex. Here, we investigated the contribution of CHD4 to TBX5 regulation of aCM gene expression. Methods: We inactivated CHD4 specifically within atrial cardiomyocytes using Cre recombinase driven by the aCM-specific Nppa promoter delivered postnatally with AAV9. Atrial remodeling, transcriptomic changes and AF susceptibility were examined. Results: We show that inactivating CHD4 in postnatal aCMs resulted in fibrotic atrial remodeling, spontaneous AF in a subset of mice, and increased AF susceptibility and burden after electrical pacing. Transcriptomic profiling of CHD4 aCM KO atria revealed gene expression changes that were highly correlated to altered gene expression downstream of TBX5 inactivation. Mechanistically, we show TBX5 recruits CHD4 to 33,170 genomic regions that are important for aCM gene regulation, whose expression is concordantly regulated downstream of TBX5 and CHD4. Together, these data reveal that CHD4 is required by TBX5 to promote accessibility of enhancers and activate many TBX5 target genes to maintain atrial rhythm.

Abstract Or106: DNA Double-stranded Breaks, a Hallmark of Aging, Defined at the Nucleotide Resolution, Are Increased and Associated with Transcription in the Cardiac Myocyte Genome in LMNA-Cardiomyopathy

Background: An intrinsic feature of gene transcription is the formation of DNA superhelices near the transcription bubble, which are resolved upon induction of transient double-stranded breaks (DSBs) by topoisomerases. Unrepaired DSBs are pathogenic as they lead to cell cycle arrest, senescence, inflammation, and organ dysfunction. Hypothesis and Aims: We posit that DSBs would be more prevalent at the genomic sites associated with gene expression. As Lamin-A (LMNA) regulates active transcription through Lamin associated domains (LADS) we predict for DSBs to be more common in non-LAD regions. The objectives were to identify and characterize genome-wide DSBs at the nucleotide resolution and determine the association of DSBs with transcription in cardiac myocytes. Methods: We identified the genome-wide DSBs in ~ 1 million cardiac myocytes per heart in 3 wild-type and 3 myocyte-specific LMNA-deficient ( Myh6-Cre:Lmna F/F ) mice by END-Sequencing, the LADs were defined in cardiac myocytes via LMNA CUT&RUN assay (N=5). Results: The prevalence of DSBs was 0.8% and 2.2% in the wild-type and Myh6-Cre:Lmna F/F myocytes, respectively. The END-Seq signals for 8 DSBs were enriched in the wild-type and 6,764 in Myh6-Cre:Lmna F/F myocytes, respectively (q<0.05). The DSBs were preferentially localized to the gene regions, transcription initiation sites, cardiac transcription factor motifs, and the G quadruplex forming structures. Because LMNA regulates transcription through the LADs, we defined the LADs in cardiac myocytes and found an average of 818 LADs per myocyte. LADs that were shared by at least three mice (N= 2,572) were considered constitutive LADs (cLADs), which comprised about a third of the mouse cardiac myocyte genomes. Transcript levels of the protein-coding genes located at the cLADs (N=3,975) were ~ 16-fold lower than those at the non-LAD regions (N=~17,778). The prevalence of DSBs was higher in the non-LAD as compared to the cLAD regions. Likewise, DSBs were more common in the loss-of LAD regions, defined as the juxtaposed genomic regions in the Myh6-Cre:Lmna F/F that corresponded to the LAD regions in the wild-type myocytes. Conclusions: To our knowledge, this is the first identification of DSBs, at the nucleotide resolution in the cardiovascular system. The prevalence of DSBs was higher in the genomic regions associated with transcription. Because transcription is pervasive, DSBs are expected to be common and pathogenic in various disease states and aging.

Abstract Mo084: Stem Cell-Based Functional Genomics Unravel A Novel Vascular Smooth Muscle Cell State Induced By The 9p21 Coronary Artery Disease Risk Locus

Background: Genome-wide association studies have identified common genetic variants (Single Nucleotide Polymorphisms – SNPs) at ~300 human genomic loci linked to coronary artery disease (CAD) susceptibility. Among these genomic regions, the most impactful is the 9p21.3 CAD risk locus, which spans a 60 kb gene desert and encompasses ~80 SNPs in high linkage disequilibrium. We have previously generated isogenic induced Pluripotent Stem Cells (iPSCs) lines from risk and non-risk donors at 9p21.3 and performed a complete deletion of the locus. By using iPSC-derived vascular smooth muscle cells (iPSC-VSMCs) we demonstrated that the risk haplotype at 9p21.3 causes altered expression of numerous genes essential for muscle function, including contraction and adhesion. Notably, deletion of the risk haplotype restores the non-risk phenotype, suggesting a gain of function effect. Here, we aimed to identify the transcriptomic signature and state trajectories of VSMCs induced by the 9p21.3 risk haplotype and validate them through functional assays and ex-vivo analysis of human coronary arteries. Methods: We have used iPSCs from individuals carrying the risk and non-risk haplotype at 9p21.3 and isogenic knockout lines. After VSMCs differentiation induction, mature VSMCs were used for single cell RNA sequencing using 10X Chromium platform and functional assays. Results: Our analysis revealed that the 9p21.3 risk haplotype prompts VSMCs to acquire a novel cellular state showing reduced plasticity and divergent from other VSMC states previously linked to CAD. We found markers of alternative lineages, including osteogenic and neuronal, coupled with altered expression of genes within other CAD loci, and functional defects. We identified a set of new molecular markers crucial to define risk-associated VSMCs and uncovered their upstream regulators. Conclusions: Our study provides insights into CAD pathogenesis driven by the 9p21.3 risk locus and identified new gene regulatory networks essential for maintaining the normal functionality of the muscle layer of the arteries. Leveraging the power of iPSCs we present novel concepts about the biology of VSMCs and shed light on the impact of the strongest CAD genetic risk factor in preserving a healthy cellular state within the vasculature.

TF-High-Evolutionary: In Vivo Mutagenesis of Gene Regulatory Networks for the Study of the Genetics and Evolution of the Drosophila Regulatory Genome.

Understanding the evolutionary potential of mutations in gene regulatory networks is essential to furthering the study of evolution and development. However, in multicellular systems, genetic manipulation of regulatory networks in a targeted and high-throughput way remains challenging. In this study, we designed TF-High-Evolutionary (HighEvo), a transcription factor (TF) fused with a base editor (activation-induced deaminase), to continuously induce germline mutations at TF-binding sites across regulatory networks in Drosophila. Populations of flies expressing TF-HighEvo in their germlines accumulated mutations at rates an order of magnitude higher than natural populations. Importantly, these mutations accumulated around the targeted TF-binding sites across the genome, leading to distinct morphological phenotypes consistent with the developmental roles of the tagged TFs. As such, this TF-HighEvo method allows the interrogation of the mutational space of gene regulatory networks at scale and can serve as a powerful reagent for experimental evolution and genetic screens focused on the regulatory genome.

Genes required for phosphosphingolipid formation in Caulobacter crescentus contribute to bacterial virulence.

Sphingolipids are ubiquitous in membranes of eukaryotes and are associated with important cellular functions. Although sphingolipids occur scarcely in bacteria, for some of them they are essential and, in other bacteria, they contribute to fitness and stability of the outer membrane, such as in the well-studied α-proteobacterium Caulobacter crescentus. We previously defined five structural genes for ceramide synthesis in C. crescentus, among them the gene for serine palmitoyltransferase, the enzyme that catalyzes the committed step of sphingolipid biosynthesis. Other mutants affected in genes of this same genomic region show cofitness with a mutant deficient in serine palmitoyltransferase. Here we show that at least two phosphosphingolipids are produced in C. crescentus and that at least another six gene products are needed for the decoration of ceramide upon phosphosphingolipid formation. All eleven genes participating in phosphosphingolipid formation are also required in C. crescentus for membrane stability and for displaying sensitivity towards the antibiotic polymyxin B. The genes for the formation of complex phosphosphingolipids are also required for C. crescentus virulence on Galleria mellonella insect larvae.

Mapping of a novel locus Ra conferring extreme resistance against potato virus A in cultivated potato (Solanum tuberosum L.).

The Ra extreme resistance against potato virus A was mapped to the upper of chromosome 4 in tetraploid potato. Potato virus A (PVA) is one of the major viruses affecting potato worldwide and can cause serious disease symptoms and yield losses. Previously, we determined that potato cultivar Barbara harbors Rysto (genotype: Ryryryry) and Ra (genotype: Rararara) that each independently confer extreme resistance to PVA. In this study, employing a combination of next-generation sequencing and bulked-segregant analysis, we further located this novel Ra on chromosome 4 using a tetraploid BC1 potato population derived from a Ry-free progeny (Rararararyryryry) of Barbara (RarararaRyryryry) × F58050 (rararararyryryry). Using 29 insertion-deletion (InDel) markers spanning chromosome 4, Ra was delimited by the InDel markers M8-83 and M10-8 within a genetic interval of 1.46cM, corresponding to a 1.86-Mb genomic region in the potato DM reference genome. The InDel marker M10-8, which is closely linked with the resistance against PVA in the Ry-free segregating populations, was then used to screen 43 selected Rysto-free tetraploid potato breeding clones. The phenotype to PVA was significantly correlated with the present/absent of the marker, albeit with a 9.3% false positive rate and a 14.0% false negative rate. These findings are of importance in furthering the cloning of Ra and employing the marker-assisted selection for PVA resistance.

Identification of piRNA Expression Profiles in Vitiligo

Objective: Vitiligo is a systemic dermatological disorder characterized by skin depigmentation due to melanocyte damage and dysfunction. The treatment of vitiligo remains a challenging aspect in dermatological practice. PIWI-interacting RNAs (piRNAs), approximately 24–32 nucleotides in length, are crucial in the epigenetic regulation of developmental processes. However, research on piRNA expression profiling in vitiligo is limited. This study aimed to identify and characterize piRNAs in skin tissues from vitiligo patients and healthy controls. Methods: We compared and analyzed piRNA expression profiles of human vitiligo skin by RNA sequencing. Results and Discussion: Our findings revealed a pronounced preference for uridine (U) at the 5′ end of piRNAs, predominantly originating from repeat and other genomic regions. Notably, we detected 73 differentially expressed piRNAs between the patient and control groups, with 40 piRNAs upregulated and 33 downregulated in vitiligo patients. Bioinformatics analysis indicated that the genes responsible for piRNA production were more prevalent in metabolic processes and implicated in regulating the PI3K-Akt and AMPK signaling pathways. Conclusions: The elucidation of piRNA expression profiles enhances our understanding of their functional roles in the pathogenesis of vitiligo.

Selection against domestication alleles in introduced rabbit populations.

Humans have moved domestic animals around the globe for thousands of years. These have occasionally established feral populations in nature, often with devastating ecological consequences. To understand how natural selection shapes re-adaptation into the wild, we investigated one of the most successful colonizers in history, the European rabbit. By sequencing the genomes of 297 rabbits across three continents, we show that introduced populations exhibit a mixed wild-domestic ancestry. We show that alleles that increased in frequency during domestication were preferentially selected against in novel natural environments. Interestingly, causative mutations for common domestication traits sometimes segregate at considerable frequencies if associated with less drastic phenotypes (for example, coat colour dilution), whereas mutations that are probably strongly maladaptive in nature are absent. Whereas natural selection largely targeted different genomic regions in each introduced population, some of the strongest signals of parallelism overlap genes associated with neuronal or brain function. This limited parallelism is probably explained by extensive standing genetic variation resulting from domestication together with the complex mixed ancestry of introduced populations. Our findings shed light on the selective and molecular mechanisms that enable domestic animals to re-adapt to the wild and provide important insights for the mitigation and management of invasive populations.

Recent secondary contact, genome-wide admixture, and asymmetric introgression of neo-sex chromosomes between two Pacific island bird species.

Secondary contact between closely related taxa represents a "moment of truth" for speciation-an opportunity to test the efficacy of reproductive isolation that evolved in allopatry and to identify the genetic, behavioral, and/or ecological barriers that separate species in sympatry. Sex chromosomes are known to rapidly accumulate differences between species, an effect that may be exacerbated for neo-sex chromosomes that are transitioning from autosomal to sex-specific inheritance. Here we report that, in the Solomon Islands, two closely related bird species in the honeyeater family-Myzomela cardinalis and Myzomela tristrami-carry neo-sex chromosomes and have come into recent secondary contact after ~1.1 my of geographic isolation. Hybrids of the two species were first observed in sympatry ~100 years ago. To determine the genetic consequences of hybridization, we use population genomic analyses of individuals sampled in allopatry and in sympatry to characterize gene flow in the contact zone. Using genome-wide estimates of diversity, differentiation, and divergence, we find that the degree and direction of introgression varies dramatically across the genome. For sympatric birds, autosomal introgression is bidirectional, with phenotypic hybrids and phenotypic parentals of both species showing admixed ancestry. In other regions of the genome, however, the story is different. While introgression on the Z/neo-Z-linked sequence is limited, introgression of W/neo-W regions and mitochondrial sequence (mtDNA) is highly asymmetric, moving only from the invading M. cardinalis to the resident M. tristrami. The recent hybridization between these species has thus enabled gene flow in some genomic regions but the interaction of admixture, asymmetric mate choice, and/or natural selection has led to the variation in the amount and direction of gene flow at sex-linked regions of the genome.