Dissecting the genetic architecture of quantitative traits using genome-wide identity-by-descent sharing.

Additive and dominance genetic variances underlying the expression of quantitative traits are important quantities for predicting short-term responses to selection, but they are notoriously challenging to estimate in most non-model wild populations. Specifically, large-sized or panmictic populations may be characterized by low variance in genetic relatedness among individuals which, in turn, can prevent accurate estimation of quantitative genetic parameters. We used estimates of genome-wide identity-by-descent (IBD) sharing from autosomal SNP loci to estimate quantitative genetic parameters for ecologically important traits in nine-spined sticklebacks (Pungitius pungitius) from a large, outbred population. Using empirical and simulated datasets, with varying sample sizes and pedigree complexity, we assessed the performance of different crossing schemes in estimating additive genetic variance and heritability for all traits. We found that low variance in relatedness characteristic of wild outbred populations with high migration rate can impair the estimation of quantitative genetic parameters and bias heritability estimates downwards. On the other hand, the use of a half-sib/full-sib design allowed precise estimation of genetic variance components and revealed significant additive variance and heritability for all measured traits, with negligible dominance contributions. Genome-partitioning and QTL mapping analyses revealed that most traits had a polygenic basis and were controlled by genes at multiple chromosomes. Furthermore, different QTL contributed to variation in the same traits in different populations suggesting heterogeneous underpinnings of parallel evolution at the phenotypic level. Our results provide important guidelines for future studies aimed at estimating adaptive potential in the wild, particularly for those conducted in outbred large-sized populations.

Quantitative genomics-enabled selection for simultaneous improvement of lint yield and seed traits in cotton (Gossypium hirsutum L.).

KEY MESSAGE: A Bayesian linkage disequilibrium-based multiple-locus mixed model identified QTLs for fibre, seed and oil traits and predicted breeding worthiness of test lines, enabling their simultaneous improvement in cotton. Improving cotton seed and oil yields has become increasingly important while continuing to breed for higher lint yield. In this study, a novel Bayesian linkage disequilibrium-based multiple-locus mixed model was developed for QTL identification and genomic prediction (GP). A multi-parent population consisting of 256 recombinant inbred lines, derived from four elite cultivars with distinct combinations of traits, was used in the analysis of QTLs for lint percentage, seed index, lint index and seed oil content and their interrelations. All four traits were moderately heritable and correlated but with no large influence of genotype × environment interactions across multiple seasons. Seven to ten major QTLs were identified for each trait with many being adjacent or overlapping for different trait pairs. A fivefold cross-validation of the model indicated prediction accuracies of 0.46-0.62. GP results based on any two-season phenotypes were strongly correlated with phenotypic means of a pooled analysis of three-season experiments (r = 0.83-0.92). When used for selection of improvement in lint, seed and oil yields, GP captured 40-100% of individuals with comparable lint yields of those selected based on the three-season phenotypic results. Thus, this quantitative genomics-enabled approach can not only decipher the genomic variation underlying lint, seed and seed oil traits and their interrelations, but can provide predictions for their simultaneous improvement. We discuss future breeding strategies in cotton that will enhance the entire value of the crop, not just its fibre.

Faecalibacterium prausnitzii-derived extracellular vesicles alleviate chronic colitis-related intestinal fibrosis by macrophage metabolic reprogramming.

Faecalibacterium prausnitzii (F. prausnitzii) has been recognized for its various intestinal and extraintestinal benefits to human. And reduction of F. prausnitzii has been linked to an increased risk of intestinal fibrosis in patients of Crohn's disease (CD). In this study, oral administration of either live F. prausnitzii or its extracellular vesicles (FEVs) can markedly mitigate the severity of fibrosis in mice induced by repetitive administration of DSS. In vitro experiment revealed that FEVs were capable of directing the polarization of peripheral blood mononuclear cells (PBMCs) towards an M2b macrophage phenotype, which has been associated with anti-fibrotic activities. This effect of FEV was found to be stable under various conditions that promote the development of pro-fibrotic M1/M2a/M2c macrophages. Proteomics and RNA sequencing were performed to uncover the molecular modulation of macrophages by FEVs. Notably, we found that FEVs reprogramed every metabolism of macrophages by damaging the mitochondria, and inhibited oxidative phosphorylation and glycolysis. Moreover, FEV-treated macrophages showed a decreased expression of PPARγ and an altered lipid processing phenotype characterized by decreased cholesterol efflux, which may promote energy reprogramming. Taken together, these findings identify FEV as a driver of macrophage reprogramming, suggesting that triggering M2b macrophage polarization by oral admiration of FEV may serve as strategy to alleviate hyperfibrotic intestine conditions in CD.

Minocycline-loaded nHAP/PLGA microspheres for prevention of injury-related corneal angiogenesis.

BACKGROUND: Corneal neovascularization (CoNV) threatens vision by disrupting corneal avascularity, however, current treatments, including pharmacotherapy and surgery, are hindered by limitations in efficacy and adverse effects. Minocycline, known for its anti-inflammatory properties, could suppress CoNV but faces challenges in effective delivery due to the cornea's unique structure. Therefore, in this study a novel drug delivery system using minocycline-loaded nano-hydroxyapatite/poly (lactic-co-glycolic acid) (nHAP/PLGA) nanoparticles was developed to improve treatment outcomes for CoNV. RESULTS: Ultra-small nHAP was synthesized using high gravity technology, then encapsulated in PLGA by a double emulsion method to form nHAP/PLGA microspheres, attenuating the acidic by-products of PLGA degradation. The MINO@PLGA nanocomplex, featuring sustained release and permeation properties, demonstrated an efficient delivery system for minocycline that significantly inhibited the CoNV area in an alkali-burn model without exhibiting apparent cytotoxicity. On day 14, the in vivo microscope examination and ex vivo CD31 staining corroborated the inhibition of neovascularization, with the significantly smaller CoNV area (29.40% ± 6.55%) in the MINO@PLGA Tid group (three times daily) than that of the control group (86.81% ± 15.71%), the MINO group (72.42% ± 30.15%), and the PLGA group (86.87% ± 14.94%) (p < 0.05). Fluorescein sodium staining show MINO@PLGA treatments, administered once daily (Qd) and three times daily (Tid) demonstrated rapid corneal epithelial healing while the Alkali injury group and the DEX group showed longer healing times (p < 0.05). Additionally, compared to the control group, treatments with dexamethasone, MINO, and MINO@PLGA were associated with an increased expression of TGF-ß as evidenced by immunofluorescence, while the levels of pro-inflammatory cytokines IL-1ß and TNF-α demonstrated a significant decrease following alkali burn. Safety evaluations, including assessments of renal and hepatic biomarkers, along with H&E staining of major organs, revealed no significant cytotoxicity of the MINO@PLGA nanocomplex in vivo. CONCLUSIONS: The novel MINO@PLGA nanocomplex, comprising minocycline-loaded nHAP/PLGA microspheres, has shown a substantial capacity for preventing CoNV. This study confirms the complex's ability to downregulate inflammatory pathways, significantly reducing CoNV with minimal cytotoxicity and high biosafety in vivo. Given these findings, MINO@PLGA stands as a highly promising candidate for ocular conditions characterized by CoNV.

Genetic Mapping and Characterization of Verticillium Wilt Resistance in a Recombinant Inbred Population of Upland Cotton.

Verticillium wilt (VW) is an important and widespread disease of cotton and once established is long-lived and difficult to manage. In Australia, the non-defoliating pathotype of Verticillium dahliae is the most common, and extremely virulent. Breeding cotton varieties with increased VW resistance is the most economical and effective method of controlling this disease and is greatly aided by understanding the genetics of resistance. This study aimed to investigate VW resistance in 240 F7 recombinant inbred lines (RIL) derived from a cross between MCU-5, which has good resistance, and Siokra 1-4, which is susceptible. Using a controlled environment bioassay, we found that resistance based on plant survival or shoot biomass was complex but with major contributions from chromosomes D03 and D09, with genomic prediction analysis estimating a prediction accuracy of 0.73 based on survival scores compared to 0.36 for shoot biomass. Transcriptome analysis of MCU-5 and Siokra 1-4 roots uninfected or infected with V. dahliae revealed that the two cultivars displayed very different root transcriptomes and responded differently to V. dahliae infection. Ninety-nine differentially expressed genes were located in the two mapped resistance regions and so are potential candidates for further identifying the genes responsible for VW resistance.

Efficacy of cell-free DNA methylation-based blood test for colorectal cancer screening in high-risk population: a prospective cohort study.

BACKGROUND: Although colonoscopy is the standard screening test for colorectal cancer (CRC), its use is limited by a poor compliance rate, the need for extensive bowel preparation, and the risk of complications. As an alternative, an FDA-approved stool-based DNA test, Cologuard, has demonstrated satisfactory detection performance for CRC, but its compliance rate remains suboptimal, primarily attributable to individuals' reluctance to provide stool samples. METHODS: We developed a noninvasive blood-based CRC test, ColonSecure, based on cell-free DNA containing cancer-specific CpG island methylation patterns. We initially screened publicly available datasets for differentially methylated CpG sites in CRC with prediction potential. Subsequently, we performed two sequential bisulfite-free methylation sequencing on blood samples obtained from CRC patients and non-cancer controls. Through rigorous evaluation of each marker and machine learning-assisted feature selection, we identified 149 hypermethylated markers from over 193,000 CpG sites. These markers were then utilized to construct the ColonSecure model, enabling accurate CRC detection. RESULTS: We validated the efficacy of our cell-free DNA methylation-based blood test for CRC screening with 3493 high-risk individuals identified from 114,136 urban residents. The ColonSecure test identified 89 out of 103 CRC patients diagnosed by the follow-up colonoscopy, outperforming CEA, CRP, and CA19-9 (with a sensitivity of 86.4% compared to 45.6%, 39.8%, and 25.2% for CEA, CRP, and CA19-9 respectively; an AUROC of 0.956 compared to an AUROC of < 0.77 for other methods). CONCLUSION: Our observations emphasize the potential of our multiple cfDNA methylation marker-based test for CRC screening in high-risk populations.

Proteomics data analysis using multiple statistical approaches identified proteins and metabolic networks associated with sucrose accumulation in sugarcane.

BACKGROUND: Sugarcane is the most important sugar crop, contributing > 80% of global sugar production. High sucrose content is a key target of sugarcane breeding, yet sucrose improvement in sugarcane remains extremely slow for decades. Molecular breeding has the potential to break through the genetic bottleneck of sucrose improvement. Dissecting the molecular mechanism(s) and identifying the key genetic elements controlling sucrose accumulation will accelerate sucrose improvement by molecular breeding. In our previous work, a proteomics dataset based on 12 independent samples from high- and low-sugar genotypes treated with ethephon or water was established. However, in that study, employing conventional analysis, only 25 proteins involved in sugar metabolism were identified . RESULTS: In this work, the proteomics dataset used in our previous study was reanalyzed by three different statistical approaches, which include a logistic marginal regression, a penalized multiple logistic regression named Elastic net, as well as a Bayesian multiple logistic regression method named Stochastic search variable selection (SSVS) to identify more sugar metabolism-associated proteins. A total of 507 differentially abundant proteins (DAPs) were identified from this dataset, with 5 of them were validated by western blot. Among the DAPs, 49 proteins were found to participate in sugar metabolism-related processes including photosynthesis, carbon fixation as well as carbon, amino sugar, nucleotide sugar, starch and sucrose metabolism. Based on our studies, a putative network of key proteins regulating sucrose accumulation in sugarcane is proposed, with glucose-6-phosphate isomerase, 2-phospho-D-glycerate hydrolyase, malate dehydrogenase and phospho-glycerate kinase, as hub proteins. CONCLUSIONS: The sugar metabolism-related proteins identified in this work are potential candidates for sucrose improvement by molecular breeding. Further, this work provides an alternative solution for omics data processing.

Age-dependent genetic architecture across ontogeny of body size in sticklebacks.

Heritable variation in traits under natural selection is a prerequisite for evolutionary response. While it is recognized that trait heritability may vary spatially and temporally depending on which environmental conditions traits are expressed under, less is known about the possibility that genetic variance contributing to the expected selection response in a given trait may vary at different stages of ontogeny. Specifically, whether different loci underlie the expression of a trait throughout development and thus providing an additional source of variation for selection to act on in the wild, is unclear. Here we show that body size, an important life-history trait, is heritable throughout ontogeny in the nine-spined stickleback (Pungitius pungitius). Nevertheless, both analyses of quantitative trait loci and genetic correlations across ages show that different chromosomes/loci contribute to this heritability in different ontogenic time-points. This suggests that body size can respond to selection at different stages of ontogeny but that this response is determined by different loci at different points of development. Hence, our study provides important results regarding our understanding of the genetics of ontogeny and opens an interesting avenue of research for studying age-specific genetic architecture as a source of non-parallel evolution.

Mitochondrial heat-shock cognate protein 70 contributes to auxin-mediated embryo development.

In Arabidopsis thaliana, mitochondrial-localized heat-shock cognate protein 70-1 (mtHSC70-1) plays an important role in vegetativegrowth. However, whether mtHSC70-1 affects reproductive growth remains unknown. Here, we found that the mtHSC70-1 gene was expressed in the provascular cells of the embryo proper from the early heart stage onward during embryogenesis. Phenotypic analyses of mthsc70-1 mutants revealed that mtHSC70 deficiency leads to defective embryo development and that this effect is mediated by auxin. In addition to a dwarf phenotype, the mthsc70-1 mutant displayed defects in flower morphology, anther development, and embryogenesis. At early developmental stages, the mthsc70-1 embryos exhibited abnormal cell divisions in both embryo proper and suspensor cells. From heart stage onward, they displayed an abnormal shape such as with no or very small cotyledon protrusions, had aberrant number of cotyledons, or were twisted. These embryo defects were associated with reduced or ectopic expression of auxin responsive reporter DR5rev:GFP. Consistently, the expression of auxin biosynthesis and polar auxin transport genes were markedly altered in mthsc70-1. On the other hand, mitochondrial retrograde regulation (MRR) was enhanced in mthsc70-1. Treatment of wild-type plants with an inhibitor that activates mitochondrial retrograde signaling reduced the expression level of auxin biosynthesis and polar auxin transport genes and induced phenotypes similar to those of mthsc70-1. Taken together, our data reveal that loss of function of mtHSC70-1 induces MRR, which inhibits auxin biosynthesis and polar auxin transport, leading to abnormal auxin gradients and defective embryo development.

Genomic prediction of cotton fibre quality and yield traits using Bayesian regression methods.

Genomic selection or genomic prediction (GP) has increasingly become an important molecular breeding technology for crop improvement. GP aims to utilise genome-wide marker data to predict genomic breeding value for traits of economic importance. Though GP studies have been widely conducted in various crop species such as wheat and maize, its application in cotton, an essential renewable textile fibre crop, is still significantly underdeveloped. We aim to develop a new GP-based breeding system that can improve the efficiency of our cotton breeding program. This article presents a GP study on cotton fibre quality and yield traits using 1385 breeding lines from the Commonwealth Scientific and Industrial Research Organisation (CSIRO, Australia) cotton breeding program which were genotyped using a high-density SNP chip that generated 12,296 informative SNPs. The aim of this study was twofold: (1) to identify the models and data sources (i.e. genomic and pedigree) that produce the highest prediction accuracies; and (2) to assess the effectiveness of GP as a selection tool in the CSIRO cotton breeding program. The prediction analyses were conducted under various scenarios using different Bayesian predictive models. Results highlighted that the model combining genomic and pedigree information resulted in the best cross validated prediction accuracies: 0.76 for fibre length, 0.65 for fibre strength, and 0.64 for lint yield. Overall, this work represents the largest scale genomic selection studies based on cotton breeding trial data. Prediction accuracies reported in our study indicate the potential of GP as a breeding tool for cotton. The study highlighted the importance of incorporating pedigree and environmental factors in GP models to optimise the prediction performance.

How does bivalve size influence microplastics accumulation?

Microplastics (wasted plastic particles < 5 mm in diameter) are ubiquitously distributed in the marine environment. Filter-feeding and low trophic level bivalves are vulnerable to microplastics accumulation from the surrounding depositional environment, thereby threatening both ecological health and human food safety. Microplastics had been detected in lots of coastal Bivalvia species. However, the influence of biological morphology on the mechanism of microplastics accumulation is not clear. There is also a knowledge gap of which species are preferred for commercial consumption, which creates loopholes in risk identification for food safety. A survey on a commercial popular eaten but under-researched hard clam (Meretrix meretrix; Linnaeus, 1758) from a famous fishery port city in southern China was carried out to comprehensively analyze shell size influence on microplastics accumulation in bivalves and consequently, human intake risk via bivalve consumption. Detected microplastics count in per individual (MCI) was 24.64 ± 19.11 items · individual-1, and microplastics count per gram (MCG; wet weight with shell) was 0.66 ± 0.54 items · g-1. When the shell width grew by 1 mm, MCI increased by 1.01 times, but MCG decreased by 0.97 times. Dominant microplastics characteristics found in this study was fiber and fragment. Sizes ranged from 25 to 150 µm, and dark colors (black, red, and blue) were found. The mostly common polymers were polyethene (PE, 40%), polyethylene terephthalate (PET, 23%), and polypropylene (PP, 18%). Estimated annual intake (EAI) risk of microplastics via hard clam consumption by residents was 6652.26 ± 5327.28 items · year -1 · person -1. The microplastics in bivalves and EAI was relatively high. When shell width grew by 1 mm, EAI decreased by 0.97 times. Therefore, eating a fixed amount of larger hard clams with a relatively low amount of microplastics can reduce EAI risk for consumers. A systematic investigation of emission sources along main coast, where bivalve production is prominent will be useful for food safety control in this region.

Metabolic syndrome, metabolic comorbid conditions and risk of early-onset colorectal cancer.

OBJECTIVE: Factors that lead to metabolic dysregulation are associated with increased risk of early-onset colorectal cancer (CRC diagnosed under age 50). However, the association between metabolic syndrome (MetS) and early-onset CRC remains unexamined. DESIGN: We conducted a nested case-control study among participants aged 18-64 in the IBM MarketScan Commercial Database (2006-2015). Incident CRC was identified using pathologist-coded International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) codes, and controls were frequency matched. MetS was defined as presence of ≥3 conditions among obesity, hypertension, hyperlipidaemia and hyperglycaemia/type 2 diabetes, based on ICD-9-CM and use of medications. Multivariable logistic regressions were used to estimate ORs and 95% CIs. RESULTS: MetS was associated with increased risk of early-onset CRC (n=4673; multivariable adjusted OR 1.25; 95% CI 1.09 to 1.43), similar to CRC diagnosed at age 50-64 (n=14 928; OR 1.21; 95% CI 1.15 to 1.27). Compared with individuals without a metabolic comorbid condition, those with 1, 2 or ≥3 conditions had a 9% (1.09; 95% CI 1.00 to 1.17), 12% (1.12; 95% CI 1.01 to 1.24) and 31% (1.31; 95% CI 1.13 to 1.51) higher risk of early-onset CRC (ptrend <0.001). No associations were observed for one or two metabolic comorbid conditions and CRC diagnosed at age 50-64. These positive associations were driven by proximal (OR per condition 1.14; 95% CI 1.06 to 1.23) and distal colon cancer (OR 1.09; 95% CI 1.00 to 1.18), but not rectal cancer (OR 1.03; 95% CI 0.97 to 1.09). CONCLUSIONS: Metabolic dysregulation was associated with increased risk of early-onset CRC, driven by proximal and distal colon cancer, thus at least in part contribute to the rising incidence of early-onset CRC.

Genetic population structure constrains local adaptation in sticklebacks.

Repeated and independent adaptation to specific environmental conditions from standing genetic variation is common. However, if genetic variation is limited, the evolution of similar locally adapted traits may be restricted to genetically different and potentially less optimal solutions or prevented from happening altogether. Using a quantitative trait locus (QTL) mapping approach, we identified the genomic regions responsible for the repeated pelvic reduction (PR) in three crosses between nine-spined stickleback populations expressing full and reduced pelvic structures. In one cross, PR mapped to linkage group 7 (LG7) containing the gene Pitx1, known to control pelvic reduction also in the three-spined stickleback. In the two other crosses, PR was polygenic and attributed to 10 novel QTL, of which 90% were unique to specific crosses. When screening the genomes from 27 different populations for deletions in the Pitx1 regulatory element, these were only found in the population in which PR mapped to LG7, even though the morphological data indicated large-effect QTL for PR in several other populations as well. Consistent with the available theory and simulations parameterized on empirical data, we hypothesize that the observed variability in genetic architecture of PR is due to heterogeneity in the spatial distribution of standing genetic variation caused by >2× stronger population structuring among freshwater populations and >10× stronger genetic isolation by distance in the sea in nine-spined sticklebacks as compared to three-spined sticklebacks.

The BAG2 and BAG6 Genes Are Involved in Multiple Abiotic Stress Tolerances in Arabidopsis Thaliana.

The BAG proteins are a family of multi-functional co-chaperones. In plants, BAG proteins were found to play roles both in abiotic and biotic stress tolerance. However, the function of Arabidopsis BAG2 remains largely unknown, whereas BAG6 is required for plants' defense to pathogens, although it remains unknown whether BAG6 is involved in plants' tolerance to abiotic stresses. Here, we show that both BAG2 and BAG6 are expressed in various tissues and are upregulated by salt, mannitol, and heat treatments and by stress-related hormones including ABA, ethylene, and SA. Germination of bag2, bag6 and bag2 bag6 seeds is less sensitive to ABA compared to the wild type (WT), whereas BAG2 and BAG6 overexpression lines are hypersensitive to ABA. bag2, bag6, and bag2 bag6 plants show higher survival rates than WT in drought treatment but display lower survival rates in heat-stress treatment. Consistently, these mutants showed differential expression of several stress- and ABA-related genes such as RD29A, RD29B, NCED3 and ABI4 compared to the WT. Furthermore, these mutants exhibit lower levels of ROS after drought and ABA treatment but higher ROS accumulation after heat treatment than the WT. These results suggest that BAG2 and BAG6 are negatively involved in drought stress but play a positive role in heat stress in Arabidopsis.

Genome-wide association study identified novel candidate loci affecting wood formation in Norway spruce.

Norway spruce is a boreal forest tree species of significant ecological and economic importance. Hence there is a strong imperative to dissect the genetics underlying important wood quality traits in the species. We performed a functional genome-wide association study (GWAS) of 17 wood traits in Norway spruce using 178 101 single nucleotide polymorphisms (SNPs) generated from exome genotyping of 517 mother trees. The wood traits were defined using functional modelling of wood properties across annual growth rings. We applied a Least Absolute Shrinkage and Selection Operator (LASSO-based) association mapping method using a functional multilocus mapping approach that utilizes latent traits, with a stability selection probability method as the hypothesis testing approach to determine a significant quantitative trait locus. The analysis provided 52 significant SNPs from 39 candidate genes, including genes previously implicated in wood formation and tree growth in spruce and other species. Our study represents a multilocus GWAS for complex wood traits in Norway spruce. The results advance our understanding of the genetics influencing wood traits and identifies candidate genes for future functional studies.

A Gaussian process model and Bayesian variable selection for mapping function-valued quantitative traits with incomplete phenotypic data.

MOTIVATION: Recent advances in high dimensional phenotyping bring time as an extra dimension into the phenotypes. This promotes the quantitative trait locus (QTL) studies of function-valued traits such as those related to growth and development. Existing approaches for analyzing functional traits utilize either parametric methods or semi-parametric approaches based on splines and wavelets. However, very limited choices of software tools are currently available for practical implementation of functional QTL mapping and variable selection. RESULTS: We propose a Bayesian Gaussian process (GP) approach for functional QTL mapping. We use GPs to model the continuously varying coefficients which describe how the effects of molecular markers on the quantitative trait are changing over time. We use an efficient gradient based algorithm to estimate the tuning parameters of GPs. Notably, the GP approach is directly applicable to the incomplete datasets having even larger than 50% missing data rate (among phenotypes). We further develop a stepwise algorithm to search through the model space in terms of genetic variants, and use a minimal increase of Bayesian posterior probability as a stopping rule to focus on only a small set of putative QTL. We also discuss the connection between GP and penalized B-splines and wavelets. On two simulated and three real datasets, our GP approach demonstrates great flexibility for modeling different types of phenotypic trajectories with low computational cost. The proposed model selection approach finds the most likely QTL reliably in tested datasets. AVAILABILITY AND IMPLEMENTATION: Software and simulated data are available as a MATLAB package 'GPQTLmapping', and they can be downloaded from GitHub (https://github.com/jpvanhat/GPQTLmapping). Real datasets used in case studies are publicly available at QTL Archive. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The roles of climate, geography and natural selection as drivers of genetic and phenotypic differentiation in a widespread amphibian Hyla annectans (Anura: Hylidae).

The role of geological events and Pleistocene climatic fluctuations as drivers of current patterns of genetic variation in extant species has been a topic of continued interest among evolutionary biologists. Nevertheless, comprehensive studies of widely distributed species are still rare, especially from Asia. Using geographically extensive sampling of many individuals and a large number of nuclear single nucleotide polymorphisms (SNPs), we studied the phylogeography and historical demography of Hyla annectans populations in southern China. Thirty-five sampled populations were grouped into seven clearly defined genetic clusters that closely match phenotype-based subspecies classification. These lineages diverged 2.32-5.23 million years ago (Ma), a timing that closely aligns with the rapid and drastic uplifting of the Qinghai-Tibet Plateau and adjacent southwest China. Demographic analyses and species distribution models indicate that different populations of this species have responded differently to past climatic changes. In the Hengduan Mountains, most populations experienced a bottleneck, whereas the populations located outside of the Hengduan Mountains have gradually declined in size since the end of the last glaciation. In addition, the levels of phenotypic and genetic divergence were strongly correlated across major clades. These results highlight the combined effects of geological events and past climatic fluctuations, as well as natural selection, as drivers of contemporary patterns of genetic and phenotypic variation in a widely distributed anuran in Asia.

Population transcriptomics reveals weak parallel genetic basis in repeated marine and freshwater divergence in nine-spined sticklebacks.

The degree to which adaptation to similar selection pressures is underlain by parallel versus non-parallel genetic changes is a topic of broad interest in contemporary evolutionary biology. Sticklebacks provide opportunities to characterize and compare the genetic underpinnings of repeated marine-freshwater divergences at both intra- and interspecific levels. While the degree of genetic parallelism in repeated marine-freshwater divergences has been frequently studied in the three-spined stickleback (Gasterosteus aculeatus), much less is known about this in other stickleback species. Using a population transcriptomic approach, we identified both genetic and gene expression variations associated with marine-freshwater divergence in the nine-spined stickleback (Pungitius pungitius). Specifically, we used a genome-wide association study approach, and found that ~1% of the total 173,491 identified SNPs showed marine-freshwater ecotypic differentiation. A total of 861 genes were identified to have SNPs associated with marine-freshwater divergence, but only 12 of these genes have also been reported as candidates associated with marine-freshwater divergence in the three-spined stickleback. Hence, our results indicate a low degree of interspecific genetic parallelism in marine-freshwater divergence. Moreover, 1,578 genes in the brain and 1,050 genes in the liver were differentially expressed between marine and freshwater nine-spined sticklebacks, ~5% of which have also been identified as candidates associated with marine-freshwater divergence in the three-spined stickleback. However, only few of these (e.g., CLDND1) appear to have been involved in repeated marine-freshwater divergence in nine-spined sticklebacks. Taken together, the results indicate a low degree of genetic parallelism in repeated marine-freshwater divergence both at intra- and interspecific levels.

Analysis of phenotypic- and Estimated Breeding Values (EBV) to dissect the genetic architecture of complex traits in a Scots pine three-generation pedigree design.

In forest tree breeding, family-based Quantitative Trait Loci (QTL) studies are valuable as methods to dissect the complexity of a trait and as a source of candidate genes. In the field of conifer research, our study contributes to the evaluation of phenotypic and predicted breeding values for the identification of QTL linked to complex traits in a three-generation pedigree population in Scots pine (Pinus sylvestris L.). A total of 11 470 open pollinated F2-progeny trees established at three different locations, were measured for growth and adaptive traits. Breeding values were predicted for their 360 mothers, originating from a single cross of two grand-parents. A multilevel LASSO association analysis was conducted to detect QTL using genotypes of the mothers with the corresponding phenotypes and Estimated Breeding Values (EBV). Different levels of genotype-by-environment (G × E) effects among sites at different years, were detected for survival and height. Moderate-to-low narrow sense heritabilities and EBV accuracies were found for all traits and all sites. We identified 18 AFLPs and 12 SNPs to be associated with QTL for one or more traits. 62 QTL were significant with percentages of variance explained ranging from 1.7 to 18.9%. In those cases where the same marker was associated to a phenotypic or an ebvQTL, the ebvQTL always explained higher proportion of the variance, maybe due to the more accurate nature of Estimated Breeding Values (EBV). Two SNP-QTL showed pleiotropic effects for traits related with hardiness, seed, cone and flower production. Furthermore, we detected several QTL with significant effects across multiple ages, which could be considered as strong candidate loci for early selection. The lack of reproducibility of some QTL detected across sites may be due to environmental heterogeneity reflected by the genotype- and QTL-by-environment effects.