Genomic Regions Research Articles

Genome wide inherited modifications of the tomato epigenome by trans-activated bacterial CG methyltransferase

BackgroundEpigenetic variation is mediated by epigenetic marks such as DNA methylation occurring in all cytosine contexts in plants. CG methylation plays a critical role in silencing transposable elements and regulating gene expression. The establishment of CG methylation occurs via the RNA-directed DNA methylation pathway and CG methylation maintenance relies on METHYLTRANSFERASE1, the homologue of the mammalian DNMT1.PurposeHere, we examined the capacity to stably alter the tomato genome methylome by a bacterial CG-specific M.SssI methyltransferase expressed through the LhG4/pOP transactivation system.ResultsMethylome analysis of M.SssI expressing plants revealed that their euchromatic genome regions are specifically hypermethylated in the CG context, and so are most of their genes. However, changes in gene expression were observed only with a set of genes exhibiting a greater susceptibility to CG hypermethylation near their transcription start site. Unlike gene rich genomic regions, our analysis revealed that heterochromatic regions are slightly hypomethylated at CGs only. Notably, some M.SssI-induced hypermethylation persisted even without the methylase or transgenes, indicating inheritable epigenetic modification.ConclusionCollectively our findings suggest that heterologous expression of M.SssI can create new inherited epigenetic variations and changes in the methylation profiles on a genome wide scale. This open avenues for the conception of epigenetic recombinant inbred line populations with the potential to unveil agriculturally valuable tomato epialleles.

New Drosophila promoter-associated architectural protein Mzfp1 interacts with CP190 and is required for housekeeping gene expression and insulator activity.

In Drosophila, a group of zinc finger architectural proteins recruits the CP190 protein to the chromatin, an interaction that is essential for the functional activity of promoters and insulators. In this study, we describe a new architectural C2H2 protein called Madf and Zinc-Finger Protein 1 (Mzfp1) that interacts with CP190. Mzfp1 has an unusual structure that includes six C2H2 domains organized in a C-terminal cluster and two tandem MADF domains. Mzfp1 predominantly binds to housekeeping gene promoters located in both euchromatin and heterochromatin genome regions. In vivo mutagenesis studies showed that Mzfp1 is an essential protein, and both MADF domains and the CP190 interaction region are required for its functional activity. The C2H2 cluster is sufficient for the specific binding of Mzfp1 to regulatory elements, while the second MADF domain is required for Mzfp1 recruitment to heterochromatin. Mzfp1 binds to the proximal part of the Fub boundary that separates regulatory domains of the Ubx and abd-A genes in the Bithorax complex. Mzfp1 participates in Fub functions in cooperation with the architectural proteins Pita and Su(Hw). Thus, Mzfp1 is a new architectural C2H2 protein involved in the organization of active promoters and insulators in Drosophila.

Human enteric viruses’ detection in mussels (<i>Mytilus galloprovincialis</i>) farmed in the central Adriatic Sea

Human enteric viruses, such as hepatitis A virus (HAV), hepatitis E virus (HEV), and norovirus genogroups I and II (NoVGI and NoVGII), cause infections, and it has been largely demonstrated that mussels play an important role if consumed as raw or undercooked food matrices. This study aimed to investigate, through qualitative and quantitative biomolecular assays, the detection of partial genomic regions belonging to the most relevant enteropathogenic viruses for humans (HAV, HEV, NoVGI and NoVGII) in mussels (Mytilus galloprovincialis) farmed along the coasts of two Italian regions on the central Adriatic Sea: Abruzzo (Casalbordino, Chieti) and Molise (Termoli, Campobasso). A total of 425 animals were sampled, and the respective georeferentiations were registered. A total of 85 pools, each composed of five subjects/aliquot, were formed (22 from Abruzzo and 63 from Molise regions). This step was followed by homogenization and RNA extraction, and then the biomolecular assays [nested reverse transcription polymerase chain reaction (PCR) and real-time reverse transcription-quantitative PCR] were performed. 1.17% of the pool was positive for HAV RNA detection (102 copies/mL), 9.41% for HEV (102-103 copies/mL), 2.35% for NoVGI (101 copies/mL), and no pool was positive for NoVGII. This study demonstrated the human enteric viruses’ presence in mussels farmed in a low-investigated marine area. Based on a one-health point of view, this paper wants to enforce the importance of biomolecular and epidemiological screenings as surveillance systems to guarantee human, animal, and environmental health.

Identifying genomic regions associated with key agro‐morphological traits in soft red winter wheat using genome‐wide association study

AbstractAgro‐morphological traits play a significant role in the adaptation of wheat to diverse agroecosystems. Genetic understanding of these traits is crucial to develop cultivars adapted to specific environments and maximize their productivity. This is a comprehensive genome‐wide association study (GWAS) of 230 diverse lines of soft red winter wheat for identifying quantitative trait loci (QTLs) related to eight key agro‐morphological traits. The diversity panel was evaluated in two locations for three consecutive years (2020–2023). A total of 150 significant marker–trait associations were detected, including 65 for three flag leaf traits, 35 for peduncle length, 33 for plant height, 16 for heading date, and one for plant vigor using 27,466 single nucleotide polymorphism (SNP) markers. Eleven high‐confidence major‐effect QTLs explaining greater than 10% phenotypic variance were detected, of which seven were stable, and one showed an association with plant height and peduncle length. QTLs possibly allelic for important dwarfing (Rht23) and vernalization (Vrn‐B1) genes were identified. Six QTLs, QFlw.uga‐1A, QPdl.uga‐1A, QFlw.uga‐2B.2, QPdl.uga‐5A, QPdl.uga‐7A, and QPht.uga‐7B, are presumed to be novel, and nearby candidate gene(s) were identified for all except QPdl.uga‐1A. The pyramiding of favorable alleles from major‐effect QTLs was found to have significant improvement in peduncle length (shortened by 5 cm), flag leaf width (increased by 0.18 cm), and plant height (shortened by 11 cm). This study has improved our genetic understanding of important agro‐morphological traits. These results, upon further validation, can be used in breeding for desirable plant architecture to improve wheat yield potential.

Genetic response of a perennial grass to warm and wet environments interacts and is associated with trait means as well as plasticity.

Potential for rapid evolution is an important mechanism allowing species to adapt to changing climatic conditions. Although such a potential has been largely studied in various short-lived organisms, to what extent we can observe similar patterns in long-lived plant species, which often dominate natural systems, is largely unexplored. We explored potential for rapid evolution in Festuca rubra, long-lived grass with extensive clonal growth dominating in alpine grasslands. We used field sowing experiment simulating expected climate change in our model region. Specifically, we exposed seeds from 5 independent seed sources to novel climatic conditions by shifting them along a natural climatic grid and explored genetic profiles of established seedling after 3 years. Data on genetic profiles of plants selected under different novel conditions indicate that different climate shifts select significantly different pools of genotypes from common seed pools. Increasing soil moisture was more important selective pressure than increasing temperature or interaction of the two climatic factors. This can indicate negative genetic interaction in response to the combined effects, or that the effects of different climates are interactive rather than additive. The selected alleles were found in genomic regions likely affecting function of specific genes or their expression. Many of these were also linked to morphological traits (mainly to trait plasticity) suggesting these changes may have a consequence for plant performance. Overall, these data indicate that even long-lived plant species may experience strong selection by climate, and their populations thus have the potential to rapidly adapt to these novel conditions.

Abstract PR009: Exploring the functional consequences of APOBEC3-induced non-coding hotspot mutations in bladder cancer using massively parallel reporter assays and CRISPR-mediated base editing

Abstract Cancer research has traditionally focused on protein-coding genes, but approximately 98% of the human genome is non-coding, comprising gene regulatory elements such as enhancers, silencers, and promoters. Mutations within these regions can foster tumorigenesis by modulating the expression of oncogenes. Recent whole-genome sequencing studies have revealed several recurrently mutated non-coding genomic regions, including in bladder cancer, but their role in tumorigenesis remains unknown. In this study, we delved into whole-genome sequencing data from the bladder cancer TCGA cohort and an independent set of bladder tumors, uncovering hundreds of recurrent mutations. Mutational signature analysis identified over 65% of these mutations as APOBEC3-induced. We substantiated the involvement of APOBEC3 enzymes by using conventional deamination assays and novel approach dubbed as APOBEC3-mediated mutation capture assay. Predominantly, these mutations occurred in the hairpin loop, a favored motif for APOBEC3 enzymes. To systematically assess which of the mutations are functional and prioritize them for biological characterization, we developed Massively Parallel Reporter Assays (MPRA). Two types of MPRA libraries were constructed, one for evaluating promoter activity and the other for assessing enhancer activity of the mutations. Each library construction entailed 36,000 oligonucleotides, comprehensively representing all mutation hotspots and their corresponding wild-type alleles. The application of MPRA resulted in the identification of a substantial number of functionally relevant mutations. Zooming in on the top 20 genomic regions, we employed CRISPR-based approaches to pinpoint target genes and elucidate their functional consequences. Most of these regions co-localized with enhancers, as determined by analysis of ATAC-seq, ChIP-seq, and Hi-C data from bladder cancer cell lines. Perturbation of hotspots regions resulted in dysregulation of multiple genes. The co-regulation patterns of genes were further supported by strong correlation of their mRNA expression in RNA-seq datasets (TCGA, SRA, etc.). Notable examples included twin hotspot mutations upstream of the LEPROTL1 gene, elevating the expression of LEPROTL1, DCTN6 and SARAF, resulting in aggressive tumor growth in in vitro functional experiments. These three genes emerged as novel oncogenes in cancer. Another functional hotspot influenced the expression of two known long non-coding RNAs, NEAT1 and MALAT1, previously characterized as cancer genes. Additionally, some recurrently mutated regions or their target genes were associated with worse overall survival in the TCGA cohort. In conclusion, our findings propose that non-coding mutations likely drive tumorigenesis, contributing to tumor heterogeneity. Employing a combination of functional, multi-omics, and in vivo approaches, we aim to comprehensively characterize the role of hotspot non-coding mutations in bladder cancer and assess their clinical significance. Citation Format: Rouf Banday, Bilal Lone, Kelly Butler, Dhanusha Yesudhas, Alexandra Dobbins, Arup Chakraborty. Exploring the functional consequences of APOBEC3-induced non-coding hotspot mutations in bladder cancer using massively parallel reporter assays and CRISPR-mediated base editing [abstract]. In: Proceedings of the AACR Special Conference on Bladder Cancer: Transforming the Field; 2024 May 17-20; Charlotte, NC. Philadelphia (PA): AACR; Clin Cancer Res 2024;30(10_Suppl):Abstract nr PR009.

Insights into the genetic architecture of Phytophthora capsici root rot resistance in chile pepper (Capsicum spp.) from multi-locus genome-wide association study

BackgroundPhytophthora root rot, a major constraint in chile pepper production worldwide, is caused by the soil-borne oomycete, Phytophthora capsici. This study aimed to detect significant regions in the Capsicum genome linked to Phytophthora root rot resistance using a panel consisting of 157 Capsicum spp. genotypes. Multi-locus genome wide association study (GWAS) was conducted using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing (GBS). Individual plants were separately inoculated with P. capsici isolates, ‘PWB-185’, ‘PWB-186’, and ‘6347’, at the 4–8 leaf stage and were scored for disease symptoms up to 14-days post-inoculation. Disease scores were used to calculate disease parameters including disease severity index percentage, percent of resistant plants, area under disease progress curve, and estimated marginal means for each genotype.ResultsMost of the genotypes displayed root rot symptoms, whereas five accessions were completely resistant to all the isolates and displayed no symptoms of infection. A total of 55,117 SNP markers derived from GBS were used to perform multi-locus GWAS which identified 330 significant SNP markers associated with disease resistance. Of these, 56 SNP markers distributed across all the 12 chromosomes were common across the isolates, indicating association with more durable resistance. Candidate genes including nucleotide-binding site leucine-rich repeat (NBS-LRR), systemic acquired resistance (SAR8.2), and receptor-like kinase (RLKs), were identified within 0.5 Mb of the associated markers.ConclusionsResults will be used to improve resistance to Phytophthora root rot in chile pepper by the development of Kompetitive allele-specific markers (KASP®) for marker validation, genomewide selection, and marker-assisted breeding.

Abstract A013: Functional and clinical consequences of hotspot non-coding mutations upstream of the LEPROTL1 gene in bladder cancer

Abstract Recent whole-genome sequencing studies have unveiled intriguing hotspot mutations within non-coding genomic regions, yet their implications in cancer remain unknown. Here, we characterized the role of twin hotspot non-coding mutations upstream of the LEPROTL1 gene in bladder cancer: chr8: 29,952,919 G>A and chr8: 29,952,919 C>T. Either mutation is found in ~20% of bladder tumors, and the mutations do not typically co-occur. Through a comprehensive analysis integrating ChiP-seq, ATAC-seq, and Hi-C data from bladder cancer cell lines, we observed that these mutations co-localize with enhancer elements that interact with multiple nearby genes. In vitro dual luciferase reporter assays showed that both mutations are functional and likely affect gene expression in a genomic context. CRISPR-Cas9-mediated knocking out of the putative enhancer region identified LEPROTL1, DCTN6, and SARAF as target genes that are likely co-regulated, which was further supported by strong positive correlations of their mRNA expression in TCGA data. Isogenic cell clones generated by CRISPR-mediated base-editing showed that the mutations result in simultaneous upregulation of the three target genes, which were also expressed at higher levels in TCGA tumors harboring the non-coding mutations. Little is known about the function of LEPROTL1, DCTN6, and SARAF. Utilizing siRNA knockdown models complemented by live-cell imaging, viability assays, and cell cycle analysis, we found that they promote cell proliferation and cell cycle progression. Furthermore, our base-edited isogenic cell lines showed more rapid proliferation and cell cycle progression, as well as resistance to cisplatin. The mutations were also associated with shorter disease-specific survival in a TCGA cohort, and DCTN6 showed a negative association with overall survival. These results underscore the clinical significance of the non-coding mutations. Currently, we are assessing the effect of these mutations on tumor growth and cisplatin response in vivo using xenograft mouse models and evaluating the mechanisms by which LEPROTL1, DCTN6, and SARAF modulate cell cycle progression. Taken together, our data suggest that non-coding enhancer mutations near LEPROTL1 contribute to bladder carcinogenesis by increasing the expression of three co-regulated novel oncogenes, all of which were previously uncharacterized in bladder cancer. Since we found the mutations to be clinically significant, our study also raises the possibility of using the non-coding mutations and their target genes as clinical biomarkers and potential therapeutic drug targets. Citation Format: Kelly Butler, Bilal Lone, Alexandra Dobbins, Arup Chakraborty, Salah Boudajadi, Andrea Apolo, Rouf Banday. Functional and clinical consequences of hotspot non-coding mutations upstream of the LEPROTL1 gene in bladder cancer [abstract]. In: Proceedings of the AACR Special Conference on Bladder Cancer: Transforming the Field; 2024 May 17-20; Charlotte, NC. Philadelphia (PA): AACR; Clin Cancer Res 2024;30(10_Suppl):Abstract nr A013.

The Dark Side of the pollen: BSA-seq identified genomic regions linked to male sterility in globe artichoke

Globe artichoke (Cynara cardunculus var. scolymus; 2n = 2x = 34) is a food crop consumed for its immature flower heads. Traditionally, globe artichoke varietal types are vegetatively propagated. However, seed propagation makes it possible to treat the crop as annual, increasing field uniformity and reducing farmers costs, as well as pathogens diffusion. Despite globe artichoke’s significant agricultural value and the critical role of heterosis in the development of superior varieties, the production of hybrids remains challenging without a reliable system for large-scale industrial seed production. Male sterility (MS) presents a promising avenue for overcoming these challenges by simplifying the hybridization process and enabling cost-effective seed production. However, within the Cynara genus, genic male sterility has been linked to three recessive loci in globe artichoke, with no definitive genetic mechanism elucidated to date. A 250 offsprings F2 population, derived from a cross between a MS globe artichoke and a male fertile (MF) cultivated cardoon (C. cardunculus var. altilis) and fitting a monogenic segregation model (3:1), was analyzed through BSA-seq, aiming at the identification of genomic regions/genes affecting male sterility. Four QTL regions were identified on chromosomes 4, 12, and 14. By analyzing the sequence around the highest pick on chromosome 14, a cytochrome P450 (CYP703A2) was identified, carrying a deleterious substitution (R/Q) fixed in the male sterile parent. A single dCAPS marker was developed around this SNP, allowing the discrimination between MS and MF genotypes within the population, suitable for applications in plant breeding programs. A 3D model of the protein was generated by homology modeling, revealing that the mutated amino acid is part of a highly conserved motif crucial for protein folding.

P12 Transcriptional and epigenetic changes drive increased cellular plasticity in keloid fibroblasts

Abstract Introduction and aims Keloids are fibroproliferative scars that progressively grow beyond the original wound site. They occur as a result of genetic, epigenetic and environmental abnormalities. Fibroblasts are best characterized for their role in depositing and remodelling the extracellular matrix (ECM) and are believed to have a crucial role in keloid formation. However, the detailed molecular signatures that promote this pathogenic phenotype are still unknown. Methods To investigate the epigenetic and genetic changes in keloid fibroblasts, we performed assay for transposase-accessible chromatin (ATAC)-Seq and RNA-Seq on a patient-matched set of keloid and normal fibroblasts. ATAC-Seq is a high-throughput sequencing method that defines chromatin accessible regions in which the binding of sites of transcription factors that regulate gene expression can be mapped. Identified candidates are validated in keloid fibroblasts via quantitative polymerase chain reaction and their functional role is analysed in various in vitro assays. Results ATAC-Seq analysis revealed more accessibility in coding regions of keloid fibroblasts compared with normal fibroblasts. Pathways enriched in these regions include transcriptional regulation, embryonic development, cell differentiation, calcium-binding and collagen remodelling. Intriguingly, several developmental epigenetic modulators were identified that may alter the transcription and chromatin landscape of keloid fibroblasts to a more plastic state and have the capacity to influence cell differentiation choices, potentially promoting chondrogenic and adipogenic cell features. Currently, we are investigating the ability of keloid fibroblasts to be more readily differentiated into other cell types including induced pluripotent stem cells, and the contribution of those candidates on fibrotic features to the cell are also being explored. Conclusions Our analysis reveals that the chromatin of keloid fibroblasts is more accessible, especially in genomic regions associated with embryonic development and keloid fibroblasts show an increased plasticity in vitro. These findings may support the development of fibroblast-targeted therapies for keloid scars and possibly other connective tissue disorders.

Genome-wide association study identifies novel loci and candidate genes for rust resistance in wheat (Triticum aestivum L.)

BackgroundWheat rusts are important biotic stresses, development of rust resistant cultivars through molecular approaches is both economical and sustainable. Extensive phenotyping of large mapping populations under diverse production conditions and high-density genotyping would be the ideal strategy to identify major genomic regions for rust resistance in wheat. The genome-wide association study (GWAS) population of 280 genotypes was genotyped using a 35 K Axiom single nucleotide polymorphism (SNP) array and phenotyped at eight, 10, and, 10 environments, respectively for stem/black rust (SR), stripe/yellow rust (YR), and leaf/brown rust (LR).ResultsForty-one Bonferroni corrected marker-trait associations (MTAs) were identified, including 17 for SR and 24 for YR. Ten stable MTAs and their best combinations were also identified. For YR, AX-94990952 on 1A + AX-95203560 on 4A + AX-94723806 on 3D + AX-95172478 on 1A showed the best combination with an average co-efficient of infection (ACI) score of 1.36. Similarly, for SR, AX-94883961 on 7B + AX-94843704 on 1B and AX-94883961 on 7B + AX-94580041 on 3D + AX-94843704 on 1B showed the best combination with an ACI score of around 9.0. The genotype PBW827 have the best MTA combinations for both YR and SR resistance. In silico study identifies key prospective candidate genes that are located within MTA regions. Further, the expression analysis revealed that 18 transcripts were upregulated to the tune of more than 1.5 folds including 19.36 folds (TraesCS3D02G519600) and 7.23 folds (TraesCS2D02G038900) under stress conditions compared to the control conditions. Furthermore, highly expressed genes in silico under stress conditions were analyzed to find out the potential links to the rust phenotype, and all four genes were found to be associated with the rust phenotype.ConclusionThe identified novel MTAs, particularly stable and highly expressed MTAs are valuable for further validation and subsequent application in wheat rust resistance breeding. The genotypes with favorable MTA combinations can be used as prospective donors to develop elite cultivars with YR and SR resistance.

Fidelity Characterization of Highly Pathogenic Porcine Reproductive and Respiratory Syndrome Virus and NADC30-like Strain

The porcine reproductive and respiratory syndrome virus (PRRSV) has significantly impacted the global pork industry for over three decades. Its high mutation rates and frequent recombination greatly intensifies its epidemic and threat. To explore the fidelity characterization of Chinese highly pathogenic PRRSV JXwn06 and the NADC30-like strain CHsx1401, self-recombination and mutation in PAMs, MARC-145 cells, and pigs were assessed. In vitro, CHsx1401 displayed a higher frequency of recombination junctions and a greater diversity of junction types than JXwn06. In vivo, CHsx1401 exhibited fewer junction types yet maintained a higher junction frequency. Notably, JXwn06 showed more accumulation of mutations. To pinpoint the genomic regions influencing their fidelity, chimeric viruses were constructed, with the exchanged nsp9-10 regions between JXwn06 and CHsx1401. The SJn9n10 strain, which incorporates JXwn06’s nsp9-10 into the CHsx1401 genome, demonstrated reduced sensitivity to nucleotide analogs compared to CHsx1401. Conversely, compared with JXwn06, the JSn9n10 strain showed increased sensitivity to these inhibitors. The swapped nsp9-10 also influences the junction frequency and accumulated mutations as their donor strains. The results indicate a propensity for different types of genetic variations between these two strains and further highlight the nsp9-10 region as a critical determinant of their fidelity.

Characterization of extrachromosomal circular DNAs in plasma of patients with clear cell renal cell carcinoma.

Extrachromosomal circular DNAs (eccDNAs) have been recognized for their significant involvement in numerous biological processes. Nonetheless, the existence and molecular characteristics of eccDNA in the peripheral blood of patients diagnosed with clear cell renal cell carcinoma (ccRCC) have not yet been reported. Our aim was to identify potentially marked plasma eccDNAs in ccRCC patients. The detection of plasma eccDNA in ccRCC patients and healthy controls was performed using the Tn5-tagmentation and next-generation sequencing (NGS) method. Comparisons were made between ccRCC patients and healthy controls regarding the distribution of length, gene annotation, pattern of junctional nucleotide motif, and expression pattern of plasma eccDNA. We found 8,568 and 8,150 plasma eccDNAs in ccRCC patients and healthy controls, respectively. There were no statistical differences in the length distribution, gene annotation, and motif signature of plasma eccDNAs between the two groups. A total of 701 differentially expressed plasma eccDNAs were identified, and 25 plasma eccDNAs with potential diagnostic value for ccRCC have been successfully screened. These up-regulated plasma eccDNAs also be indicated to originate from the genomic region of the tumor-associated genes. This work demonstrates the characterization of plasma eccDNAs in ccRCC and suggests that the up-regulated plasma eccDNAs could be considered as a promising non-invasive biomarker in ccRCC.

Genomic loci associated with grain yield under well-watered and water-stressed conditions in multiple bi-parental maize populations

Smallholder maize farming systems in sub-Saharan Africa (SSA) are vulnerable to drought-induced yield losses, which significantly impact food security and livelihoods within these communities. Mapping and characterizing genomic regions associated with water stress tolerance in tropical maize is essential for future breeding initiatives targeting this region. In this study, three biparental F3 populations composed of 753 families were evaluated in Kenya and Zimbabwe and genotyped with high-density single nucleotide polymorphism (SNP) markers. Quantitative trait loci maping was performed on these genotypes to dissect the genetic architecture for grain yield (GY), plant height (PH), ear height (EH) and anthesis-silking interval (ASI) under well-watered (WW) and water-stressed (WS) conditions. Across the studied maize populations, mean GY exhibited a range of 4.55–8.55 t/ha under WW and 1.29–5.59 t/ha under WS, reflecting a 31–59% reduction range under WS conditions. Genotypic and genotype-by-environment (G × E) variances were significant for all traits except ASI. Overall broad sense heritabilities for GY were low to high (0.25–0.60). For GY, these genetic parameters were decreased under WS conditions. Linkage mapping revealed a significant difference in the number of QTLs detected, with 93 identified under WW conditions and 41 under WS conditions. These QTLs were distributed across all maize chromosomes. For GY, eight and two major effect QTLs (>10% phenotypic variation explained) were detected under WW and WS conditions, respectively. Under WS conditions, Joint Linkage Association Mapping (JLAM) identified several QTLs with minor effects for GY and revealed genomic region overlaps in the studied populations. Across the studied water regimes, five-fold cross-validation showed moderate to high prediction accuracies (−0.15–0.90) for GY and other agronomic traits. Our findings demonstrate the polygenic nature of WS tolerance and highlights the immense potential of using genomic selection in improving genetic gain in maize breeding.

The great tit HapMap project: A continental-scale analysis of genomic variation in a songbird.

A major aim of evolutionary biology is to understand why patterns of genomic diversity vary within taxa and space. Large-scale genomic studies of widespread species are useful for studying how environment and demography shape patterns of genomic divergence. Here, we describe one of the most geographically comprehensive surveys of genomic variation in a wild vertebrate to date; the great tit (Parus major) HapMap project. We screened ca 500,000 SNP markers across 647 individuals from 29 populations, spanning ~30 degrees of latitude and 40 degrees of longitude - almost the entire geographical range of the European subspecies. Genome-wide variation was consistent with a recent colonisation across Europe from a South-East European refugium, with bottlenecks and reduced genetic diversity in island populations. Differentiation across the genome was highly heterogeneous, with clear 'islands of differentiation', even among populations with very low levels of genome-wide differentiation. Low local recombination rates were a strong predictor of high local genomic differentiation (FST), especially in island and peripheral mainland populations, suggesting that the interplay between genetic drift and recombination causes highly heterogeneous differentiation landscapes. We also detected genomic outlier regions that were confined to one or more peripheral great tit populations, probably as a result of recent directional selection at the species' range edges. Haplotype-based measures of selection were related to recombination rate, albeit less strongly, and highlighted population-specific sweeps that likely resulted from positive selection. Our study highlights how comprehensive screens of genomic variation in wild organisms can provide unique insights into spatio-temporal evolutionary dynamics.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary.

The piRNA pathway is conserved to repress transposable element (TE) activity in the animal germline via a specialized class of small RNAs called piwi-interacting RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). While the molecular processes by which piCs function are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC origin and evolution, we use a population genomics approach to compare piC activity and sequence composition across 8 geographically distant strains of D. melanogaster with high quality long-read genome assemblies. We perform annotations of ovary piCs and genome-wide TE content in each strain. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs. Most TEs inferred to be recently active show an enrichment of insertions into old and large piCs, consistent with the previously proposed 'trap' model of piC evolution. By contrast, a small subset of active LTR families is enriched for the formation of new piCs, suggesting that these TEs have higher proclivity to form piCs. Thus, our findings uncover processes leading to the origin of piCs. We propose piC evolution begins with the emergence of piRNAs from a few specific LTR retrotransposon insertions that subsequently expand by accretion of other TE insertions during evolution to form larger 'trap' clusters. Our study shows that TEs themselves are the major force driving the rapid evolution of piCs.

Global Distribution and Molecular Evolution of Bat Coronaviruses

Bat coronaviruses cause a wide range of illnesses in humans and animals. Bats are known to harbor a wide diversity of Alphacoronaviruses and Betacoronaviruses. Betacoronaviruses have been linked to Severe Acute Respiratory Syndrome (SARS), Middle East Respiratory Syndrome (MERS), and other diseases such as gastroenteritis, bronchiolitis, and pneumonia. In the last 20 years, three betacoronaviruses emerged and caused widespread outbreaks in humans, including two deadly betacoronavirus epidemics, SARS-CoV, with mortality rate of 10%, and MERS-CoV, with mortality rate of 34.7%, and SARS-CoV-2, which caused the COVID-19 pandemic, with mortality rate of 3.4%. Studies have shown that bats are the main natural reservoirs for these viruses or their ancestral viruses. Observed variations in bat coronavirus genomes indicate that these viruses may have a potential to transmit to other hosts in close contact with humans and subsequently transmit to humans. As of today, there are no reported cases of direct coronavirus transmission from bats to humans. One reason for this might be that intermediate hosts are required for the transmission of bat coronaviruses to humans. Further studies are needed to map the amino acids and genomic regions responsible for the interactions between the spike of coronavirus and its receptors.

Genetic tapestry of Capsicum fruit colors: a comparative analysis of four cultivated species.

Genome-wide association study of color spaces across the four cultivated Capsicum spp. revealed a shared set of genes influencing fruit color, suggesting mechanisms and pathways across Capsicum species are conserved during the speciation. Notably, Cytochrome P450 of the carotenoid pathway, MYB transcription factor, and pentatricopeptide repeat-containing protein are the major genes responsible for fruit color variation across the Capsicum species. Peppers (Capsicum spp.) rank among the most widely consumed spices globally. Fruit color, serving as a determinant for use in food colorants and cosmeceuticals and an indicator of nutritional contents, significantly influences market quality and price. Cultivated Capsicum species display extensive phenotypic diversity, especially in fruit coloration. Our study leveraged the genetic variance within four Capsicum species (Capsicum baccatum, Capsicum chinense, Capsicum frutescens, and Capsicum annuum) to elucidate the genetic mechanisms driving color variation in peppers and related Solanaceae species. We analyzed color metrics and chromatic attributes (Red, Green, Blue, L*, a*, b*, Luminosity, Hue, and Chroma) on samples cultivated over six years (2015-2021). We resolved genomic regions associated with fruit color diversity through the sets of SNPs obtained from Genotyping by Sequencing (GBS) and genome-wide association study (GWAS) with a Multi-Locus Mixed Linear Model (MLMM). Significant SNPs with FDR correction were identified, within the Cytochrome P450, MYB-related genes, Pentatricopeptide repeat proteins, and ABC transporter family were the most common among the four species, indicating comparative evolution of fruit colors. We further validated the role of a pentatricopeptide repeat-containing protein (Chr01:31,205,460) and a cytochrome P450 enzyme (Chr08:45,351,919) via competitive allele-specific PCR (KASP) genotyping. Our findings advance the understanding of the genetic underpinnings of Capsicum fruit coloration, with developed KASP assays holding potential for applications in crop breeding and aligning with consumer preferences. This study provides a cornerstone for future research into exploiting Capsicum's diverse fruit color variation.

Behaviors of nucleosomes with mutant histone H4s in euchromatic domains of living human cells.

Since Robert Feulgen first stained DNA in the cell, visualizing genome chromatin has been a central issue in cell biology to uncover how chromatin is organized and behaves in the cell. To approach this issue, we have developed single-molecule imaging of nucleosomes, a basic unit of chromatin, to unveil local nucleosome behavior in living cells. In this study, we investigated behaviors of nucleosomes with various histone H4 mutants in living HeLa cells to address the role of H4 tail acetylation, including H4K16Ac and others, which are generally associated with more transcriptionally active chromatin regions. We ectopically expressed wild-type (wt) or mutated H4s (H4K16 point; H4K5,8,12,16 quadruple; and H4 tail deletion) fused with HaloTag in HeLa cells. Cells that expressed wtH4-Halo, H4K16-Halo mutants, and multiple H4-Halo mutants had euchromatin-concentrated distribution. Consistently, the genomic regions of the wtH4-Halo nucleosomes corresponded to Hi-C contact domains (or topologically associating domains, TADs)with active chromatin marks (A-compartment). Utilizing single-nucleosome imaging, we found that none of the H4 deacetylation or acetylation mimicked H4 mutants altered the overall local nucleosome motion. This finding suggests that H4 mutant nucleosomes embedded in the condensed euchromatic domains with excess endogenous H4 nucleosomes cannot cause an observable change in the local motion. Interestingly, H4 with four lysine-to-arginine mutations displayed a substantial freely diffusing fraction in the nucleoplasm, whereas H4 with a truncated N-terminal tail was incorporated in heterochromatic regions as well as euchromatin. Our study indicates the power of single-nucleosome imaging to understand individual histone/nucleosome behavior reflecting chromatin environments in living cells.

Genomic regions, candidate genes, and pleiotropic variants associated with physiological and anatomical indicators of heat stress response in lactating sows

BackgroundHeat stress (HS) poses significant threats to the sustainability of livestock production. Genetically improving heat tolerance could enhance animal welfare and minimize production losses during HS events. Measuring phenotypic indicators of HS response and understanding their genetic background are crucial steps to optimize breeding schemes for improved climatic resilience. The identification of genomic regions and candidate genes influencing the traits of interest, including variants with pleiotropic effects, enables the refinement of genotyping panels used to perform genomic prediction of breeding values and contributes to unraveling the biological mechanisms influencing heat stress response. Therefore, the main objectives of this study were to identify genomic regions, candidate genes, and potential pleiotropic variants significantly associated with indicators of HS response in lactating sows using imputed whole-genome sequence (WGS) data. Phenotypic records for 18 traits and genomic information from 1,645 lactating sows were available for the study. The genotypes from the PorcineSNP50K panel containing 50,703 single nucleotide polymorphisms (SNPs) were imputed to WGS and after quality control, 1,622 animals and 7,065,922 SNPs were included in the analyses.ResultsA total of 1,388 unique SNPs located on sixteen chromosomes were found to be associated with 11 traits. Twenty gene ontology terms and 11 biological pathways were shown to be associated with variability in ear skin temperature, shoulder skin temperature, rump skin temperature, tail skin temperature, respiration rate, panting score, vaginal temperature automatically measured every 10 min, vaginal temperature measured at 0800 h, hair density score, body condition score, and ear area. Seven, five, six, two, seven, 15, and 14 genes with potential pleiotropic effects were identified for indicators of skin temperature, vaginal temperature, animal temperature, respiration rate, thermoregulatory traits, anatomical traits, and all traits, respectively.ConclusionsPhysiological and anatomical indicators of HS response in lactating sows are heritable but highly polygenic. The candidate genes found are associated with important gene ontology terms and biological pathways related to heat shock protein activities, immune response, and cellular oxidative stress. Many of the candidate genes with pleiotropic effects are involved in catalytic activities to reduce cell damage from oxidative stress and cellular mechanisms related to immune response.