Comparative restriction enzyme analysis of methylation (CREAM) reveals methylome variability within a clonal in vitro cannabis population.

The primary focus of medicinal cannabis research is to ensure the stability of cannabis lines for consistent administration of chemically uniform products to patients. In recent years, tissue culture has emerged as a valuable technique for genetic preservation and rapid multiplication of cannabis clones. However, there is concern that the physical and chemical conditions of the growing media can induce somaclonal variation, potentially impacting the viability and uniformity of clones. To address this concern, we developed Comparative Restriction Enzyme Analysis of Methylation (CREAM), a novel method to assess DNA methylation patterns and used it to study a population of 78 cannabis clones maintained in tissue culture. Through bioinformatics analysis of the methylome, we successfully detected 2,272 polymorphic methylated regions among the clones. Remarkably, our results demonstrated that DNA methylation patterns were preserved across subcultures within the clonal population, allowing us to distinguish between two subsets of clonal lines used in this study. These findings significantly contribute to our understanding of the epigenetic variability within clonal lines in medicinal cannabis produced through tissue culture techniques. This knowledge is crucial for understanding the effects of tissue culture on DNA methylation and ensuring the consistency and reliability of medicinal cannabis products with therapeutic properties. Additionally, the CREAM method is a fast and affordable technology to get a first glimpse at methylation in a biological system. It offers a valuable tool for studying epigenetic variation in other plant species, thereby facilitating broader applications in plant biotechnology and crop improvement.

Impact of Combinatorial Histone Modifications on Acetyllysine Recognition by the ATAD2 and ATAD2B Bromodomains.

The ATPase family AAA+ domain containing 2 (ATAD2) protein and its paralog ATAD2B have a C-terminal bromodomain (BRD) that functions as a reader of acetylated lysine residues on histone proteins. Using a structure-function approach, we investigated the ability of the ATAD2/B BRDs to select acetylated lysine among multiple histone post-translational modifications. The ATAD2B BRD can bind acetylated histone ligands that also contain adjacent methylation or phosphorylation marks, while the presence of these modifications significantly weakened the acetyllysine binding activity of the ATAD2 BRD. Our structural studies provide mechanistic insights into how ATAD2/B BRD-binding pocket residues coordinate the acetyllysine group in the context of adjacent post-translational modifications. Furthermore, we investigated how sequence changes in amino acids of the histone ligands impact the recognition of an adjacent acetyllysine residue. Our study highlights how the interplay between multiple combinations of histone modifications influences the reader activity of the ATAD2/B BRDs, resulting in distinct binding modes.

Hollow Polyethyleneimine Nanoparticles with Drug Loaded DNA for Chemotherapeutic Applications.

The next generation of anticancer agents are emerging from rationally designed nanostructured materials. This work involved the synthesis and characterization of novel hollow DNA-conjugated gold nanoparticles (DNA-AuNPs) for controlled drug delivery. Polyethyleneimine (PEI) was bound to AuNPs, forming polymer-shell nanoparticles. Dissolution of the gold core via iodine formed hollow core polymeric nanoparticles (HCPNPs) and a high density (85 molecules/particle) of DNA intercalated with daunorubicin was conjugated. Particles were spherical with an average diameter of 105.7±17.3 nm and zeta potential of 20.4±3.54 mV. We hypothesize the DNA backbone electrostatically condensed to the primary amines on the surface of the particle toroidally, weaving itself within the polymer shell. During the DNA intercalation process, increasing the ionic concentration and decreasing the amine/phosphate ratio 10-fold increased drug intercalation 64 % and 61 %, respectively, allowing us to determine the optimal method of particle synthesis. As intercalation sites increased with increasing DNA strand length, drug loading increased. An average of 874±40.1 daunorubicin molecules were loaded per HCPNP. HCPNPs with drug intercalated DNA have strong potential to be clinically efficacious drug delivery vehicles due to the versatility of DNA and high drug loading capacities.

The invasive land flatworm Arthurdendyus triangulatus has repeated sequences in the mitogenome, extra-long cox2 gene and paralogous nuclear rRNA clusters.

Using a combination of short- and long-reads sequencing, we were able to sequence the complete mitochondrial genome of the invasive 'New Zealand flatworm' Arthurdendyus triangulatus (Geoplanidae, Rhynchodeminae, Caenoplanini) and its two complete paralogous nuclear rRNA gene clusters. The mitogenome has a total length of 20,309 bp and contains repetitions that includes two types of tandem-repeats that could not be solved by short-reads sequencing. We also sequenced for the first time the mitogenomes of four species of Caenoplana (Caenoplanini). A maximum likelihood phylogeny associated A. triangulatus with the other Caenoplanini but Parakontikia ventrolineata and Australopacifica atrata were rejected from the Caenoplanini and associated instead with the Rhynchodemini, with Platydemus manokwari. It was found that the mitogenomes of all species of the subfamily Rhynchodeminae share several unusual structural features, including a very long cox2 gene. This is the first time that the complete paralogous rRNA clusters, which differ in length, sequence and seemingly number of copies, were obtained for a Geoplanidae.

Genomic and Phenotypic Analysis of Salmonella enterica Bacteriophages Identifies Two Novel Phage Species.

Bacteriophages (phages) are potential alternatives to chemical antimicrobials against pathogens of public health significance. Understanding the diversity and host specificity of phages is important for developing effective phage biocontrol approaches. Here, we assessed the host range, morphology, and genetic diversity of eight Salmonella enterica phages isolated from a wastewater treatment plant. The host range analysis revealed that six out of eight phages lysed more than 81% of the 43 Salmonella enterica isolates tested. The genomic sequences of all phages were determined. Whole-genome sequencing (WGS) data revealed that phage genome sizes ranged from 41 to 114 kb, with GC contents between 39.9 and 50.0%. Two of the phages SB13 and SB28 represent new species, Epseptimavirus SB13 and genera Macdonaldcampvirus, respectively, as designated by the International Committee for the Taxonomy of Viruses (ICTV) using genome-based taxonomic classification. One phage (SB18) belonged to the Myoviridae morphotype while the remaining phages belonged to the Siphoviridae morphotype. The gene content analyses showed that none of the phages possessed virulence, toxin, antibiotic resistance, type I-VI toxin-antitoxin modules, or lysogeny genes. Three (SB3, SB15, and SB18) out of the eight phages possessed tailspike proteins. Whole-genome-based phylogeny of the eight phages with their 113 homologs revealed three clusters A, B, and C and seven subclusters (A1, A2, A3, B1, B2, C1, and C2). While cluster C1 phages were predominantly isolated from animal sources, cluster B contained phages from both wastewater and animal sources. The broad host range of these phages highlights their potential use for controlling the presence of S. enterica in foods.

High-density genetic linkage mapping in Sitka spruce advances the integration of genomic resources in conifers.

In species with large and complex genomes such as conifers, dense linkage maps are a useful resource for supporting genome assembly and laying the genomic groundwork at the structural, populational, and functional levels. However, most of the 600+ extant conifer species still lack extensive genotyping resources, which hampers the development of high-density linkage maps. In this study, we developed a linkage map relying on 21,570 single nucleotide polymorphism (SNP) markers in Sitka spruce (Picea sitchensis [Bong.] Carr.), a long-lived conifer from western North America that is widely planted for productive forestry in the British Isles. We used a single-step mapping approach to efficiently combine RAD-seq and genotyping array SNP data for 528 individuals from 2 full-sib families. As expected for spruce taxa, the saturated map contained 12 linkages groups with a total length of 2,142 cM. The positioning of 5,414 unique gene coding sequences allowed us to compare our map with that of other Pinaceae species, which provided evidence for high levels of synteny and gene order conservation in this family. We then developed an integrated map for P. sitchensis and Picea glauca based on 27,052 markers and 11,609 gene sequences. Altogether, these 2 linkage maps, the accompanying catalog of 286,159 SNPs and the genotyping chip developed, herein, open new perspectives for a variety of fundamental and more applied research objectives, such as for the improvement of spruce genome assemblies, or for marker-assisted sustainable management of genetic resources in Sitka spruce and related species.

Automated single-cell proteomics providing sufficient proteome depth to study complex biology beyond cell type classifications.

Mass spectrometry (MS)-based single-cell proteomics (SCP) has gained massive attention as a viable complement to other single cell approaches. The rapid technological and computational advances in the field have pushed the boundaries of sensitivity and throughput. However, reproducible quantification of thousands of proteins within a single cell at reasonable proteome depth to characterize biological phenomena remains a challenge. To address some of those limitations we present a combination of fully automated single cell sample preparation utilizing a dedicated chip within the picolitre dispensing robot, the cellenONE. The proteoCHIP EVO 96 can be directly interfaced with the Evosep One chromatographic system for in-line desalting and highly reproducible separation with a throughput of 80 samples per day. This, in combination with the Bruker timsTOF MS instruments, demonstrates double the identifications without manual sample handling. Moreover, relative to standard high-performance liquid chromatography, the Evosep One separation provides further 2-fold improvement in protein identifications. The implementation of the newest generation timsTOF Ultra with our proteoCHIP EVO 96-based sample preparation workflow reproducibly identifies up to 4,000 proteins per single HEK-293T without a carrier or match-between runs. Our current SCP depth spans over 4 orders of magnitude and identifies over 50 biologically relevant ubiquitin ligases. We complement our highly reproducible single-cell proteomics workflow to profile hundreds of lipopolysaccharide (LPS)-perturbed THP-1 cells and identified key regulatory proteins involved in interleukin and interferon signaling. This study demonstrates that the proteoCHIP EVO 96-based SCP sample preparation with the timsTOF Ultra provides sufficient proteome depth to study complex biology beyond cell-type classifications.

Mitochondrial genome sequence of the protist Ancyromonas sigmoides Kent, 1881 (Ancyromonadida) from the Sugluk Inlet, Hudson Strait, Nunavik, Québec.

Introduction: There is little information on evolutionarily ancient eukaryotes, which are often referred to as basal eukaryotes, in Arctic waters. Despite earlier studies being conducted in the Russian White Sea, only few have been reported. Methods: Following a shotgun sequence survey of diatom cultures from Sugluk Inlet off the Hudson Strait in Northern Québec, we obtained the complete mitochondrial genome and the operon of nuclear ribosomal RNA genes from a strain that matches that of Ancyromonas sigmoides (Kent, 1881). Results: The sequence of the mitogenome retrieved was 41,889 bp in length and encoded 38 protein-coding genes, 5 non-conserved open-reading frames, and 2 rRNA and 24 tRNA genes. The mitogenome has retained sdh2 and sdh3, two genes of the succinate dehydrogenase complex, which are sometimes found among basal eukaryotes but seemingly missing among the Malawimonadidae, a lineage sister to Ancyromonadida in some phylogenies. The phylogeny inferred from the 18S rRNA gene associated A. sigmoides from Sugluk Inlet with several other strains originating from the Arctic. The study also unveiled the presence of a metagenomic sequence ascribed to bacteria in GenBank, but it was clearly a mitochondrial genome with a gene content highly similar to that of A. sigmoides, including the non-conserved open-reading frames. Discussion: After re-annotation, a phylogeny was inferred from mitochondrial protein sequences, and it strongly associated A. sigmoides with the misidentified organism, with the two being possibly conspecific or sibling species as they are more similar to one another than to species of the genus Malawimonas. Overall our phylogeny showed that the ice associated ancryomonads were clearly distinct from more southerly strains.

Genomic diversity of wild and cultured Yesso scallop Mizuhopecten yessoensis from Japan and Canada.

The Yesso scallop Mizuhopecten yessoensis is an important aquaculture species that was introduced to Western Canada from Japan to establish an economically viable scallop farming industry. This highly fecund species has been propagated in Canadian aquaculture hatcheries for the past 40 years, raising questions about genetic diversity and genetic differences among hatchery stocks. In this study, we compare cultured Canadian and wild Japanese populations of Yesso scallop using double-digest restriction site-associated DNA (ddRAD) sequencing to genotype 21,048 variants in 71 wild-caught scallops from Japan, 65 scallops from the Vancouver Island University breeding population, and 37 scallops obtained from a commercial farm off Vancouver Island, British Columbia. The wild scallops are largely comprised of equally unrelated individuals, whereas cultured scallops are comprised of multiple families of related individuals. The polymorphism rate estimated in wild scallops was 1.7%, whereas in the cultured strains, it ranged between 1.35 and 1.07%. Interestingly, heterozygosity rates were highest in the cultured populations, which is likely due to shellfish hatchery practices of crossing divergent strains to gain benefits of heterosis and to avoid inbreeding. Evidence of founder effects and drift was observed in the cultured strains, including high genetic differentiation between cultured populations and between cultured populations and the wild population. Cultured populations had effective population sizes ranging from 9 to 26 individuals whereas the wild population was estimated at 25,048-56,291 individuals. Further, a depletion of low-frequency variants was observed in the cultured populations. These results indicate significant genetic diversity losses in cultured scallops in Canadian breeding programs.

The mitochondrial genome of the bioluminescent fish Malacosteusniger Ayres, 1848 (Stomiidae, Actinopterygii) is large and complex, and contains an inverted-repeat structure.

We determined the complete mitogenome sequence of the bioluminescent fish Malacosteusniger using long-read sequencing technologies. The 21,263 bp mitogenome features a complex structure with two copies of a 1198-bp inverted-repeat and a region of 2616-bp containing alternating copies of 16 and 26 bp repeat elements. Whole mitogenome phylogenies inferred from both nucleotide and amino-acid datasets place M.niger among Melanostomiinae. The need for additional complete mitogenome sequences from the subfamily Malacosteinae is discussed.

Single nucleotide variants in Pseudomonas aeruginosa populations from sputum correlate with baseline lung function and predict disease progression in individuals with cystic fibrosis.

The severity and progression of lung disease are highly variable across individuals with cystic fibrosis (CF) and are imperfectly predicted by mutations in the human gene CFTR, lung microbiome variation or other clinical factors. The opportunistic pathogen Pseudomonas aeruginosa (Pa) dominates airway infections in most CF adults. Here we hypothesized that within-host genetic variation of Pa populations would be associated with lung disease severity. To quantify Pa genetic variation within CF sputum samples, we used deep amplicon sequencing (AmpliSeq) of 209 Pa genes previously associated with pathogenesis or adaptation to the CF lung. We trained machine learning models using Pa single nucleotide variants (SNVs), microbiome diversity data and clinical factors to classify lung disease severity at the time of sputum sampling, and to predict lung function decline after 5 years in a cohort of 54 adult CF patients with chronic Pa infection. Models using Pa SNVs alone classified lung disease severity with good sensitivity and specificity (area under the receiver operating characteristic curve: AUROC=0.87). Models were less predictive of lung function decline after 5 years (AUROC=0.74) but still significantly better than random. The addition of clinical data, but not sputum microbiome diversity data, yielded only modest improvements in classifying baseline lung function (AUROC=0.92) and predicting lung function decline (AUROC=0.79), suggesting that Pa AmpliSeq data account for most of the predictive value. Our work provides a proof of principle that Pa genetic variation in sputum tracks lung disease severity, moderately predicts lung function decline and could serve as a disease biomarker among CF patients with chronic Pa infections.

3D-GBS: a universal genotyping-by-sequencing approach for genomic selection and other high-throughput low-cost applications in species with small to medium-sized genomes.

Despite the increased efficiency of sequencing technologies and the development of reduced-representation sequencing (RRS) approaches allowing high-throughput sequencing (HTS) of multiplexed samples, the per-sample genotyping cost remains the most limiting factor in the context of large-scale studies. For example, in the context of genomic selection (GS), breeders need genome-wide markers to predict the breeding value of large cohorts of progenies, requiring the genotyping of thousands candidates. Here, we introduce 3D-GBS, an optimized GBS procedure, to provide an ultra-high-throughput and ultra-low-cost genotyping solution for species with small to medium-sized genome and illustrate its use in soybean. Using a combination of three restriction enzymes (PstI/NsiI/MspI), the portion of the genome that is captured was reduced fourfold (compared to a "standard" ApeKI-based protocol) while reducing the number of markers by only 40%. By better focusing the sequencing effort on limited set of restriction fragments, fourfold more samples can be genotyped at the same minimal depth of coverage. This GBS protocol also resulted in a lower proportion of missing data and provided a more uniform distribution of SNPs across the genome. Moreover, we investigated the optimal number of reads per sample needed to obtain an adequate number of markers for GS and QTL mapping (500-1000 markers per biparental cross). This optimization allows sequencing costs to be decreased by ~ 92% and ~ 86% for GS and QTL mapping studies, respectively, compared to previously published work. Overall, 3D-GBS represents a unique and affordable solution for applications requiring extremely high-throughput genotyping where cost remains the most limiting factor.

Comparison of genotyping by sequencing procedures to determine population genetic structure.

Advancements in technology over the past few decades have resulted in the development of genome sequencing at lower costs. Protocols, costs, and the amount of data produced by different sequencing technologies are highly variable. Ion Torrent and Illumina sequencing instruments are two sequencing technologies which use very similar library preparation procedures. Enzymatic combinations can be changed in genotyping by sequencing (GbS) library protocols without significant adjustments. To compare the outputs from two different GbS procedures, we sequenced samples of two sister species of yellow-nosed albatross collected at multiple geographic locations. The data sets involving different sequencing instruments and enzymatic combinations were analysed using the Stacks pipeline and aligned to the same reference genome. Both procedures identified the same genetic clusters separating Atlantic and Indian yellow-nosed albatross and substructure within Indian yellow-nosed albatross.

Spruce giga-genomes: structurally similar yet distinctive with differentially expanding gene families and rapidly evolving genes.

Spruces (Picea spp.) are coniferous trees widespread in boreal and mountainous forests of the northern hemisphere, with large economic significance and enormous contributions to global carbon sequestration. Spruces harbor very large genomes with high repetitiveness, hampering their comparative analysis. Here, we present and compare the genomes of four different North American spruces: the genome assemblies for Engelmann spruce (Picea engelmannii) and Sitka spruce (Picea sitchensis) together with improved and more contiguous genome assemblies for white spruce (Picea glauca) and for a naturally occurring introgress of these three species known as interior spruce (P. engelmannii × glauca × sitchensis). The genomes were structurally similar, and a large part of scaffolds could be anchored to a genetic map. The composition of the interior spruce genome indicated asymmetric contributions from the three ancestral genomes. Phylogenetic analysis of the nuclear and organelle genomes revealed a topology indicative of ancient reticulation. Different patterns of expansion of gene families among genomes were observed and related with presumed diversifying ecological adaptations. We identified rapidly evolving genes that harbored high rates of non-synonymous polymorphisms relative to synonymous ones, indicative of positive selection and its hitchhiking effects. These gene sets were mostly distinct between the genomes of ecologically contrasted species, and signatures of convergent balancing selection were detected. Stress and stimulus response was identified as the most frequent function assigned to expanding gene families and rapidly evolving genes. These two aspects of genomic evolution were complementary in their contribution to divergent evolution of presumed adaptive nature. These more contiguous spruce giga-genome sequences should strengthen our understanding of conifer genome structure and evolution, as their comparison offers clues into the genetic basis of adaptation and ecology of conifers at the genomic level. They will also provide tools to better monitor natural genetic diversity and improve the management of conifer forests. The genomes of four closely related North American spruces indicate that their high similarity at the morphological level is paralleled by the high conservation of their physical genome structure. Yet, the evidence of divergent evolution is apparent in their rapidly evolving genomes, supported by differential expansion of key gene families and large sets of genes under positive selection, largely in relation to stimulus and environmental stress response.

The Spruce Budworm Genome: Reconstructing the Evolutionary History of Antifreeze Proteins.

Insects have developed various adaptations to survive harsh winter conditions. Among freeze-intolerant species, some produce "antifreeze proteins" (AFPs) that bind to nascent ice crystals and inhibit further ice growth. Such is the case of the spruce budworm, Choristoneura fumiferana (Lepidoptera: Tortricidae), a destructive North American conifer pest that can withstand temperatures below -30°C. Despite the potential importance of AFPs in the adaptive diversification of Choristoneura, genomic tools to explore their origins have until now been limited. Here we present a chromosome-scale genome assembly for C. fumiferana, which we used to conduct comparative genomic analyses aimed at reconstructing the evolutionary history of tortricid AFPs. The budworm genome features 16 genes homologous to previously reported C. fumiferana AFPs (CfAFPs), 15 of which map to a single region on chromosome 18. Fourteen of these were also detected in five congeneric species, indicating Choristoneura AFP diversification occurred before the speciation event that led to C. fumiferana. Although budworm AFPs were previously considered unique to the genus Choristoneura, a search for homologs targeting recently sequenced tortricid genomes identified seven CfAFP-like genes in the distantly related Notocelia uddmanniana. High structural similarity between Notocelia and Choristoneura AFPs suggests a common origin, despite the absence of homologs in three related tortricids. Interestingly, one Notocelia AFP formed the C-terminus of a "zonadhesin-like" protein, possibly representing the ancestral condition from which tortricid AFPs evolved. Future work should clarify the evolutionary path of AFPs between Notocelia and Choristoneura and assess the role of the "zonadhesin-like" protein as precursor of tortricid AFPs.

Neurocognitive effects of repeated ketamine infusion treatments in patients with treatment resistant depression: a retrospective chart review.

BACKGROUND: Ketamine has emerged as a rapid-acting antidepressant in treatment-resistant depression (TRD) increasingly used in non-research, clinical settings. Few studies, however, have examined neurocognitive effects of repeated racemic ketamine infusion treatments in patients with TRD. In an effort to identify potential effects after serial infusions, we conducted a retrospective chart review to identify statistically significant changes in cognition in patient undergoing serial intravenous infusions; concomitantly, we examined baseline cognition as potential predictor of anti-depressant potential. METHODS: Twenty-two patients with TRD were examined after they finished the induction phase of 8-10 repeated intravenous ketamine infusions and completed the assessments of their depressive symptoms (measured by the 16-item Quick Inventory of Depressive Symptomatology-Self Report Scale: QIDS-SR16) and cognitive function (measured by the Montreal Cognitive Assessment: MoCA) before the first and the last ketamine treatments. RESULTS: Repeated ketamine infusions administered through an escalating dose protocol with 8-10 infusion sessions produced a 47.2% reduction response in depression; there was no evidence of impairment as reflected in MoCA testing. There was a moderate association between baseline cognition and antidepressant response with a Pearson correlation of 0.453. CONCLUSION: In this naturalistic sample of patients with TRD in our clinical service, repeated ketamine infusions significantly decreased depression symptoms without impairing cognitive performance. The baseline cognition may positively predict antidepressant responses of repeated ketamine treatment.

Salmonella enterica subsp. enterica virulence potential can be linked to higher survival within a dynamic in vitro human gastrointestinal model.

Salmonella enterica subsp. enterica is one of the leading causes of human foodborne infections and several outbreaks are now associated with the consumption of fresh fruit and vegetables. This study aims at evaluating whether Salmonella virulence can be linked to an enhanced ability to survive successive digestive environments. Thirteen S. enterica strains were selected according to high and low virulence phenotypes. Lettuce inoculated separately with each S. enterica strain was used as food matrix in the TNO gastrointestinal model (TIM-1) of the human upper gastrointestinal tract. During the passage in the stomach, counts determined using PMA-qPCR were 2-5 logs higher than the cultivable counts for all strains indicating the presence of viable but non-cultivable cells. Bacterial growth was observed in the duodenum compartment after 180 min for all but one strain and growth continued into the ileal compartment. After passage through the simulated gastrointestinal tract, both virulent and avirulent S. enterica strains survived but high virulence strains had a significantly (p = 0.004) better average survival rate (1003 %-3753 %) than low virulence strains (from 25 % to 3730%). The survival rates of S. enterica strains could be linked to the presence of genes associated with acid and bile resistance and their predicted products. The presence of single nucleotide polymorphisms may also impact the function of virulence associated genes and play a role in the resulting phenotype. These data provide an understanding of the relationship between measured virulence potential and survival of S. enterica during dynamic simulated gastrointestinal transit.

GWAS identifies an ortholog of the rice D11 gene as a candidate gene for grain size in an international collection of hexaploid wheat.

Grain size is a key agronomic trait that contributes to grain yield in hexaploid wheat. Grain length and width were evaluated in an international collection of 157 wheat accessions. These accessions were genetically characterized using a genotyping-by-sequencing (GBS) protocol that produced 73,784 single nucleotide polymorphism (SNP) markers. GBS-derived genotype calls obtained on Chinese Spring proved extremely accurate when compared to the reference (> 99.9%) and showed > 95% agreement with calls made at SNP loci shared with the 90 K SNP array on a subset of 71 Canadian wheat accessions for which both types of data were available. This indicates that GBS can yield a large amount of highly accurate SNP data in hexaploid wheat. The genetic diversity analysis performed using this set of SNP markers revealed the presence of six distinct groups within this collection. A GWAS was conducted to uncover genomic regions controlling variation for grain length and width. In total, seven SNPs were found to be associated with one or both traits, identifying three quantitative trait loci (QTLs) located on chromosomes 1D, 2D and 4A. In the vicinity of the peak SNP on chromosome 2D, we found a promising candidate gene (TraesCS2D01G331100), whose rice ortholog (D11) had previously been reported to be involved in the regulation of grain size. These markers will be useful in breeding for enhanced wheat productivity.

Combining QTL Mapping and Transcriptomics to Decipher the Genetic Architecture of Phenolic Compounds Metabolism in the Conifer White Spruce.

Conifer forests worldwide are becoming increasingly vulnerable to the effects of climate change. Although the production of phenolic compounds (PCs) has been shown to be modulated by biotic and abiotic stresses, the genetic basis underlying the variation in their constitutive production level remains poorly documented in conifers. We used QTL mapping and RNA-Seq to explore the complex polygenic network underlying the constitutive production of PCs in a white spruce (Picea glauca) full-sib family for 2 years. QTL detection was performed for nine PCs and differentially expressed genes (DEGs) were identified between individuals with high and low PC contents for five PCs exhibiting stable QTLs across time. A total of 17 QTLs were detected for eight metabolites, including one major QTL explaining up to 91.3% of the neolignan-2 variance. The RNA-Seq analysis highlighted 50 DEGs associated with phenylpropanoid biosynthesis, several key transcription factors, and a subset of 137 genes showing opposite expression patterns in individuals with high levels of the flavonoids gallocatechin and taxifolin glucoside. A total of 19 DEGs co-localized with QTLs. Our findings represent a significant step toward resolving the genomic architecture of PC production in spruce and facilitate the functional characterization of genes and transcriptional networks responsible for differences in constitutive production of PCs in conifers.

Legacy and emerging contaminants in common guillemot Uria aalge eggs in Ireland and Wales.

Guillemot eggs from multiple Irish colonies and one Welsh colony were analysed for legacy pollutants such as polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs) and other organochlorine compounds (OCs), as well as metals. Stable isotope ratios of carbon (δ13C) and nitrogen (δ15N) were measured in eggs to understand the influence of diet on contaminant levels detected. Wide-scope target and suspect screening techniques were used on a single guillemot egg, providing novel information on contaminants of emerging concern. Stable isotope ratio analysis showed that guillemots from Great Saltee Island and Lambay Island (Ireland's east coast) had a similar carbon source (δ13C) and fed at similar trophic levels (δ15N), pollutant levels were higher in eggs from Lambay Island near Dublin, Ireland's industrialised capital city. Guillemot eggs from Aughris Head (Atlantic west coast of Ireland), and Skomer Island (Wales) had differing isotopic niches to other colonies. Egg samples from Aughris Head had the lowest levels of pollutants in this study (with the exception of mercury) and amongst the lowest levels reported worldwide. In contrast, Skomer Island had the highest level of pollutants with higher concentrations of Σ16PCB, Σ6PBDE and HCB than Irish colonies, most likely a result of its proximity to historically industrial areas. Levels of PCBs, p,p' -DDE and mercury in guillemot eggs have decreased over time according to this study, in concurrence with worldwide trends. Levels of pollutants in guillemot eggs, in this study, fall below existing thresholds for adverse effects in other species, with the exception of mercury.