QTL-seq for the identification of candidate genes for days to flowering and leaf shape in pigeonpea

To identify genomic segments associated with days to flowering (DF) and leaf shape in pigeonpea, QTL-seq approach has been used in the present study. Genome-wide SNP profiling of extreme phenotypic bulks was conducted for both the traits from the segregating population (F2) derived from the cross combination- ICP 5529 × ICP 11605. A total of 126.63 million paired-end (PE) whole-genome resequencing data were generated for five samples, including one parent ICP 5529 (obcordate leaf and late-flowering plant), early and late flowering pools (EF and LF) and obcordate and lanceolate leaf shape pools (OLF and LLS). The QTL-seq identified two significant genomic regions, one on CcLG03 (1.58 Mb region spanned from 19.22 to 20.80 Mb interval) for days to flowering (LF and EF pools) and another on CcLG08 (2.19 Mb region spanned from 6.69 to 8.88 Mb interval) for OLF and LLF pools, respectively. Analysis of genomic regions associated SNPs with days to flowering and leaf shape revealed 5 genic SNPs present in the unique regions. The identified genomic regions for days to flowering were also validated with the genotyping-by-sequencing based classical QTL mapping method. A comparative analysis of the identified seven genes associated with days to flowering on 12 Fabaceae genomes, showed synteny with 9 genomes. A total of 153 genes were identified through the synteny analysis ranging from 13 to 36. This study demonstrates the usefulness of QTL-seq approach in precise identification of candidate gene(s) for days to flowering and leaf shape which can be deployed for pigeonpea improvement.

Blossom-end rot: a century-old problem in tomato (Solanum lycopersicum L.) and other vegetables

Blossom-end rot (BER) is a devastating physiological disorder affecting vegetable production worldwide. Extensive research into the physiological aspects of the disorder has demonstrated that the underlying causes of BER are associated with perturbed calcium (Ca2+) homeostasis and irregular watering conditions in predominantly cultivated accessions. Further, Reactive Oxygen Species (ROS) are critical players in BER development which, combined with unbalanced Ca2+ concentrations, greatly affect the severity of the disorder. The availability of a high-quality reference tomato genome as well as the whole genome resequencing of many accessions has recently permitted the genetic dissection of BER in segregating populations derived from crosses between cultivated tomato accessions. This has led to the identification of five loci contributing to BER from several studies. The eventual cloning of the genes contributing to BER would result in a deeper understanding of the molecular bases of the disorder. This will undoubtedly create crop improvement strategies for tomato as well as many other vegetables that suffer from BER.

Contrasting whole-genome and reduced representation sequencing for population demographic inference: an alpine mammal example

Genomic approaches to the study of population demography rely on accurate SNP calling and by-proxy the site frequency spectrum (SFS). Two main questions for the design of such studies remain poorly investigated: do reduced genomic sequencing summary statistics reflect that of whole genome, and how do sequencing strategies and derived summary statistics impact demographic inferences? To address those questions, we applied the ddRAD sequencing approach to 254 individuals and whole genome resequencing approach to 35 mountain goat (Oreamnos americanus) individuals across the species range with a known demographic history. We identified SNPs with 5 different variant callers and used ANGSD to estimate the genotype likelihoods (GLs). We tested combinations of SNP filtering by linkage disequilibrium (LD), minor allele frequency (MAF) and the genomic region. We compared the resulting suite of summary statistics reflective of the SFS and quantified the relationship to demographic inferences by estimating the contemporary effective population size (Ne), isolation-by-distance and population structure, FST, and explicit modelling of the demographic history with δaδi. Filtering had a larger effect than sequencing strategy, with the former strongly influencing summary statistics. Estimates of contemporary Ne and isolation-by-distance patterns were largely robust to the choice of sequencing, pipeline, and filtering. Despite the high variance in summary statistics, whole genome and reduced representation approaches were overall similar in supporting a glacial induced vicariance and low Ne in mountain goats. We discuss why whole genome resequencing data is preferable, and reiterate support the use of GLs, in part because it limits user-determined filters.

Construction and Genetic Analysis of a Worldwide Core Collection of Rapeseed

Rapeseed (Brassica napus) is an important oilseed crop, which is widely planted in the world. In a previous study, we collected 991 accessions of rapeseed from the worldwide germplasm and revealed genetic polymorphisms within these germplasm by whole-genome resequencing. However, management of such a large amount of accessions is time-consuming, laborious and costly. Therefore, we constructed a core collection of rapeseed consisting of 300 worldwide accessions based on their genetic diversity. Compared with 991 accessions, the worldwide core collection showed similar geographic distribution, the proportion of three ecotypes, nucleotide diversity and the associated SNPs of flowering time. Besides, we identified FT ortholog (BnaA02g12130D) and FLC ortholog (BnaA10g22080D) responsible for flowering time and ecotype differentiation through selective sweep analysis and genome-wide association analysis (GWAS) of flowering time using the rapeseed core collection. FT and FLC are two well-known genes regulating flowering time in Arabidopsis. These results indicate that the worldwide core collection can represent the genetic diversity of 991 worldwide accessions, which could be more efficiently used for phenotypic and genetic studies in rapeseed.

Whole-Genome Resequencing of Near-Isogenic Lines Reveals a Genomic Region Associated with High Trans-Lycopene Contents in Watermelon

Trans-lycopene is a functional phytochemical abundant in red-fleshed watermelons, and its contents vary among cultivars. In this study, the genetic basis of high trans-lycopene contents in scarlet red flesh was evaluated. Three near-isogenic lines (NILs) with high trans-lycopene contents were derived from the scarlet red-fleshed donor parent DRD and three coral red-fleshed (low trans-lycopene contents) recurrent parents. The lycopene contents of DRD (589.4 ± 71.8 µg/g) were two times higher than that of the recurrent parents, and values for NILs were intermediate between those of the parents. Coral red-fleshed lines and F1 cultivars showed low trans-lycopene contents (135.7 ± 18.0 µg/g to 213.7 ± 39.5 µg/g). Whole-genome resequencing of two NILs and their parents and an analysis of genome-wide single-nucleotide polymorphisms revealed three common introgressed regions (CIRs) on chromosomes 6, 9, and 10. Twenty-eight gene-based cleaved amplified polymorphic sequence (CAPS) markers were developed from the CIRs. The CAPS markers derived from CIR6 on chromosome 6, spanning approximately 1 Mb, were associated (R2 = 0.45–0.72) with the trans-lycopene contents, particularly CIR6-M1 and CIR6-M4. Our results imply that CIR6 is a major genomic region associated with variation in the trans-lycopene contents in red-fleshed watermelon, and CIR6-M1 and CIR6-M4 may be useful for marker-assisted selection.

SNP-Based Breeding for Broiler Resistance to Ascites and Evaluation of Correlated Production Traits

Background: The goal of this study was to evaluate marker-assisted selection (MAS) in broiler chickens using previously mapped gene regions associated with ascites syndrome incidence. The second-generation MAS products were assessed for impact on ascites phenotype and whether there were associated changes in important production traits. Previously, we used whole genome resequencing (WGR) to fine-map 28 chromosomal regions as associated with ascites phenotype in our experimental ascites broiler line (Relaxed, REL) based on a hypobaric chamber challenge. Genotypes for single nucleotide polymorphisms (SNPs) in mapped regions on chromosomes 2 and 22, were used for MAS in our REL line. After two generations, birds homozygous for the genotypes associated with resistance for both chromosomal regions were established. The MAS F 2 generation was then compared to the REL line for ascites susceptibility and 25 production traits. Results: Selection based on SNPs in the carboxypeptidase Q (CPQ, Gga2) and leucine rich repeat transmembrane neuronal 4 (LRRTM4, Gga22) gene regions resulted in a sex- and simulated altitude- dependent reduction of ascites incidence in two F 2 cohorts of the MAS line. Comparisons of the F 2 MAS and REL lines for production traits when reared at ambient pressure found no significant negative impacts for feed intake (FI), feed conversion ratio (FCR), or deboned part yields for either sex for two F 2 cohorts. There were, however, improvements in the MAS for full-trial body weight gain (BWG), FCR, absolute and relative tender weights, and relative drumstick weight. Conclusions: These results validate the mapping of the 28 chromosomal regions and demonstrate that fine mapping by WGR is an effective strategy for addressing a complex trait; it also stands as the first successful SNP-based selection program against a complex disease trait, such as ascites. The MAS line is comparable and, in some instances, superior, in growth performance to the REL control while being more resistant to ascites. This study indicates that MAS based on WGR can provide significant breeding potential in agricultural systems.

Whole genome resequencing and comparative genome analysis of three Puccinia striiformis f. sp. tritici pathotypes prevalent in India

Stripe rust disease of wheat, caused by Puccinia striiformis f. sp. tritici, ( Pst ) is one of the most serious diseases of wheat worldwide. In India, virulent stripe rust races have been constantly evolving in the North-Western Plains Zone leading to the failure of some of the most widely grown resistant varieties in the region. With the goal of studying the recent evolution of virulent races in this region, we conducted whole-genome sequencing of three prevalent Indian Pst pathotypes Pst46S119, Pst78S84 and Pst110S119. We assembled 58.62, 58.33 and 55.78 Mb of Pst110S119, Pst46S119 and Pst78S84 genome, respectively. Pathotypes were found to be highly heterozygous. Comparative phylogenetic analysis indicated the recent evolution of pathotypes Pst110S119 and Pst78S84 from Pst46S119. Pathogenicity-related genes classes (CAZyme, proteases, effectors, and secretome proteins) were identified and found to be under positive selection. Higher rate of gene family expansion was also observed in the three pathotypes. A strong association between the effector genes and transposable elements may be the source of the rapid evolution of these strains. Phylogenetic analysis differentiated the Indian races in this study from other known US, European, African and Asian races. Diagnostic markers developed for the identification of different Pst pathotypes will help tracking of yellow rust at farmers’ field and strategizing resistance gene deployment.

Whole-Genome Resequencing of Wild and Cultivated Cannabis Reveals the Genetic Structure and Adaptive Selection of Important Traits

Background: Cannabis is an important industrial crop, whose bast fiber, seed, flowers and leaves are widely used by humans，especially cannabinoids extracted from plants as medicine is a hot spot in recent years. China is one of the origins of cannabis, where it has been cultivated and utilized for more than 6000 years, with the largest planting area of industrial hemp at present. China is rich in cannabis germplasm resources covering different latitudes (23 to 51°N) and is one of the few countries with wild cannabis resources. However, the genetic structure of Chinese cannabis populations and the adaptive selection of important traits remain unclear.Results: We identified the main morphological and physiological characteristics of wild cannabis and defined the genetic structure and relationships among wild and cultivated Chinese cannabis accessions and foreign representatives. This suggested that wild resources in Xinjiang have played an important role in the process of cannabis domestication. Adaptive selection analysis revealed that cultivated cannabis has undergone selective evolution or adaptation in flowering, growth and stress tolerance, and many functional genes were identified. Flowering characteristics analysis implied that wild cannabis is native to high latitudes and possesses the typical characteristic of early flowering, while cultivated cannabis has undergone a process of adaptive evolution to adjust to natural photoperiod conditions in different latitudes through regulation of FT-like expression.Conclusion: This study clarifies the genetic structure of Chinese cannabis and provides insight into adaptive selection and breeding in cannabis.