The genomics of linkage drag in inbred lines of sunflower.

Crop wild relatives represent valuable sources of alleles for crop improvement, including adaptation to climate change and emerging diseases. However, introgressions from wild relatives might have deleterious effects on desirable traits, including yield, due to linkage drag. Here, we analyzed the genomic and phenotypic impacts of wild introgressions in inbred lines of cultivated sunflower to estimate the impacts of linkage drag. First, we generated reference sequences for seven cultivated and one wild sunflower genotype, as well as improved assemblies for two additional cultivars. Next, relying on previously generated sequences from wild donor species, we identified introgressions in the cultivated reference sequences, as well as the sequence and structural variants they contain. We then used a ridge-regression best linear unbiased prediction (BLUP) model to test the effects of the introgressions on phenotypic traits in the cultivated sunflower association mapping population. We found that introgression has introduced substantial sequence and structural variation into the cultivated sunflower gene pool, including >3,000 new genes. While introgressions reduced genetic load at protein-coding sequences, they mostly had negative impacts on yield and quality traits. Introgressions found at high frequency in the cultivated gene pool had larger effects than low-frequency introgressions, suggesting that the former likely were targeted by artificial selection. Also, introgressions from more distantly related species were more likely to be maladaptive than those from the wild progenitor of cultivated sunflower. Thus, breeding efforts should focus, as far as possible, on closely related and fully compatible wild relatives.

Multi-genome comparisons reveal gain-and-loss evolution of anti-Mullerian hormone receptor type 2 as a candidate master sex-determining gene in Percidae.

BACKGROUND: The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii, and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. RESULTS: We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex-determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplicates (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been likely lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome 18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variations (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex-determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in the testis than in the ovary. CONCLUSIONS: Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

A cluster of putative resistance genes is associated with a dominant resistance to sunflower broomrape.

KEY MESSAGE: The HaOr5 resistance gene is located in a large genomic insertion containing putative resistance genes and provides resistance to O. cumana, preventing successful connection to the sunflower root vascular system. Orobanche cumana (sunflower broomrape) is a parasitic plant that is part of the Orobanchaceae family and specifically infests sunflower crops. This weed is an obligate parasitic plant that does not carry out photosynthetic activity or develop roots and is fully dependent on its host for its development. It produces thousands of dust-like seeds per plant. It possesses a high spreading ability and has been shown to quickly overcome resistance genes successively introduced by selection in cultivated sunflower varieties. The first part of its life cycle occurs underground. The connection to the sunflower vascular system is essential for parasitic plant survival and development. The HaOr5 gene provides resistance to sunflower broomrape race E by preventing the connection of O. cumana to the root vascular system. We mapped a single position of the HaOr5 gene by quantitative trait locus mapping using two segregating populations. The same location of the HaOr5 gene was identified by genome-wide association. Using a large population of thousands of F2 plants, we restricted the location of the HaOr5 gene to a genomic region of 193 kb. By sequencing the whole genome of the resistant line harboring the major resistance gene HaOr5, we identified a large insertion of a complex genomic region containing a cluster of putative resistance genes.

Whole-genome sequencing identifies interferon-induced protein IFI6/IFI27-like as a strong candidate gene for VNN resistance in European sea bass.

BACKGROUND: Viral nervous necrosis (VNN) is a major disease that affects European sea bass, and understanding the biological mechanisms that underlie VNN resistance is important for the welfare of farmed fish and sustainability of production systems. The aim of this study was to identify genomic regions and genes that are associated with VNN resistance in sea bass. RESULTS: We generated a dataset of 838,451 single nucleotide polymorphisms (SNPs) identified from whole-genome sequencing (WGS) in the parental generation of two commercial populations (A: 2371 individuals and B: 3428 individuals) of European sea bass with phenotypic records for binary survival in a VNN challenge. For each population, three cohorts were submitted to a red-spotted grouper nervous necrosis virus (RGNNV) challenge by immersion and genotyped on a 57K SNP chip. After imputation of WGS SNPs from their parents, quantitative trait loci (QTL) were mapped using a Bayesian sparse linear mixed model (BSLMM). We found several QTL regions that were specific to one of the populations on different linkage groups (LG), and one 127-kb QTL region on LG12 that was shared by both populations and included the genes ZDHHC14, which encodes a palmitoyltransferase, and IFI6/IFI27-like, which encodes an interferon-alpha induced protein. The most significant SNP in this QTL region was only 1.9 kb downstream of the coding sequence of the IFI6/IFI27-like gene. An unrelated population of four large families was used to validate the effect of the QTL. Survival rates of susceptible genotypes were 40.6% and 45.4% in populations A and B, respectively, while that of the resistant genotype was 66.2% in population B and 78% in population A. CONCLUSIONS: We have identified a genomic region that carries a major QTL for resistance to VNN and includes the ZDHHC14 and IFI6/IFI27-like genes. The potential involvement of the interferon pathway, a well-known anti-viral defense mechanism in several organisms (chicken, human, or fish), in survival to VNN infection is of particular interest. Our results can lead to major improvements for sea bass breeding programs through marker-assisted genomic selection to obtain more resistant fish.

SNP markers for early identification of high molecular weight glutenin subunits (HMW-GSs) in bread wheat.

KEY MESSAGE: A set of eight SNP markers was developed to facilitate the early selection of HMW-GS alleles in breeding programmes. In bread wheat (Triticum aestivum), the high molecular weight glutenin subunits (HMW-GSs) are the most important determinants of technological quality. Known to be very diverse, HMW-GSs are encoded by the tightly linked genes Glu-1-1 and Glu-1-2. Alleles that improve the quality of dough have been identified. Up to now, sodium dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE) of grain proteins is the most widely used for their identification. To facilitate the early selection of HMW-GS alleles in breeding programmes, we developed DNA-based molecular markers. For each accession of a core collection (n = 364 lines) representative of worldwide bread wheat diversity, HMW-GSs were characterized by both genotyping and SDS-PAGE. Based on electrophoresis, we observed at least 8, 22 and 9 different alleles at the Glu-A1, Glu-B1 and Glu-D1 loci, respectively, including new variants. We designed a set of 17 single-nucleotide polymorphism (SNP) markers that were representative of the most frequent SDS-PAGE alleles at each locus. At Glu-A1 and Glu-D1, two and three marker-based haplotypes, respectively, captured the diversity of the SDS-PAGE alleles rather well. Discrepancies were found mainly for the Glu-B1 locus. However, statistical tests revealed that two markers at each Glu-B1 gene and their corresponding haplotypes were more significantly associated with the rheological properties of the dough than were the relevant SDS-PAGE alleles. To conclude, this study demonstrates that the SNP markers developed provide additional information on HMW-GS diversity. Two markers at Glu-A1, four at Glu-B1 and two at Glu-D1 constitute a useful toolbox for breeding wheat to improve end-use value.

Amplifying recombination genome-wide and reshaping crossover landscapes in Brassicas.

Meiotic recombination by crossovers (COs) is tightly regulated, limiting its key role in producing genetic diversity. However, while COs are usually restricted in number and not homogenously distributed along chromosomes, we show here how to disrupt these rules in Brassica species by using allotriploid hybrids (AAC, 2n = 3x = 29), resulting from the cross between the allotetraploid rapeseed (B. napus, AACC, 2n = 4x = 38) and one of its diploid progenitors (B. rapa, AA, 2n = 2x = 20). We produced mapping populations from different genotypes of both diploid AA and triploid AAC hybrids, used as female and/or as male. Each population revealed nearly 3,000 COs that we studied with SNP markers well distributed along the A genome (on average 1 SNP per 1.25 Mbp). Compared to the case of diploids, allotriploid hybrids showed 1.7 to 3.4 times more overall COs depending on the sex of meiosis and the genetic background. Most surprisingly, we found that such a rise was always associated with (i) dramatic changes in the shape of recombination landscapes and (ii) a strong decrease of CO interference. Hybrids carrying an additional C genome exhibited COs all along the A chromosomes, even in the vicinity of centromeres that are deprived of COs in diploids as well as in most studied species. Moreover, in male allotriploid hybrids we found that Class I COs are mostly responsible for the changes of CO rates, landscapes and interference. These results offer the opportunity for geneticists and plant breeders to dramatically enhance the generation of diversity in Brassica species by disrupting the linkage drag coming from limits on number and distribution of COs.

Allele-specific expression and genetic determinants of transcriptomic variations in response to mild water deficit in tomato.

Characterizing the natural diversity of gene expression across environments is an important step in understanding how genotype-by-environment interactions shape phenotypes. Here, we analyzed the impact of water deficit onto gene expression levels in tomato at the genome-wide scale. We sequenced the transcriptome of growing leaves and fruit pericarps at cell expansion stage in a cherry and a large fruited accession and their F1 hybrid grown under two watering regimes. Gene expression levels were steadily affected by the genotype and the watering regime. Whereas phenotypes showed mostly additive inheritance, ~80% of the genes displayed non-additive inheritance. By comparing allele-specific expression (ASE) in the F1 hybrid to the allelic expression in both parental lines, respectively, 3005 genes in leaf and 2857 genes in fruit deviated from 1:1 ratio independently of the watering regime. Among these genes, ~55% were controlled by cis factors, ~25% by trans factors and ~20% by a combination of both types of factors. A total of 328 genes in leaf and 113 in fruit exhibited significant ASE-by-watering regime interaction, among which ~80% presented trans-by-watering regime interaction, suggesting a response to water deficit mediated through a majority of trans-acting loci in tomato. We cross-validated the expression levels of 274 transcripts in fruit and leaves of 124 recombinant inbred lines (RILs) and identified 163 expression quantitative trait loci (eQTLs) mostly confirming the divergences identified by ASE. Combining phenotypic and expression data, we observed a complex network of variation between genes encoding enzymes involved in the sugar metabolism.

Genetic architecture of a plant adaptive trait: QTL mapping of intraspecific variation for tolerance to metal pollution in Arabidopsis halleri.

Anthropogenic activities are among the main drivers of global change and result in drastic habitat modifications, which represent strong evolutionary challenges for biological species that can either migrate, adapt, or disappear. In this context, understanding the genetics of adaptive traits is a prerequisite to enable long-term maintenance of populations under strong environmental constraints. To examine these processes, a QTL approach was developed here using the pseudometallophyte Arabidopsis halleri, which displays among-population adaptive divergence for tolerance to metallic pollution in soils. An F2 progeny was obtained by crossing individuals from metallicolous and non-metallicolous populations from Italian Alps, where intense metallurgic activities have created strong landscape heterogeneity. Then, we combined genome de novo assembly and genome resequencing of parental genotypes to obtain single-nucleotide polymorphism markers and achieve high-throughput genotyping of the progeny. QTL analysis was performed using growth parameters and photosynthetic yield to assess zinc tolerance levels. One major QTL was identified for photosynthetic yield. It explained about 27% of the phenotypic variance. Functional annotation of the QTL and gene expression analyses highlighted putative candidate genes. Our study represents a successful approach combining evolutionary genetics and advanced molecular tools, helping to better understand how a species can face new selective pressures of anthropogenic origin.

Genome-wide estimates of genetic diversity, inbreeding and effective size of experimental and commercial rainbow trout lines undergoing selective breeding.

BACKGROUND: Selective breeding is a relatively recent practice in aquaculture species compared to terrestrial livestock. Nevertheless, the genetic variability of farmed salmonid lines, which have been selected for several generations, should be assessed. Indeed, a significant decrease in genetic variability due to high selection intensity could have occurred, potentially jeopardizing the long-term genetic progress as well as the adaptive capacities of populations facing change(s) in the environment. Thus, it is important to evaluate the impact of selection practices on genetic diversity to limit future inbreeding. The current study presents an analysis of genetic diversity within and between six French rainbow trout (Oncorhynchus mykiss) experimental or commercial lines based on a medium-density single nucleotide polymorphism (SNP) chip and various molecular genetic indicators: fixation index (FST), linkage disequilibrium (LD), effective population size (Ne) and inbreeding coefficient derived from runs of homozygosity (ROH). RESULTS: Our results showed a moderate level of genetic differentiation between selected lines (FST ranging from 0.08 to 0.15). LD declined rapidly over the first 100 kb, but then remained quite high at long distances, leading to low estimates of Ne in the last generation ranging from 24 to 68 depending on the line and methodology considered. These results were consistent with inbreeding estimates that varied from 10.0% in an unselected experimental line to 19.5% in a commercial line, and which are clearly higher than corresponding estimates in ruminants or pigs. In addition, strong variations in LD and inbreeding were observed along the genome that may be due to differences in local rates of recombination or due to key genes that tended to have fixed favorable alleles for domestication or production. CONCLUSIONS: This is the first report on ROH for any aquaculture species. Inbreeding appeared to be moderate to high in the six French rainbow trout lines, due to founder effects at the start of the breeding programs, but also likely to sweepstakes reproductive success in addition to selection for the selected lines. Efficient management of inbreeding is a major goal in breeding programs to ensure that populations can adapt to future breeding objectives and SNP information can be used to manage the rate at which inbreeding builds up in the fish genome.

Development and validation of the Axiom(®) Apple480K SNP genotyping array.

Cultivated apple (Malus × domestica Borkh.) is one of the most important fruit crops in temperate regions, and has great economic and cultural value. The apple genome is highly heterozygous and has undergone a recent duplication which, combined with a rapid linkage disequilibrium decay, makes it difficult to perform genome-wide association (GWA) studies. Single nucleotide polymorphism arrays offer highly multiplexed assays at a relatively low cost per data point and can be a valid tool for the identification of the markers associated with traits of interest. Here, we describe the development and validation of a 487K SNP Affymetrix Axiom(®) genotyping array for apple and discuss its potential applications. The array has been built from the high-depth resequencing of 63 different cultivars covering most of the genetic diversity in cultivated apple. The SNPs were chosen by applying a focal points approach to enrich genic regions, but also to reach a uniform coverage of non-genic regions. A total of 1324 apple accessions, including the 92 progenies of two mapping populations, have been genotyped with the Axiom(®) Apple480K to assess the effectiveness of the array. A large majority of SNPs (359 994 or 74%) fell in the stringent class of poly high resolution polymorphisms. We also devised a filtering procedure to identify a subset of 275K very robust markers that can be safely used for germplasm surveys in apple. The Axiom(®) Apple480K has now been commercially released both for public and proprietary use and will likely be a reference tool for GWA studies in apple.

Analysis of the genetic diversity and structure across a wide range of germplasm reveals prominent gene flow in apple at the European level.

BACKGROUND: The amount and structure of genetic diversity in dessert apple germplasm conserved at a European level is mostly unknown, since all diversity studies conducted in Europe until now have been performed on regional or national collections. Here, we applied a common set of 16 SSR markers to genotype more than 2,400 accessions across 14 collections representing three broad European geographic regions (North + East, West and South) with the aim to analyze the extent, distribution and structure of variation in the apple genetic resources in Europe. RESULTS: A Bayesian model-based clustering approach showed that diversity was organized in three groups, although these were only moderately differentiated (FST = 0.031). A nested Bayesian clustering approach allowed identification of subgroups which revealed internal patterns of substructure within the groups, allowing a finer delineation of the variation into eight subgroups (FST = 0.044). The first level of stratification revealed an asymmetric division of the germplasm among the three groups, and a clear association was found with the geographical regions of origin of the cultivars. The substructure revealed clear partitioning of genetic groups among countries, but also interesting associations between subgroups and breeding purposes of recent cultivars or particular usage such as cider production. Additional parentage analyses allowed us to identify both putative parents of more than 40 old and/or local cultivars giving interesting insights in the pedigree of some emblematic cultivars. CONCLUSIONS: The variation found at group and subgroup levels may reflect a combination of historical processes of migration/selection and adaptive factors to diverse agricultural environments that, together with genetic drift, have resulted in extensive genetic variation but limited population structure. The European dessert apple germplasm represents an important source of genetic diversity with a strong historical and patrimonial value. The present work thus constitutes a decisive step in the field of conservation genetics. Moreover, the obtained data can be used for defining a European apple core collection useful for further identification of genomic regions associated with commercially important horticultural traits in apple through genome-wide association studies.

FRIZZY PANICLE drives supernumerary spikelets in bread wheat.

Bread wheat (Triticum aestivum) inflorescences, or spikes, are characteristically unbranched and normally bear one spikelet per rachis node. Wheat mutants on which supernumerary spikelets (SSs) develop are particularly useful resources for work towards understanding the genetic mechanisms underlying wheat inflorescence architecture and, ultimately, yield components. Here, we report the characterization of genetically unrelated mutants leading to the identification of the wheat FRIZZY PANICLE (FZP) gene, encoding a member of the APETALA2/Ethylene Response Factor transcription factor family, which drives the SS trait in bread wheat. Structural and functional characterization of the three wheat FZP homoeologous genes (WFZP) revealed that coding mutations of WFZP-D cause the SS phenotype, with the most severe effect when WFZP-D lesions are combined with a frameshift mutation in WFZP-A. We provide WFZP-based resources that may be useful for genetic manipulations with the aim of improving bread wheat yield by increasing grain number.

Multi-genome comparisons reveal gain-and-loss evolution of the anti-Mullerian hormone receptor type 2 gene, an old master sex determining gene, in Percidae.

The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplications (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome-18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variants (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in testis than ovary. Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Deciphering the genomic structure, function and evolution of carotenogenesis related phytoene synthases in grasses.

BACKGROUND: Carotenoids are isoprenoid pigments, essential for photosynthesis and photoprotection in plants. The enzyme phytoene synthase (PSY) plays an essential role in mediating condensation of two geranylgeranyl diphosphate molecules, the first committed step in carotenogenesis. PSY are nuclear enzymes encoded by a small gene family consisting of three paralogous genes (PSY1-3) that have been widely characterized in rice, maize and sorghum. RESULTS: In wheat, for which yellow pigment content is extremely important for flour colour, only PSY1 has been extensively studied because of its association with QTLs reported for yellow pigment whereas PSY2 has been partially characterized. Here, we report the isolation of bread wheat PSY3 genes from a Renan BAC library using Brachypodium as a model genome for the Triticeae to develop Conserved Orthologous Set markers prior to gene cloning and sequencing. Wheat PSY3 homoeologous genes were sequenced and annotated, unravelling their novel structure associated with intron-loss events and consequent exonic fusions. A wheat PSY3 promoter region was also investigated for the presence of cis-acting elements involved in the response to abscisic acid (ABA), since carotenoids also play an important role as precursors of signalling molecules devoted to plant development and biotic/abiotic stress responses. Expression of wheat PSYs in leaves and roots was investigated during ABA treatment to confirm the up-regulation of PSY3 during abiotic stress. CONCLUSIONS: We investigated the structural and functional determinisms of PSY genes in wheat. More generally, among eudicots and monocots, the PSY gene family was found to be associated with differences in gene copy numbers, allowing us to propose an evolutionary model for the entire PSY gene family in Grasses.

Isolation and characterisation of 11 polymorphic microsatellite markers in Papaver rhoeas L. (Corn Poppy), a major annual plant species from cultivated areas.

Papaver rhoeas, an annual plant species in the Papaveraceae family, is part of the biodiversity of agricultural ecosystems and also a noxious agronomic weed. We developed microsatellite markers to study the genetic diversity of P. rhoeas, using an enriched microsatellite library coupled with 454 next-generation sequencing. A total of 13,825 sequences were obtained that yielded 1795 microsatellite loci. After discarding loci with less than six repeats of the microsatellite motif, automated primer design was successful for 598 loci. We tested 74 of these loci for amplification with a total of 97 primer pairs. Thirty loci passed our tests and were subsequently tested for polymorphism using 384 P. rhoeas plants originating from 12 populations from France. Of the 30 loci, 11 showed reliable polymorphism not affected by the presence of null alleles. The number of alleles and the expected heterozygosity ranged from 3 to 7.4 and from 0.27 to 0.73, respectively. A low but significant genetic differentiation among populations was observed (F(ST) = 0.04; p < 0.001). The 11 validated polymorphic microsatellite markers developed in this work will be useful in studies of genetic diversity and population structure of P. rhoeas, assisting in designing management strategies for the control or the conservation of this species.

Development of a High-Density 665 K SNP Array for Rainbow Trout Genome-Wide Genotyping.

Single nucleotide polymorphism (SNP) arrays, also named « SNP chips ¼, enable very large numbers of individuals to be genotyped at a targeted set of thousands of genome-wide identified markers. We used preexisting variant datasets from USDA, a French commercial line and 30X-coverage whole genome sequencing of INRAE isogenic lines to develop an Affymetrix 665 K SNP array (HD chip) for rainbow trout. In total, we identified 32,372,492 SNPs that were polymorphic in the USDA or INRAE databases. A subset of identified SNPs were selected for inclusion on the chip, prioritizing SNPs whose flanking sequence uniquely aligned to the Swanson reference genome, with homogenous repartition over the genome and the highest Minimum Allele Frequency in both USDA and French databases. Of the 664,531 SNPs which passed the Affymetrix quality filters and were manufactured on the HD chip, 65.3% and 60.9% passed filtering metrics and were polymorphic in two other distinct French commercial populations in which, respectively, 288 and 175 sampled fish were genotyped. Only 576,118 SNPs mapped uniquely on both Swanson and Arlee reference genomes, and 12,071 SNPs did not map at all on the Arlee reference genome. Among those 576,118 SNPs, 38,948 SNPs were kept from the commercially available medium-density 57 K SNP chip. We demonstrate the utility of the HD chip by describing the high rates of linkage disequilibrium at 2-10 kb in the rainbow trout genome in comparison to the linkage disequilibrium observed at 50-100 kb which are usual distances between markers of the medium-density chip.

Cantaloupe melon genome reveals 3D chromatin features and structural relationship with the ancestral cucurbitaceae karyotype.

Cucumis melo displays a large diversity of horticultural groups with cantaloupe melon the most cultivated type. Using a combination of single-molecule sequencing, 10X Genomics link-reads, high-density optical and genetic maps, and chromosome conformation capture (Hi-C), we assembled a chromosome scale C. melo var. cantalupensis Charentais mono genome. Integration of RNA-seq, MeDip-seq, ChIP-seq, and Hi-C data revealed a widespread compartmentalization of the melon genome, segregating constitutive heterochromatin and euchromatin. Genome-wide comparative and evolutionary analysis between melon botanical groups identified Charentais mono genome increasingly more divergent from Harukei-3 (reticulatus), Payzawat (inodorus), and HS (ssp. agrestis) genomes. To assess the paleohistory of the Cucurbitaceae, we reconstructed the ancestral Cucurbitaceae karyotype and compared it to sequenced cucurbit genomes. In contrast to other species that experienced massive chromosome shuffling, melon has retained the ancestral genome structure. We provide comprehensive genomic resources and new insights in the diversity of melon horticultural groups and evolution of cucurbits.

An ancient truncated duplication of the anti-Müllerian hormone receptor type 2 gene is a potential conserved master sex determinant in the Pangasiidae catfish family.

The evolution of sex determination (SD) in teleosts is amazingly dynamic, as reflected by the variety of different master sex-determining genes identified. Pangasiids are economically important catfishes in South Asian countries, but little is known about their SD system. Here, we generated novel genomic resources for 12 Pangasiids and characterized their SD system. Based on a Pangasianodon hypophthalmus chromosome-scale genome assembly, we identified an anti-Müllerian hormone receptor type Ⅱ gene (amhr2) duplication, which was further characterized as being sex-linked in males and expressed only in testes. These results point to a Y chromosome male-specific duplication (amhr2by) of the autosomal amhr2a. Sequence annotation revealed that the P. hypophthalmus Amhr2by is truncated in its N-terminal domain, lacking the cysteine-rich extracellular part of the receptor that is crucial for ligand binding, suggesting a potential route for its neofunctionalization. Reference-guided assembly of 11 additional Pangasiids, along with sex-linkage studies, revealed that this truncated amhr2by duplication is a male-specific conserved gene in Pangasiids. Reconstructions of the amhr2 phylogeny suggested that amhr2by arose from an ancient duplication/insertion event at the root of the Siluroidei radiation that is dated to ~100 million years ago. Together these results bring multiple lines of evidence supporting that amhr2by is an ancient and conserved master sex-determining gene in Pangasiids, a finding that highlights the recurrent use of the transforming growth factor ß pathway, which is often used for the recruitment of teleost master SD genes, and provides another empirical case towards firther understanding of dynamics of SD systems.

Development of microsatellite markers in Capsella rubella and Capsella bursa-pastoris (Brassicaceae).

PREMISE OF THE STUDY: We developed microsatellite markers to investigate genetic diversity within and among populations of Capsella rubella and Capsella bursa-pastoris and between these two species. METHODS AND RESULTS: Fourteen polymorphic microsatellite loci were identified in the two species and one more polymorphic microsatellite locus only in C. rubella. Samples from different European localities were genotyped. Up to six alleles per locus were observed in C. rubella, and up to 22 alleles per locus in C. bursa-pastoris. Observed heterozygosities were low, indicating high selfing rates, especially in C. rubella. CONCLUSIONS: The results provide valuable information on genetic diversity for future studies of population genetics in C. rubella and C. bursa-pastoris.

Harnessing the power of next-generation sequencing technologies to the purpose of high-throughput pesticide resistance diagnosis.

BACKGROUND: Next Generation Sequencing (NGS) technologies offer tremendous possibilities for high-throughput pesticide resistance diagnosis via massive genotyping-by-sequencing. Herein, we used Illumina sequencing combined with a simple, non-commercial bioinformatics pipe-line to seek mutations involved in herbicide resistance in two weeds. RESULTS: DNA was extracted from 96 pools of 50 plants for each species. Three amplicons encompassing 15 ALS (acetolactate-synthase) codons crucial for herbicide resistance were amplified from each DNA extract. Above 18 and 20 million quality 250-nucleotide sequence reads were obtained for groundsel (Senecio vulgaris, tetraploid) and ragweed (Ambrosia artemisiifolia, diploid), respectively. Herbicide resistance-endowing mutations were identified in 45 groundsel and in eight ragweed field populations. The mutations detected and their frequencies assessed by NGS were checked by individual plant genotyping or Sanger sequencing. NGS results were fully confirmed, except in three instances out of 12 where mutations present at a frequency of 1% were detected below the threshold set for reliable mutation detection. CONCLUSION: Analyzing 9600 plants requested 192 DNA extractions followed by 1728 PCRs and two Illumina runs. Equivalent results obtained by individual analysis would have necessitated 9600 individual DNA extractions followed by 216 000 genotyping PCRs, or by 121 500 PCRs and 40 500 Sanger sequence runs. This clearly demonstrates the interest and power of NGS-based detection of pesticide resistance from pools of individuals for diagnosing resistance in massive numbers of individuals. © 2019 Society of Chemical Industry.