Construction of a high-density DArTseq SNP-based genetic map and identification of genomic regions with segregation distortion in a genetic population derived from a cross between feral and cultivated-type watermelon

Quantitative Trait Loci Mapping and Comparative Transcriptome Analysis of Fruit Weight (FW) in Watermelon (Citrullus lanatus L.)

The watermelon (Citrullus lanatus L.) holds substantial economic value as a globally cultivated horticultural crop. However, the genetic architecture of watermelon fruit weight (FW) remains poorly understood. In this study, we used sh14-11 with small fruit and N14 with big fruit to construct 100 recombinant inbred lines (RILs). Based on whole-genome resequencing (WGR), 218,127 single nucleotide polymorphisms (SNPs) were detected to construct a high-quality genetic map. After quantitative trait loci (QTL) mapping, a candidate interval of 31-38 Mb on chromosome 2 was identified for FW. Simultaneously, the bulked segregant analysis (BSA) in the F2 population corroborated the identification of the same interval, encompassing the homologous gene linked to the known FW-related gene fas. Additionally, RNA-seq was carried out across 11 tissues from sh14-11 and N14, revealing expression profiles that identified 1695 new genes and corrected the annotation of 2941 genes. Subsequent differential expression analysis unveiled 8969 differentially expressed genes (DEGs), with 354 of these genes exhibiting significant differences across four key developmental stages. The integration of QTL mapping and differential expression analysis facilitated the identification of 14 FW-related genes, including annotated TGA and NAC transcription factors implicated in fruit development. This combined approach offers valuable insights into the genetic basis of FW, providing crucial resources for enhancing watermelon cultivation.

Preserving Pure Siamese Crocodile Populations: A Comprehensive Approach Using Multi-Genetic Tools

Hybrids between the critically endangered Siamese crocodile (Crocodylus siamensis) and least-concern saltwater crocodile (C. porosus) in captive populations represent a serious challenge for conservation and reintroduction programs due to the impact of anthropogenic activities. A previous study used microsatellite and mitochondrial DNA data to establish the criteria for identifying species and their hybrids; however, the results may have been influenced by biased allelic frequencies and genetic drift within the examined population. To overcome these limitations and identify the true signals of selection, alternative DNA markers and a diverse set of populations should be employed. Therefore, this study used DArT sequencing to identify genome-wide single nucleotide polymorphisms (SNPs) in both species and confirm the genetic scenario of the parental species and their hybrids. A population of saltwater crocodiles from Australia was used to compare the distribution of species-diagnostic SNPs. Different analytical approaches were compared to diagnose the level of hybridization when an admixture was present, wherein three individuals had potential backcrossing. Approximately 17.00-26.00% of loci were conserved between the Siamese and saltwater crocodile genomes. Species-diagnostic SNP loci for Siamese and saltwater crocodiles were identified as 8051 loci and 1288 loci, respectively. To validate the species-diagnostic SNP loci, a PCR-based approach was used by selecting 20 SNP loci for PCR primer design, among which 3 loci were successfully able to differentiate the actual species and different hybridization levels. Mitochondrial and nuclear genetic information, including microsatellite genotyping and species-diagnostic DNA markers, were combined as a novel method that can compensate for the limitations of each method. This method enables conservation prioritization before release into the wild, thereby ensuring sustainable genetic integrity for long-term species survival through reintroduction and management programs.

Widespread genetic heterogeneity and genotypic grouping associated with fungicide resistance among barley spot form net blotch isolates in Australia

Spot form net blotch, caused by Pyrenophora teres f. maculata, is a major foliar disease of barley worldwide. Knowledge of the pathogen's genetic diversity and population structure is critical for a better understanding of inherent evolutionary capacity and for the development of sustainable disease management strategies. Genome-wide, single nucleotide polymorphism data of 254 Australian isolates revealed genotypic diversity and an absence of population structure, either between states, or between fields and cultivars in different agro-ecological zones. This indicates there is little geographical isolation or cultivar directional selection and that the pathogen is highly mobile across the continent. However, two cryptic genotypic groups were found only in Western Australia, predominantly associated with genes involved in fungicide resistance. The findings in this study are discussed in the context of current cultivar resistance and the pathogen's adaptive potential.

Development of a new AgriSeq 4K mid-density SNP genotyping panel and its utility in pearl millet breeding

Pearl millet is a crucial nutrient-rich staple food in Asia and Africa and adapted to the climate of semi-arid topics. Since the genomic resources in pearl millet are very limited, we have developed a brand-new mid-density 4K SNP panel and demonstrated its utility in genetic studies. A set of 4K SNPs were mined from 925 whole-genome sequences through a comprehensive in-silico pipeline. Three hundred and seventy-three genetically diverse pearl millet inbreds were genotyped using the newly-developed 4K SNPs through the AgriSeq Targeted Genotyping by Sequencing technology. The 4K SNPs were uniformly distributed across the pearl millet genome and showed considerable polymorphism information content (0.23), genetic diversity (0.29), expected heterozygosity (0.29), and observed heterozygosity (0.03). The SNP panel successfully differentiated the accessions into two major groups, namely B and R lines, through genetic diversity, PCA, and structure models as per their pedigree. The linkage disequilibrium (LD) analysis showed Chr3 had higher LD regions while Chr1 and Chr2 had more low LD regions. The genetic divergence between the B- and R-line populations was 13%, and within the sub-population variability was 87%. In this experiment, we have mined 4K SNPs and optimized the genotyping protocol through AgriSeq technology for routine use, which is cost-effective, fast, and highly reproducible. The newly developed 4K mid-density SNP panel will be useful in genomics and molecular breeding experiments such as assessing the genetic diversity, trait mapping, backcross breeding, and genomic selection in pearl millet.

Identification of ancestral sex chromosomes in the frog Glandirana rugosa bearing XX-XY and ZZ-ZW sex-determining systems

Sex chromosomes constantly exist in a dynamic state of evolution: rapid turnover and change of heterogametic sex during homomorphic state, and often stepping out to a heteromorphic state followed by chromosomal decaying. However, the forces driving these different trajectories of sex chromosome evolution are still unclear. The Japanese frog Glandirana rugosa is one taxon well suited to the study on these driving forces. The species has two different heteromorphic sex chromosome systems, XX-XY and ZZ-ZW, which are separated in different geographic populations. Both XX-XY and ZZ-ZW sex chromosomes are represented by chromosome 7 (2n = 26). Phylogenetically, these two systems arose via hybridization between two ancestral lineages of West Japan and East Japan populations, of which sex chromosomes are homomorphic in both sexes and to date have not yet been identified. Identification of the sex chromosomes will give us important insight into the mechanisms of sex chromosome evolution in this species. Here, we used a high-throughput genomic approach to identify the homomorphic XX-XY sex chromosomes in both ancestral populations. Sex-linked DNA markers of West Japan were aligned to chromosome 1, whereas those of East Japan were aligned to chromosome 3. These results reveal that at least two turnovers across three different sex chromosomes 1, 3 and 7 occurred during evolution of this species. This finding raises the possibility that cohabitation of the two different sex chromosomes from ancestral lineages induced turnover to another new one in their hybrids, involving transition of heterogametic sex and evolution from homomorphy to heteromorphy.

Something Fishy about Siamese Fighting Fish (Betta splendens) Sex: Polygenic Sex Determination or a Newly Emerged Sex-Determining Region?

Fishes provide a unique and intriguing model system for studying the genomic origin and evolutionary mechanisms underlying sex determination and high sex-chromosome turnover. In this study, the mode of sex determination was investigated in Siamese fighting fish, a species of commercial importance. Genome-wide SNP analyses were performed on 75 individuals (40 males and 35 females) across commercial populations to determine candidate sex-specific/sex-linked loci. In total, 73 male-specific loci were identified and mapped to a 5.6 kb region on chromosome 9, suggesting a putative male-determining region (pMDR) containing localized dmrt1 and znrf3 functional sex developmental genes. Repeat annotations of the pMDR revealed an abundance of transposable elements, particularly Ty3/Gypsy and novel repeats. Remarkably, two out of the 73 male-specific loci were located on chromosomes 7 and 19, implying the existence of polygenic sex determination. Besides male-specific loci, five female-specific loci on chromosome 9 were also observed in certain populations, indicating the possibility of a female-determining region and the polygenic nature of sex determination. An alternative explanation is that male-specific loci derived from other chromosomes or female-specific loci in Siamese fighting fish recently emerged as new sex-determining loci during domestication and repeated hybridization.

Genome-Wide SNP Analysis of Hybrid Clariid Fish Reflects the Existence of Polygenic Sex-Determination in the Lineage

The African catfish (Clarias gariepinus) may exhibit the co-existence of XX/XY and ZZ/ZW sex-determination systems (SDSs). However, the SDS of African catfish might be influenced by a polygenic sex-determination (PSD) system, comprising multiple independently segregating sex “switch” loci to determine sex within a species. Here, we aimed to detect the existence of PSD using hybrid. The hybrid produced by crossing male African catfish with female bighead catfish (C. macrocephalus, XX/XY) is a good animal model to study SDSs. Determining the SDS of hybrid catfish can help in understanding the interactions between these two complex SDS systems. Using the genotyping-by-sequencing “DART-seq” approach, we detected seven moderately male-linked loci and seventeen female-linked loci across all the examined hybrid specimens. Most of these loci were not sex-linked in the parental species, suggesting that the hybrid exhibits a combination of different alleles. Annotation of the identified sex-linked loci revealed the presence of one female-linked locus homologous with the B4GALNT1 gene, which is involved in the spermatogenesis pathway and hatchability. However, this locus was not sex-linked in the parental species, and the African catfish might also exhibit PSD.

The Snakeskin Gourami (Trichopodus pectoralis) Tends to Exhibit XX/XY Sex Determination

The snakeskin gourami (Trichopodus pectoralis) has a high meat yield and is one of the top five aquaculture freshwater fishes in Thailand. The species is not externally sexually dimorphic, and its sex determination system is unknown. Understanding the sex determination system of this species will contribute to its full-scale commercialization. In this study, a cytogenetic analysis did not reveal any between-sex differences in chromosomal patterns. However, we used genotyping-by-sequencing to identify 4 male-linked loci and 1 female-linked locus, indicating that the snakeskin gourami tends to exhibit an XX/XY sex determination system. However, we did not find any male-specific loci after filtering the loci for a ratio of 100:0 ratio of males:females. This suggests that the putative Y chromosome is young and that the sex determination region is cryptic. This approach provides solid information that can help identify the sex determination mechanism and potential sex determination regions in the snakeskin gourami, allowing further investigation of genetic improvements in the species.

The effects of weather variability on patterns of genetic diversity in Tasmanian bettongs

While the effects of climate (long-term, prevailing weather) on species abundance, range and genetic diversity have been widely studied, short-term, localized variations in atmospheric conditions (i.e., weather) can also rapidly alter species' geographical ranges and population sizes, but little is known about how they affect genetic diversity. We investigated the relationship between weather and range-wide genetic diversity in a marsupial, Bettongia gaimardi, using dynamic species distribution models (SDMs). Genetic diversity was lower in parts of the range where the weather-based SDM predicted high variability in probability of B. gaimardi occurrence during 1950-2009. This is probably an effect of lower population sizes and extinction-recolonization cycles in places with highly variable weather. Spatial variation in genetic diversity was also better predicted by mean probabilities of B. gaimardi occurrence from weather- than climate-based SDMs. Our results illustrate the importance of weather in driving population dynamics and species distributions on decadal timescales and thereby in affecting genetic diversity. Modelling the links between changing weather patterns, species distributions and genetic diversity will allow researchers to better forecast biological impacts of climate change.

An Investigation of ZZ/ZW and XX/XY Sex Determination Systems in North African Catfish (Clarias gariepinus, )

An investigation of sex-specific loci may provide important insights into fish sex determination strategies. This may be useful for biotechnological purposes, for example, to produce all-male or all-female fish for commercial breeding. The North African catfish species, Clarias gariepinus, has been widely adopted for aquaculture because its superior growth and disease resistance render the species suitable for hybridization with other catfish to improve the productivity and quality of fish meat. This species has either a ZZ/ZW or XX/XY sex determination system. Here, we investigate and characterize these systems using high-throughput genome complexity reduction sequencing as Diversity Arrays Technology. This approach was effective in identifying moderately sex-linked loci with both single-nucleotide polymorphisms (SNPs) and restriction fragment presence/absence (PA) markers in 30 perfectly sexed individuals of C. gariepinus. However, SNPs based markers were not found in this study. In total, 41 loci met the criteria for being moderately male-linked (with male vs. female ratios 80:20 and 70:30), while 25 loci were found to be moderately linked to female sex. No strictly male- or female-linked loci were detected. Seven moderately male-linked loci were partially homologous to some classes of transposable elements and three moderately male-linked loci were partially homologous to functional genes. Our data showed that the male heterogametic XX/XY sex determination system should co-exist with the ZZ/ZW system in C. gariepinus. Our finding of the co-existence of XX/XY and ZZ/ZW systems can be applied to benefit commercial breeding of this species in Thailand. This approach using moderately sex-linked loci provides a solid baseline for revealing sex determination mechanisms and identify potential sex determination regions in catfish, allowing further investigation of genetic improvements in breeding programs.

Genome-Wide Association Studies Detect Multiple QTLs for Productivity in Mesoamerican Diversity Panel of Common Bean Under Drought Stress

Drought stress is an important abiotic factor limiting common bean yield, with great impact on the production worldwide. Understanding the genetic basis regulating beans' yield and seed weight (SW) is a fundamental prerequisite for the development of superior cultivars. The main objectives of this work were to conduct genome-wide marker discovery by genotyping a Mesoamerican panel of common bean germplasm, containing cultivated and landrace accessions of broad origin, followed by the identification of genomic regions associated with productivity under two water regimes using different genome-wide association study (GWAS) approaches. A total of 11,870 markers were genotyped for the 339 genotypes, of which 3,213 were SilicoDArT and 8,657 SNPs derived from DArT and CaptureSeq. The estimated linkage disequilibrium extension, corrected for structure and relatedness (r ² _sv ), was 98.63 and 124.18 kb for landraces and breeding lines, respectively. Germplasm was structured into landraces and lines/cultivars. We carried out GWASs for 100-SW and yield in field environments with and without water stress for 3 consecutive years, using single-, segment-, and gene-based models. Higher number of associations at high stringency was identified for the SW trait under irrigation, totaling ∼185 QTLs for both single- and segment-based, whereas gene-based GWASs showed ∼220 genomic regions containing ∼650 genes. For SW under drought, 18 QTLs were identified for single- and segment-based and 35 genes by gene-based GWASs. For yield, under irrigation, 25 associations were identified, whereas under drought the total was 10 using both approaches. In addition to the consistent associations detected across experiments, these GWAS approaches provided important complementary QTL information (∼221 QTLs; 650 genes; r ² from 0.01% to 32%). Several QTLs were mined within or near candidate genes playing significant role in productivity, providing better understanding of the genetic mechanisms underlying these traits and making available molecular tools to be used in marker-assisted breeding. The findings also allowed the identification of genetic material (germplasm) with better yield performance under drought, promising to a common bean breeding program. Finally, the availability of this highly diverse Mesoamerican panel is of great scientific value for the analysis of any relevant traits in common bean.

Genome Complexity Reduction High-Throughput Genome Sequencing of Green Iguana (Iguana iguana) Reveal a Paradigm Shift in Understanding Sex-Chromosomal Linkages on Homomorphic X and Y Sex Chromosomes

The majority of lizards classified in the superfamily Iguanoidea have an XX/XY sex-determination system in which sex-chromosomal linkage shows homology with chicken (Gallus gallus) chromosome 15 (GGA15). However, the genomics of sex chromosomes remain largely unexplored owing to the presence of homomorphic sex chromosomes in majority of the species. Recent advances in high-throughput genome complexity reduction sequencing provide an effective approach to the identification of sex-specific loci with both single-nucleotide polymorphisms (SNPs) and restriction fragment presence/absence (PA), and a better understanding of sex chromosome dynamics in Iguanoidea. In this study, we applied Diversity Arrays Technology (DArTseq^TM) in 29 phenotypic sex assignments (14 males and 15 females) of green iguana (Iguana iguana). We confirmed a male heterogametic (XX/XY) sex determination mode in this species, identifying 29 perfectly sex-linked SNP/PA loci and 164 moderately sex-linked SNP/PA loci, providing evidence probably indicative of XY recombination. Three loci from among the perfectly sex-linked SNP/PA loci showed partial homology with several amniote sex chromosomal linkages. The results support the hypothesis of an ancestral super-sex chromosome with overlaps of partial sex-chromosomal linkages. However, only one locus among the moderately sex-linked loci showed homology with GGA15, which suggests that the specific region homologous to GGA15 was located outside the non-recombination region but in close proximity to this region of the sex chromosome in green iguana. Therefore, the location of GGA15 might be further from the putative sex-determination locus in green iguana. This is a paradigm shift in understanding linkages on homomorphic X and Y sex chromosomes. The DArTseq platform provides an easy-to-use strategy for future research on the evolution of sex chromosomes in Iguanoidea, particularly for non-model species with homomorphic or highly cryptic sex chromosomes.

Recurrent breakdown and rebalance of segregation distortion in the genomes: battle for the transmission advantage

Mendel's laws state that each of the two alleles would segregate during gamete formation and show the same transmission ratio in the next generation. However, an unexpected biased allele transmission was first detected in Drosophila a century ago, and was subsequently observed in other animals, plants, and microorganisms. Such segregation distortion (SD) shows substantial effects in population structure and fitness of the progenies, which would ultimately lead to reproductive isolation and speciation. Here, we trace the early investigations on the violation of Mendelian genetic principle, which appears as a wide-existence phenomenon rather than a case of exception. The occurence of SD in the whole genome was observed in a number of plant species at the single- and multi-locus level. Biased transmission ratio might occur at meiosis stage due to asymmetric movement of the chromosome; transmission ratio advantage is also caused by interaction and battle between the alleles from respective genomes at the genetic and molecular level. The origin of a SD system is likely to be determined by coevolution of the killer and protector via recurrent breakdown or rebalance loop. These updated understandings also promote genetic improvement of hybrid crops.

Genome-wide SNP analysis of Siamese cobra (Naja kaouthia) reveals the molecular basis of transitions between Z and W sex chromosomes and supports the presence of an ancestral super-sex chromosome in amniotes

Elucidation of the process of sex chromosome differentiation is necessary to understand the dynamics of evolutionary mechanisms in organisms. The W sex chromosome of the Siamese cobra (Naja kaouthia) contains a large number of repeats and shares amniote sex chromosomal linkages. Diversity Arrays Technology provides an effective approach to identify sex-specific loci that are epoch-making, to understand the dynamics of molecular transitions between the Z and W sex chromosomes in a snake lineage. From a total of 543 sex-specific loci, 90 showed partial homology with sex chromosomes of several amniotes and 89 loci were homologous to transposable elements. Two loci were confirmed as W-specific nucleotides after PCR amplification. These loci might result from a sex chromosome differentiation process and involve putative sex-determination regions in the Siamese cobra. Sex-specific loci shared linkage homologies among amniote sex chromosomes, supporting an ancestral super-sex chromosome.

Genome-wide genetic diversity and population structure of Garcinia kola (Heckel) in Benin using DArT-Seq technology

Garcinia kola (Heckel) is a versatile tree indigenous to West and Central Africa. All parts of the tree have value in traditional medicine. Natural populations of the species have declined over the years due to overexploitation. Assessment of genetic diversity and population structure of G. kola is important for its management and conservation. The present study investigates the genetic diversity and population structure of G. kola populations in Benin using ultra-high-throughput diversity array technology (DArT) single nucleotide polymorphism (SNP) markers. From the 102 accessions sampled, two were excluded from the final dataset owing to poor genotyping coverage. A total of 43,736 SNPs were reported, of which 12,585 were used for analyses after screening with quality control parameters including Minor allele frequency (≥ 0.05), call rate (≥ 80%), reproducibility (≥ 95%), and polymorphic information content (≥ 1%). Analysis revealed low genetic diversity with expected heterozygosity per population ranging from 0.196 to 0.228. Pairwise F-statistics (FST) revealed low levels of genetic differentiation between populations while an Analysis of molecular variance (AMOVA) indicated that the majority of variation (97.86%) was within populations. Population structure analysis through clustering and discriminant analysis on principal component revealed two admixed clusters, implying little genetic structure. However, the model-based maximum likelihood in Admixture indicated only one genetic cluster. The present study indicated low genetic diversity of G. kola, and interventions are needed to be tailored towards its conservation.

Evaluation of the Minimum Sampling Design for Population Genomic and Microsatellite Studies: An Analysis Based on Wild Maize

Massive parallel sequencing (MPS) is revolutionizing the field of molecular ecology by allowing us to understand better the evolutionary history of populations and species, and to detect genomic regions that could be under selection. However, the economic and computational resources needed generate a tradeoff between the amount of loci that can be obtained and the number of populations or individuals that can be sequenced. In this work, we analyzed and compared two simulated genomic datasets fitting a hierarchical structure, two extensive empirical genomic datasets, and a dataset comprising microsatellite information. For all datasets, we generated different subsampling designs by changing the number of loci, individuals, populations, and individuals per population to test for deviations in classic population genetics parameters (H_S, F_IS, F_ST). For the empirical datasets we also analyzed the effect of sampling design on landscape genetic tests (isolation by distance and environment, central abundance hypothesis). We also tested the effect of sampling a different number of populations in the detection of outlier SNPs. We found that the microsatellite dataset is very sensitive to the number of individuals sampled when obtaining summary statistics. F_IS was particularly sensitive to a low sampling of individuals in the simulated, genomic, and microsatellite datasets. For the empirical and simulated genomic datasets, we found that as long as many populations are sampled, few individuals and loci are needed. For the empirical datasets, we found that increasing the number of populations sampled was important in obtaining precise landscape genetic estimates. Finally, we corroborated that outlier tests are sensitive to the number of populations sampled. We conclude by proposing different sampling designs depending on the objectives.

Discovery of thermophilic Bacillales using reduced-representation genotyping for identification

Background

This study demonstrates the use of reduced-representation genotyping to provide preliminary identifications for thermophilic bacterial isolates. The approach combines restriction enzyme digestion and PCR with next-generation sequencing to provide thousands of short-read sequences from across the bacterial genomes. Isolates were obtained from compost, hot water systems, and artesian bores of the Great Artesian Basin. Genomic DNA was double-digested with two combinations of restriction enzymes followed by PCR amplification, using a commercial provider of DArTseq™, Diversity Arrays Technology Pty Ltd. (Canberra, Australia). The resulting fragments which formed a reduced-representation of approximately 2.3% of the genome were sequenced. The sequence tags obtained were aligned against all available RefSeq bacterial genome assemblies by BLASTn to identify the nearest reference genome.

Results

Based on the preliminary identifications, a total of 99 bacterial isolates were identified to species level, from which 8 isolates were selected for whole-genome sequencing to assess the identification results. Novel species and strains were discovered within this set of isolates. The preliminary identifications obtained by reduced-representation genotyping, as well as identifications obtained by BLASTn alignment of the 16S rRNA gene sequence, were compared with those derived from the whole-genome sequence data, using the same RefSeq sequence database for the three methods. Identifications obtained with reduced-representation sequencing agreed with the identifications provided by whole-genome sequencing in 100% of cases. The identifications produced by BLASTn alignment of 16S rRNA gene sequence to the same database differed from those provided by whole-genome sequencing in 37.5% of cases, and produced ambiguous identifications in 50% of cases.

Conclusions

Previously, this method has been successfully demonstrated for use in bacterial identification for medical microbiology. This study demonstrates the first successful use of DArTseq™ for preliminary identification of thermophilic bacterial isolates, providing results in complete agreement with those obtained from whole-genome sequencing of the same isolates. The growing database of bacterial genome sequences provides an excellent resource for alignment of reduced-representation sequence data for identification purposes, and as the available sequenced genomes continue to grow, the technique will become more effective.

Whole-genome diversity, population structure and linkage disequilibrium analysis of globally diverse wheat genotypes using genotyping-by-sequencing DArTseq platform

In this study, 129 wheat genotypes from globally diverse origins were genotyped using DArTseq (SilicoDArT and SNP) markers. After filtering markers for quality-filtering, 14,270 SilicoDArTs and 6484 SNPs were retained and used for genetic diversity, population structure and linkage disequilibrium analyses. The highest number of SilicoDArT and SNP markers mapped on genome A and B compared to genome D. In both marker types, polymorphism information content (PIC) values ranged from 0.1 to 0.5, while > 0.80% of SilicoDArTs and > 0.44% SNPs showed PIC value more than median (0.25%). Un-weighted Neighbor Joining cluster analysis and Bayesian-based model population structure grouped wheat genotypes into three and four clusters, respectively. Principal component analysis and discriminant analysis of principal component results showed highly match with cluster and population structure analysis. Linkage disequilibrium (LD) was more extensive in both marker types, while graphical display of LD decay for both marker types showed that LD declined in the region close to 15 kbp, where r 2-values corresponded to r 2 = 0.16. Overall, our genetic diversity analysis showed high level of variation in studied wheat genotypes, even though there was no relationship between wheat grouping and origins. This might be attributed to admixture level that occurred during long-term natural selection of wheat genotypes in different parts of the world. Highly diverse wheat genotypes used in this study may possess unique genes and are useful sources in breeding programs to improve grain yield and quality.

Genetic Mapping and Discovery of the Candidate Gene for Black Seed Coat Color in Watermelon (Citrullus lanatus)

Seed coat color is an important trait highly affecting the seed quality and flesh appearance of watermelon (Citrullus lanatus). However, the molecular regulation mechanism of seed coat color in watermelon is still unclear. In the present study, genetic analysis was performed by evaluating F₁, F₂ and BC₁ populations derived from two parental lines (9904 with light yellow seeds and Handel with black seeds), suggesting that a single dominant gene controls the black seed coat. The initial mapping result revealed a region of interest spanning 370 kb on chromosome 3. Genetic mapping with CAPS and SNP markers narrowed down the candidate region to 70.2 kb. Sequence alignment of the three putative genes in the candidate region suggested that there was a single-nucleotide insertion in the coding region of Cla019481 in 9904, resulting in a frameshift mutation and premature stop codon. The results indicated that Cla019481 named ClCS1 was the candidate gene for black seed coat color in watermelon. In addition, gene annotation revealed that Cla019481 encoded a polyphenol oxidase (PPO), which involved in the oxidation step of the melanin biosynthesis. This research finding will facilitate maker-assisted selection in watermelon and provide evidence for the study of black seed coat coloration in plants.

Identification and Molecular Mapping of a Gummy Stem Blight Resistance Gene in Wild Watermelon (Citrullus amarus) Germplasm PI 189225

Gummy stem blight (GSB), caused by Stagonosporopsis cucurbitacearum (syn. Didymella bryoniae), is a destructive foliar disease of watermelon in areas with hot and humid climates. The wild watermelon germplasm PI 189225 is a known source of resistance to GSB. The identification and use of molecular markers linked to resistance genes in the wild-type germplasm will speed up the introgression of GSB resistance into new watermelon varieties. An F₂ segregating population was obtained from a cross between the resistant wild watermelon genotype PI 189225 and the susceptible genotype K3. The F₂-derived F₃ families were inoculated with a single isolate of S. cucurbitacearum (JS002) from Jiangsu Academy of Agricultural Sciences. The results of the genetic analysis demonstrated that GSB resistance in PI 189225 was controlled by a major quantitative trait locus (QTL), temporarily designated Qgsb8.1. Based on the results of bulk sergeant analysis and sequencing, one associated region spanning 5.7 Mb (10,358,659 to 16,101,517) on chromosome 8 was identified as responsible for the resistance to GSB using the Δ(single-nucleotide polymorphism [SNP]-index) method. The result of a QTL linkage analysis with Kompetitive allele-specific PCR (KASP) SNP markers further mapped the GSB resistance locus between the SNP markers KASP_JS9383 and KASP_JS9168 in a region of 571.27 kb on chromosome 8. According to the watermelon gene annotation database, the region contains approximately 19 annotated genes and, of these 19 genes, 2 are disease resistance gene analogs: Cla001017 (coiled-coil nucleotide-binding site leucine-rich repeat resistance protein) and Cla001019 (pathogenesis related). Reverse-transcription quantitative PCR demonstrated that the expression of the two genes changed following S. cucurbitacearum infection, suggesting that they play important roles in GSB resistance in watermelon. This result will facilitate fine mapping and cloning of the Qgsb8.1 locus, and the linked markers will further provide a useful tool for marker-assisted selection of this locus in watermelon breeding programs.

Multiple Small-Effect Alleles of Indica Origin Enhance High Iron-Associated Stress Tolerance in Rice Under Field Conditions in West Africa

Understanding the genetics of field-based tolerance to high iron-associated (HIA) stress in rice can accelerate the development of new varieties with enhanced yield performance in West African lowland ecosystems. To date, few field-based studies have been undertaken to rigorously evaluate rice yield performance under HIA stress conditions. In this study, two NERICA × O. sativa bi-parental rice populations and one O.sativa diversity panel consisting of 296 rice accessions were evaluated for grain yield and leaf bronzing symptoms over multiple years in four West African HIA stress and control sites. Mapping of these traits identified a large number of QTLs and single nucleotide polymorphisms (SNPs) associated with stress tolerance in the field. Favorable alleles associated with tolerance to high levels of iron in anaerobic rice soils were rare and almost exclusively derived from the indica subpopulation, including the most favorable alleles identified in NERICA varieties. These findings highlight the complex genetic architecture underlying rice response to HIA stress and suggest that a recurrent selection program focusing on an expanded indica genepool could be productively used in combination with genomic selection to increase the efficiency of selection in breeding programs designed to enhance tolerance to this prevalent abiotic stress in West Africa.

Identifying sex-linked markers in Litoria aurea: a novel approach to understanding sex chromosome evolution in an amphibian

Few taxa exhibit the variability of sex-determining modes as amphibians. However, due to the presence of homomorphic sex chromosomes in many species, this phenomenon has been difficult to study. The Australian frog, Litoria aurea, has been relatively well studied over the past 20 years due to widespread declines largely attributable to chytrid fungus. However, it has been subject to few molecular studies and its mode of sex determination remained unknown. We applied DArTseq™ to develop sex-linked single nucleotide polymorphisms (SNPs) and restriction fragment presence/absence (PA) markers in 44 phenotypically sexed L. aurea individuals from the Molonglo River in NSW, Australia. We conclusively identified a male heterogametic (XX-XY) sex determination mode in this species, identifying 11 perfectly sex-linked SNP and six strongly sex-linked PA markers. We identified a further 47 moderately sex-linked SNP loci, likely serving as evidence indicative of XY recombination. Furthermore, within these 47 loci, a group of nine males were found to have a feminised Y chromosome that significantly differed to all other males. We postulate ancestral sex-reversal as a means for the evolution of this now pseudoautosomal region on the Y chromosome. Our findings present new evidence for the ‘fountain of youth’ hypothesis for the retention of homomorphic sex chromosomes in amphibians and describe a novel approach for the study of sex chromosome evolution in amphibia.

Whole genome resequencing of watermelons to identify single nucleotide polymorphisms related to flesh color and lycopene content

Cultivated watermelon (Citrullus lanatus) is one of the most important food crops in the Cucurbitaceae family. Diversification after domestication has led cultivated watermelons to exhibit diverse fruit flesh colors, including red, yellow, and orange. Recently, there has been increased interest in red-fleshed watermelons because they contain the antioxidant cis-isomeric lycopene. We performed whole genome resequencing (WGRS) of 24 watermelons with different flesh colors to identify single-nucleotide polymorphisms (SNPs) related to high lycopene content. The resequencing data revealed 203,894-279,412 SNPs from read mapping between inbred lines and the 97103 reference genome. In total, 295,065 filtered SNPs were identified, which had an average polymorphism information content of 0.297. Most of these SNPs were intergenic (90.1%) and possessed a transversion (Tv) rate of 31.64%. Overall, 2,369 SNPs were chosen at 0.5 Mb physical intervals to analyze genetic diversity across the 24 inbred lines. A neighbor-joining dendrogram and principal coordinate analysis (PCA) based on the 2,369 SNPs revealed that the 24 inbred lines could be grouped into high and low lycopene-type watermelons. In addition, we analyzed SNPs that could discriminate high lycopene content, red-fleshed watermelon from low lycopene, yellow or orange watermelon inbred lines. For validation, 19 SNPs (designated as WMHL1-19) were chosen randomly, and cleavage amplified polymorphic sequence (CAPS) markers were designed. Genotyping of the above 24 lines and 12 additional commercial cultivars using WMHL1-19 CAPS markers resulted in match rates of over 0.92 for most validated markers in correlation with the flesh color phenotypes. Our results provide valuable genomic information regarding the high lycopene content phenotype of red-fleshed cultivated watermelons, and the identified SNPs will be useful for the development of molecular markers in the marker-assisted breeding of watermelons with high lycopene content.

Genetic mapping and development of molecular markers for a candidate gene locus controlling rind color in watermelon

ClCG08G017810 (ClCGMenG) encoding a 2-phytyl-1,4-beta-naphthoquinone methyltransferase protein is associated with formation of dark green versus light green rind color in watermelon. Rind color is an important agronomic trait in watermelon [Citrullus lanatus (Thunb.) Matsum. and Nakai], but the underlying molecular mechanism for this trait is not fully known. In the present study, we identified a single locus on chromosome 8 accounting for watermelon rind color (dark green vs. light green). Genetic analysis of F₁, F₂, and BC₁ populations derived from two parental lines (9904 with dark green rind and Handel with light green rind) revealed that the watermelon rind color (dark green vs. light green) is controlled by a single locus, and dark green is dominant to light green rind. Initial mapping revealed a region of interest spanning 2.07 Mb on chromosome 8. Genetic mapping with CAPS and SNP markers narrowed down the candidate region to 31.4 kb. Gene annotation of the corresponding region in the reference genome revealed the ClCG08G017810 gene sequence encoding the 2-phytyl-1,4-beta-naphthoquinone methyltransferase protein. The sequence alignment of the candidate gene with the two parental lines suggested a nonsynonymous SNP mutation in the coding region of ClCG08G017810, converting an arginine (R) to glycine (G). The SNP might be associated with rind color of 103 watermelon germplasm lines investigated in this study. The qRT-PCR analysis revealed higher expression of ClCG08G017810 in dark green rind than in light green rind. Therefore, ClCG08G017810 is a candidate gene associated with watermelon rind color. The present study facilitates marker-assisted selection useful for the development of cultivars with desirable rind color.

Whole Genome Diversity, Population Structure, and Linkage Disequilibrium Analysis of Chickpea (Cicer arietinum L.) Genotypes Using Genome-Wide DArTseq-Based SNP Markers

Characterization of genetic diversity, population structure, and linkage disequilibrium is a prerequisite for proper management of breeding programs and conservation of genetic resources. In this study, 186 chickpea genotypes, including advanced "Kabuli" breeding lines and Iranian landrace "Desi" chickpea genotypes, were genotyped using DArTseq-Based single nucleotide polymorphism (SNP) markers. Out of 3339 SNPs, 1152 markers with known chromosomal position were selected for genome diversity analysis. The number of mapped SNP markers varied from 52 (LG8) to 378 (LG4), with an average of 144 SNPs per linkage group. The chromosome size that was covered by SNPs varied from 16,236.36 kbp (LG8) to 67,923.99 kbp (LG5), while LG4 showed a higher number of SNPs, with an average of 6.56 SNPs per Mbp. Polymorphism information content (PIC) value of SNP markers ranged from 0.05 to 0.50, with an average of 0.32, while the markers on LG4, LG6, and LG8 showed higher mean PIC value than average. Unweighted neighbor joining cluster analysis and Bayesian-based model population structure grouped chickpea genotypes into four distinct clusters. Principal component analysis (PCoA) and discriminant analysis of principal component (DAPC) results were consistent with that of the cluster and population structure analysis. Linkage disequilibrium (LD) was extensive and LD decay in chickpea germplasm was relatively low. A few markers showed r2 ≥ 0.8, while 2961 pairs of markers showed complete LD (r2 = 1), and a huge LD block was observed on LG4. High genetic diversity and low kinship value between pairs of genotypes suggest the presence of a high genetic diversity among the studied chickpea genotypes. This study also demonstrates the efficiency of DArTseq-based SNP genotyping for large-scale genome analysis in chickpea. The genotypic markers provided in this study are useful for various association mapping studies when combined with phenotypic data of different traits, such as seed yield, abiotic, and biotic stresses, and therefore can be efficiently used in breeding programs to improve chickpea.

Indications of strong adaptive population genetic structure in albacore tuna (Thunnus alalunga) in the southwest and central Pacific Ocean

Albacore tuna (Thunnus alalunga) has a distinctly complex life history in which juveniles and adults separate geographically but at times inhabit the same spaces sequentially. The species also migrates long distances and presumably experiences varied regimes of physical stress over a lifetime. There are, therefore, many opportunities for population structure to arise based on stochastic differences or environmental factors that promote local adaptation. However, with the extent of mobility consistently demonstrated by tagged individuals, there is also a strong argument for panmixia within an ocean basin. It is important to confirm such assumptions from a population genetics standpoint for this species in particular because albacore is one of the principal market tuna species that sustains massive global fisheries and yet is also a slow-growing temperate tuna. Consequently, we used 1,837 neutral SNP loci and 89 loci under potential selection to analyze population genetic structure among five sample groups collected from the western and central South Pacific. We found no evidence to challenge panmixia at neutral loci, but strong indications of structuring at adaptive loci. One population sample, from French Polynesia in 2004, was particularly differentiated. Unfortunately, the current study cannot infer whether the divergence is geographic or temporal, or possibly caused by sample distribution. We encourage future studies to include potentially adaptive loci and to continue fine-scale observations within an ocean basin, and not to assume genome-wide panmixia.

Short-read fastA files dataset from complexity-reduced genotyping by sequencing data of bacterial isolates from a public hospital in Australia

This data article contains short-read sequences (length 30-69 bp) obtained from complexity-reduced genotyping by sequencing (GBS) of 165 samples bacterial isolates from hospital patients in the Australian Capital Territory, between 2013 and 2015. These samples represented 14 bacterial species. Data format is shown as filtered fastA files obtained from an Illumina HiSeq2500 sequencer. The experimental factors of this research used three complexity reduction methods with three combinations of restriction enzymes: PstI with MseI, PstI with HpaII and MseI with HpaII.

Divergence with gene flow is driven by local adaptation to temperature and soil phosphorus concentration in teosinte subspecies (Zea mays parviglumis and Zea mays mexicana)

Patterns of genomic divergence between hybridizing taxa can be heterogeneous along the genome. Both differential introgression and local adaptation may contribute to this pattern. Here, we analysed two teosinte subspecies, Zea mays ssp. parviglumis and ssp. mexicana, to test whether their divergence has occurred in the face of gene flow and to infer which environmental variables have been important drivers of their ecological differentiation. We generated 9,780 DArTseqTM SNPs for 47 populations, and used an additional data set containing 33,454 MaizeSNP50 SNPs for 49 populations. With these data, we inferred features of demographic history and performed genome wide scans to determine the number of outlier SNPs associated with climate and soil variables. The two data sets indicate that divergence has occurred or been maintained despite continuous gene flow and/or secondary contact. Most of the significant SNP associations were to temperature and to phosphorus concentration in the soil. A large proportion of these candidate SNPs were located in regions of high differentiation that had been identified previously as putative inversions. We therefore propose that genomic differentiation in teosintes has occurred by a process of adaptive divergence, with putative inversions contributing to reduced gene flow between locally adapted populations.

Identification of bacterial isolates from a public hospital in Australia using complexity-reduced genotyping

Bacterial identification methods used in routine identification of pathogens in medical microbiology include a combination approach of biochemical tests, mass spectrometry or molecular biology techniques. Extensive publicly-available databases of DNA sequence data from pathogenic bacteria have been amassed in recent years; this provides an opportunity for using bacterial genome sequencing for identification purposes. Whole genome sequencing is increasing in popularity, although at present it remains a relatively expensive approach to bacterial identification and typing. Complexity-reduced bacterial genome sequencing provides an alternative. We evaluate genomic complexity-reduction using restriction enzymes and sequencing to identify bacterial isolates. A total of 165 bacterial isolates from hospital patients in the Australian Capital Territory, between 2013 and 2015 were used in this study. They were identified and typed by the Microbiology Department of Canberra Public Hospital, and represented 14 bacterial species. DNA extractions from these samples were processed using a combination of the restriction enzymes PstI with MseI, PstI with HpaII and MseI with HpaII. The resulting sequences (length 30-69 bp) were aligned against publicly available bacterial genome and plasmid sequences. Results of the alignment were processed using a bioinformatics pipeline developed for this project, Currito3.1 DNA Fragment Analysis Software. All 165 samples were correctly identified to genus and species by each of the three combinations of restriction enzymes. A further 35 samples typed to the level of strain identified and compared for consistency with MLST typing data and in silico MLST data derived from the nearest sequenced candidate reference. The high level of agreement between bacterial identification using complexity-reduced genome sequencing and standard hospital identifications indicating that this new approach is a viable alternative for identification of bacterial isolates derived from pathology specimens. The effectiveness of species identification and in particular, strain typing, depends on access to a comprehensive and taxonomically accurate bacterial genome sequence database containing relevant bacterial species and strains.

Construction of A High-Density Genetic Map and Mapping of Fruit Traits in Watermelon (Citrullus Lanatus L.) Based on Whole-Genome Resequencing

Watermelon (Citrullus lanatus L.) is an important horticultural crop that is grown worldwide and has a high economic value. To dissect the loci associated with important horticultural traits and to analyze the genetic and genomic information of this species, a high-density genetic map was constructed based on whole-genome resequencing (WGR), a powerful high-resolution method for single-nucleotide polymorphism (SNP) marker development, genetic map construction, and gene mapping. Resequencing of both parental lines and 126 recombinant inbred lines (RIL) resulted in the detection of 178,762 single-nucleotide polymorphism (SNP) markers in the parental lines at a sequencing depth greater than four-fold. Additionally, 2132 recombination bin markers comprising 103,029 SNP markers were mapped onto 11 linkage groups (LGs). Substantially more SNP markers were mapped to the genetic map compared with other recent studies. The total length of the linkage map was 1508.94 cM, with an average distance of 0.74 cM between adjacent bin markers. Based on this genetic map, one locus for fruit bitterness, one locus for rind color, and one locus for seed coat color with high LOD scores (58.361, 18.353, 26.852) were identified on chromosome 1, chromosome 8, and chromosome 3, respectively. These prominent loci were identified in a region of 6.16 Mb, 2.07 Mb, and 0.37 Mb, respectively. On the basis of current research, the high-density map and mapping results will provide a valuable tool for identifying candidate genes, map-based gene cloning, comparative mapping, and marker-assisted selection (MAS) in watermelon breeding.

Molecular Markers Associated to Two Non-allelic Genic Male Sterility Genes in Peppers (Capsicum annuum L.)

Male sterility is of high importance in hybrid seed production of hot and sweet peppers. Genic (or nuclear) male sterility (GMS) is a simply inherited (usually monogenic recessive) and highly stable trait. However, one major disadvantage of using GMS is 1:1 segregation of male sterile to male fertile plants in every subsequent generation. Molecular markers tightly linked to genic male sterility (ms) genes would facilitate an efficient and rapid transfer of ms genes into different genetic backgrounds through marker-assisted backcrossing. The two non-allelic genic male sterility genes ms3 and ms w in hot and sweet pepper backgrounds, respectively, are monogenic recessive. Genotyping by sequencing (GBS) in an F2 population segregating for ms3 gene in hot pepper and in an F6 inbred near-isogenic line (NIL) population segregating for ms w gene in sweet pepper yielded 9,713 and 7,453 single nucleotide polymorphism markers, respectively. Four candidate SNPs co-segregating with ms3 gene and one co-segregating with ms w gene were identified by bulk segregant analysis and physically mapped to chromosomes 1 and 5, respectively. In hot pepper, two markers [HPGMS2 (CAPS) and HPGMS3 (dCAPS)] located 3.83 cM away from the ms3 gene and in sweet pepper the dCAPS marker SPGMS1 co-segregated (completely linked) with the ms w gene were developed. These markers will increase the efficacy of the male sterility genes for pepper breeding, as they can be useful in developing the genic male sterile lines in parental inbred lines of commercial hybrids through marker-assisted backcrossing, hybrid seed production, and genetic purity testing of hybrid seeds.

Conservation of Sex-Linked Markers among Conspecific Populations of a Viviparous Skink, Niveoscincus ocellatus, Exhibiting Genetic and Temperature-Dependent Sex Determination

Sex determination systems are exceptionally diverse and have undergone multiple and independent evolutionary transitions among species, particularly reptiles. However, the mechanisms underlying these transitions have not been established. Here, we tested for differences in sex-linked markers in the only known reptile that is polymorphic for sex determination system, the spotted snow skink, Niveoscincus ocellatus, to quantify the genomic differences that have accompanied this transition. In a highland population, sex is determined genetically, whereas in a lowland population, offspring sex ratio is influenced by temperature. We found a similar number of sex-linked loci in each population, including shared loci, with genotypes consistent with male heterogamety (XY). However, population-specific linkage disequilibrium suggests greater differentiation of sex chromosomes in the highland population. Our results suggest that transitions between sex determination systems can be facilitated by subtle genetic differences.

Effective BAC clone anchoring with genotyping-by-sequencing and Diversity Arrays Technology in a large genome cereal rye

Identification of bacterial artificial chromosome (BAC) clones containing specific sequences is a prerequisite for many applications, such as physical map anchoring or gene cloning. Existing BAC library screening strategies are either low-throughput or require a considerable initial input of resources for platform establishment. We describe a high-throughput, reliable, and cost-effective BAC library screening approach deploying genotyping platforms which are independent from the availability of sequence information: a genotyping-by-sequencing (GBS) method DArTSeq and the microarray-based Diversity Arrays Technology (DArT). The performance of these methods was tested in a very large and complex rye genome. The DArTseq approach delivered superior results: a several fold higher efficiency of addressing genetic markers to BAC clones and anchoring of BAC clones to genetic map and also a higher reliability. Considering the sequence independence of the platform, the DArTseq-based library screening can be proposed as an attractive method to speed up genomics research in resource poor species.

Populations of doubled haploids for genetic mapping in hexaploid winter triticale

To create a framework for genetic dissection of hexaploid triticale, six populations of doubled haploid (DH) lines were developed from pairwise hybrids of high-yielding winter triticale cultivars. The six populations comprise between 97 and 231 genotyped DH lines each, totaling 957 DH lines. A consensus genetic map spans 4593.9 cM is composed of 1576 unique DArT markers. The maps reveal several structural rearrangements in triticale genomes. In preliminary tests of the populations and maps, markers specific to wheat segments of the engineered rye chromosome 1R (RM1B) were identified. Example QTL mapping of days to heading in cv. Krakowiak revealed loci on chromosomes 2BL and 2R responsible for extended vernalization requirement, and candidate genes were identified. The material is available to all parties interested in triticale genetics.

Genome-wide comparisons reveal evidence for a species complex in the black-lip pearl oyster Pinctada margaritifera (Bivalvia: Pteriidae)

Evolutionary relationships in the black-lip pearl oyster Pinctada margaritifera which is highly valued for pearl production remain poorly understood. This species possesses an 18,000 km Indo-Pacific natural distribution, and its current description includes six subspecies defined exclusively on morphological characters. To evaluate its taxonomic identity using molecular data, 14 populations in both the Indian and Pacific Oceans (n = 69), and the congeneric taxa P. maxima and P. mazatlanica (n = 29 and n = 10, respectively) were sampled. Phylogenomic reconstruction was carried out using both 8,308 genome-wide SNPs and 10,000 dominant loci (DArTseq PAVs). Reconstructions using neighbour-joining (Nei's 1972 distances), maximum likelihood and Bayesian approaches all indicate that the taxonomy of P. margaritifera is quite complex, with distinct evolutionary significant units (ESUs) identified within Tanzanian and Iranian populations. Contrastingly, phylogenies generated for Pacific Ocean oysters resolved a large monophyletic clade, suggesting little support for two current morphological subspecies classifications. Furthermore, P. mazatlanica formed a basal clade closest to French Polynesian P. margaritifera, suggesting it may be conspecific. Collectively, these findings provide evidence that P. margaritifera comprises a species complex, perhaps as a result of population fragmentation and increased divergence at range limits.

A High-Density Integrated DArTseq SNP-Based Genetic Map of Pisum fulvum and Identification of QTLs Controlling Rust Resistance

Pisum fulvum, a wild relative of pea is an important source of allelic diversity to improve the genetic resistance of cultivated species against fungal diseases of economic importance like the pea rust caused by Uromyces pisi. To unravel the genetic control underlying resistance to this fungal disease, a recombinant inbred line (RIL) population was generated from a cross between two P. fulvum accessions, IFPI3260 and IFPI3251, and genotyped using Diversity Arrays Technology. A total of 9,569 high-quality DArT-Seq and 8,514 SNPs markers were generated. Finally, a total of 12,058 markers were assembled into seven linkage groups, equivalent to the number of haploid chromosomes of P. fulvum and P. sativum. The newly constructed integrated genetic linkage map of P. fulvum covered an accumulated distance of 1,877.45 cM, an average density of 1.19 markers cM^-1 and an average distance between adjacent markers of 1.85 cM. The composite interval mapping revealed three QTLs distributed over two linkage groups that were associated with the percentage of rust disease severity (DS%). QTLs UpDSII and UpDSIV were located in the LGs II and IV respectively and were consistently identified both in adult plants over 3 years at the field (Córdoba, Spain) and in seedling plants under controlled conditions. Whenever they were detected, their contribution to the total phenotypic variance varied between 19.8 and 29.2. A third QTL (UpDSIV.2) was also located in the LGIVand was environmentally specific as was only detected for DS % in seedlings under controlled conditions. It accounted more than 14% of the phenotypic variation studied. Taking together the data obtained in the study, it could be concluded that the expression of resistance to fungal diseases in P. fulvum originates from the resistant parent IFPI3260.

Assessment of Genetic Variation and Population Structure of Diverse Rice Genotypes Adapted to Lowland and Upland Ecologies in Africa Using SNPs

Using interspecific crosses involving Oryza glaberrima Steud. as donor and O. sativa L. as recurrent parents, rice breeders at the Africa Rice Center developed several 'New Rice for Africa (NERICA)' improved varieties. A smaller number of interspecific and intraspecific varieties have also been released as 'Advanced Rice for Africa (ARICA)'. The objective of the present study was to investigate the genetic variation, relatedness, and population structure of 330 widely used rice genotypes in Africa using DArTseq-based single nucleotide polymorphisms (SNPs). A sample of 11 ARICAs, 85 NERICAs, 62 O. sativa spp. japonica, and 172 O. sativa spp. indica genotypes were genotyped with 27,560 SNPs using diversity array technology (DArT)-based sequencing (DArTseq) platform. Nearly 66% of the SNPs were polymorphic, of which 15,020 SNPs were mapped to the 12 rice chromosomes. Genetic distance between pairs of genotypes that belong to indica, japonica, ARICA, and NERICA varied from 0.016 to 0.623, from 0.020 to 0.692, from 0.075 to 0.763, and from 0.014 to 0.644, respectively. The proportion of pairs of genotypes with genetic distance > 0.400 was the largest within NERICAs (35.1% of the pairs) followed by ARICAs (18.2%), japonica (17.4%), and indica (5.6%). We found one pair of japonica, 11 pairs of indica, and 35 pairs of NERICA genotypes differing by <2% of the total scored alleles, which was due to 26 pairs of genotypes with identical pedigrees. Cluster analysis, principal component analysis, and the model-based population structure analysis all revealed two distinct groups corresponding to the lowland (primarily indica and lowland NERICAs) and upland (japonica and upland NERICAs) growing ecologies. Most of the interspecific lowland NERICAs formed a sub-group, likely caused by differences in the O. glaberrima genome as compared with the indica genotypes. Analysis of molecular variance revealed very great genetic differentiation (F_ST = 0.688) between the lowland and upland ecologies, and 31.2% of variation attributable to differences within cluster groups. About 8% (1,197 of 15,020) of the 15,020 SNPs were significantly (P < 0.05) different between the lowland and upland ecologies and formed contrasting haplotypes that could clearly discriminate lowland from upland genotypes. This is the first study using high density markers that characterized NERICA and ARICA varieties in comparison with indica and japonica varieties widely used in Africa, which could aid rice breeders on parent selection for developing new improved rice germplasm.

High density SNP and DArT-based genetic linkage maps of two closely related oil palm populations

Oil palm (Elaeis guineensis Jacq.) is an outbreeding perennial tree crop with long breeding cycles, typically 12 years. Molecular marker technologies can greatly improve the breeding efficiency of oil palm. This study reports the first use of the DArTseq platform to genotype two closely related self-pollinated oil palm populations, namely AA0768 and AA0769 with 48 and 58 progeny respectively. Genetic maps were constructed using the DArT and SNP markers generated in combination with anchor SSR markers. Both maps consisted of 16 major independent linkage groups (2n = 2× = 32) with 1399 and 1466 mapped markers for the AA0768 and AA0769 populations, respectively, including the morphological trait "shell-thickness" (Sh). The map lengths were 1873.7 and 1720.6 cM with an average marker density of 1.34 and 1.17 cM, respectively. The integrated map was 1803.1 cM long with 2066 mapped markers and average marker density of 0.87 cM. A total of 82% of the DArTseq marker sequence tags identified a single site in the published genome sequence, suggesting preferential targeting of gene-rich regions by DArTseq markers. Map integration of higher density focused around the Sh region identified closely linked markers to the Sh, with D.15322 marker 0.24 cM away from the morphological trait and 5071 bp from the transcriptional start of the published SHELL gene. Identification of the Sh marker demonstrates the robustness of using the DArTseq platform to generate high density genetic maps of oil palm with good genome coverage. Both genetic maps and integrated maps will be useful for quantitative trait loci analysis of important yield traits as well as potentially assisting the anchoring of genetic maps to genomic sequences.

In-depth genome characterization of a Brazilian common bean core collection using DArTseq high-density SNP genotyping

BACKGROUND

Common bean is a legume of social and nutritional importance as a food crop, cultivated worldwide especially in developing countries, accounting for an important source of income for small farmers. The availability of the complete sequences of the two common bean genomes has dramatically accelerated and has enabled new experimental strategies to be applied for genetic research. DArTseq has been widely used as a method of SNP genotyping allowing comprehensive genome coverage with genetic applications in common bean breeding programs.

RESULTS

Using this technology, 6286 SNPs (1 SNP/86.5 Kbp) were genotyped in genic (43.3%) and non-genic regions (56.7%). Genetic subdivision associated to the common bean gene pools (K = 2) and related to grain types (K = 3 and K = 5) were reported. A total of 83% and 91% of all SNPs were polymorphic within the Andean and Mesoamerican gene pools, respectively, and 26% were able to differentiate the gene pools. Genetic diversity analysis revealed an average H E of 0.442 for the whole collection, 0.102 for Andean and 0.168 for Mesoamerican gene pools (F ST  = 0.747 between gene pools), 0.440 for the group of cultivars and lines, and 0.448 for the group of landrace accessions (F ST  = 0.002 between cultivar/line and landrace groups). The SNP effects were predicted with predominance of impact on non-coding regions (77.8%). SNPs under selection were identified within gene pools comparing landrace and cultivar/line germplasm groups (Andean: 18; Mesoamerican: 69) and between the gene pools (59 SNPs), predominantly on chromosomes 1 and 9. The LD extension estimate corrected for population structure and relatedness (r2SV) was ~ 88 kbp, while for the Andean gene pool was ~ 395 kbp, and for the Mesoamerican was ~ 130 kbp.

CONCLUSIONS

For common bean, DArTseq provides an efficient and cost-effective strategy of generating SNPs for large-scale genome-wide studies. The DArTseq resulted in an operational panel of 560 polymorphic SNPs in linkage equilibrium, providing high genome coverage. This SNP set could be used in genotyping platforms with many applications, such as population genetics, phylogeny relation between common bean varieties and support to molecular breeding approaches.

Pentaploid Wheat Hybrids: Applications, Characterisation, and Challenges

Interspecific hybridisation between hexaploid and tetraploid wheat species leads to the development of F1 pentaploid hybrids with unique chromosomal constitutions. Pentaploid hybrids derived from bread wheat (Triticum aestivum L.) and durum wheat (Triticum turgidum spp. durum Desf.) crosses can improve the genetic background of either parent by transferring traits of interest. The genetic variability derived from bread and durum wheat and transferred into pentaploid hybrids has the potential to improve disease resistance, abiotic tolerance, and grain quality, and to enhance agronomic characters. Nonetheless, pentaploid wheat hybrids have not been fully exploited in breeding programs aimed at improving crops. There are several potential barriers for efficient pentaploid wheat production, such as low pollen compatibility, poor seed set, failed seedling establishment, and frequent sterility in F1 hybrids. However, most of the barriers can be overcome by careful selection of the parental genotypes and by employing the higher ploidy level genotype as the maternal parent. In this review, we summarize the current research on pentaploid wheat hybrids and analyze the advantages and pitfalls of current methods used to assess pentaploid-derived lines. Furthermore, we discuss current and potential applications in commercial breeding programs and future directions for research into pentaploid wheat.

Swept away: ocean currents and seascape features influence genetic structure across the 18,000 Km Indo-Pacific distribution of a marine invertebrate, the black-lip pearl oyster Pinctada margaritifera

BACKGROUND

Genetic structure in many widely-distributed broadcast spawning marine invertebrates remains poorly understood, posing substantial challenges for their fishery management, conservation and aquaculture. Under the Core-Periphery Hypothesis (CPH), genetic diversity is expected to be highest at the centre of a species' distribution, progressively decreasing with increased differentiation towards outer range limits, as populations become increasingly isolated, fragmented and locally adapted. The unique life history characteristics of many marine invertebrates such as high dispersal rates, stochastic survival and variable recruitment are also likely to influence how populations are organised. To examine the microevolutionary forces influencing population structure, connectivity and adaptive variation in a highly-dispersive bivalve, populations of the black-lip pearl oyster Pinctada margaritifera were examined across its ~18,000 km Indo-Pacific distribution.

RESULTS

Analyses utilising 9,624 genome-wide SNPs and 580 oysters, discovered differing patterns of significant and substantial broad-scale genetic structure between the Indian and Pacific Ocean basins. Indian Ocean populations were markedly divergent (F st = 0.2534-0.4177, p < 0.001), compared to Pacific Ocean oysters, where basin-wide gene flow was much higher (F st = 0.0007-0.1090, p < 0.001). Partitioning of genetic diversity (hierarchical AMOVA) attributed 18.1% of variance between ocean basins, whereas greater proportions were resolved within samples and populations (45.8% and 35.7% respectively). Visualisation of population structure at selectively neutral loci resolved three and five discrete genetic clusters for the Indian and Pacific Oceans respectively. Evaluation of genetic structure at adaptive loci for Pacific populations (89 SNPs under directional selection; F st = 0.1012-0.4371, FDR = 0.05), revealed five clusters identical to those detected at neutral SNPs, suggesting environmental heterogeneity within the Pacific. Patterns of structure and connectivity were supported by Mantel tests of isolation by distance (IBD) and independent hydrodynamic particle dispersal simulations.

CONCLUSIONS

It is evident that genetic structure and connectivity across the natural range of P. margaritifera is highly complex, and produced by the interaction of ocean currents, IBD and seascape features at a broad scale, together with habitat geomorphology and local adaptation at regional levels. Overall population organisation is far more elaborate than generalised CPH predictions, however valuable insights for regional fishery management, and a greater understanding of range-wide genetic structure in a highly-dispersive marine invertebrate have been gained.

Genetic Variation and Population Structure of Oryza glaberrima and Development of a Mini-Core Collection Using DArTseq

The sequence variation present in accessions conserved in genebanks can best be used in plant improvement when it is properly characterized and published. Using low cost and high density single nucleotide polymorphism (SNP) assays, the genetic diversity, population structure, and relatedness between pairs of accessions can be quickly assessed. This information is relevant for different purposes, including creating core and mini-core sets that represent the maximum possible genetic variation contained in the whole collection. Here, we studied the genetic variation and population structure of 2,179 Oryza glaberrima Steud. accessions conserved at the AfricaRice genebank using 27,560 DArTseq-based SNPs. Only 14% (3,834 of 27,560) of the SNPs were polymorphic across the 2,179 accessions, which is much lower than diversity reported in other Oryza species. Genetic distance between pairs of accessions varied from 0.005 to 0.306, with 1.5% of the pairs nearly identical, 8.0% of the pairs similar, 78.1% of the pairs moderately distant, and 12.4% of the pairs very distant. The number of redundant accessions that contribute little or no new genetic variation to the O. glaberrima collection was very low. Using the maximum length sub-tree method, we propose a subset of 1,330 and 350 accessions to represent a core and mini-core collection, respectively. The core and mini-core sets accounted for ~61 and 16%, respectively, of the whole collection, and captured 97-99% of the SNP polymorphism and nearly all allele and genotype frequencies observed in the whole O. glaberrima collection available at the AfricaRice genebank. Cluster, principal component and model-based population structure analyses all divided the 2,179 accessions into five groups, based roughly on country of origin but less so on ecology. The first, third and fourth groups consisted of accessions primarily from Liberia, Nigeria, and Mali, respectively; the second group consisted primarily of accessions from Togo and Nigeria; and the fifth and smallest group was a mixture of accessions from multiple countries. Analysis of molecular variance showed between 10.8 and 28.9% of the variation among groups with the remaining 71.1-89.2% attributable to differences within groups.

Identification by the DArTseq method of the genetic origin of the Coffea canephora cultivated in Vietnam and Mexico

BACKGROUND

The coffee species Coffea canephora is commercially identified as "Conilon" when produced in Brazil, or "Robusta" when produced elsewhere in the world. It represents approximately 40 % of coffee production worldwide. While the genetic diversity of wild C. canephora has been well studied in the past, only few studies have addressed the genetic diversity of currently cultivated varieties around the globe. Vietnam is the largest Robusta producer in the world, while Mexico is the only Latin American country, besides Brazil, that has a significant Robusta production. Knowledge of the genetic origin of Robusta cultivated varieties in countries as important as Vietnam and Mexico is therefore of high interest.

RESULTS

Through the use of Sequencing-based diversity array technology-DArTseq method-on a collection of C. canephora composed of known accessions and accessions cultivated in Vietnam and Mexico, 4,021 polymorphic SNPs were identified. We used a multivariate analysis using SNP data from reference accessions in order to confirm and further fine-tune the genetic diversity of C. canephora. Also, by interpolating the data obtained for the varieties from Vietnam and Mexico, we determined that they are closely related to each other, and identified that their genetic origin is the Robusta Congo - Uganda group.

CONCLUSIONS

The genetic characterization based on SNP markers of the varieties grown throughout the world, increased our knowledge on the genetic diversity of C. canephora, and contributed to the understanding of the genetic background of varieties from very important coffee producers. Given the common genetic origin of the Robusta varieties cultivated in Vietnam, Mexico and Uganda, and the similar characteristics of climatic areas and relatively high altitude where they are grown, we can state that the Vietnamese and the Mexican Robusta have the same genetic potential to produce good cup quality.

The Development of Quality Control Genotyping Approaches: A Case Study Using Elite Maize Lines

Quality control (QC) of germplasm identity and purity is a critical component of breeding and conservation activities. SNP genotyping technologies and increased availability of markers provide the opportunity to employ genotyping as a low-cost and robust component of this QC. In the public sector available low-cost SNP QC genotyping methods have been developed from a very limited panel of markers of 1,000 to 1,500 markers without broad selection of the most informative SNPs. Selection of optimal SNPs and definition of appropriate germplasm sampling in addition to platform section impact on logistical and resource-use considerations for breeding and conservation applications when mainstreaming QC. In order to address these issues, we evaluated the selection and use of SNPs for QC applications from large DArTSeq data sets generated from CIMMYT maize inbred lines (CMLs). Two QC genotyping strategies were developed, the first is a “rapid QC”, employing a small number of SNPs to identify potential mislabeling of seed packages or plots, the second is a “broad QC”, employing a larger number of SNP, used to identify each germplasm entry and to measure heterogeneity. The optimal marker selection strategies combined the selection of markers with high minor allele frequency, sampling of clustered SNP in proportion to marker cluster distance and selecting markers that maintain a uniform genomic distribution. The rapid and broad QC SNP panels selected using this approach were further validated using blind test assessments of related re-generation samples. The influence of sampling within each line was evaluated. Sampling 192 individuals would result in close to 100% possibility of detecting a 5% contamination in the entry, and approximately a 98% probability to detect a 2% contamination of the line. These results provide a framework for the establishment of QC genotyping. A comparison of financial and time costs for use of these approaches across different platforms is discussed providing a framework for institutions involved in maize conservation and breeding to assess the resource use effectiveness of QC genotyping. Application of these research findings, in combination with existing QC approaches, will ensure the regeneration, distribution and use in breeding of true to type inbred germplasm. These findings also provide an effective approach to optimize SNP selection for QC genotyping in other species.

The application of GBS markers for extending the dense genetic map of rye (Secale cereale L.) and the localization of the Rfc1 gene restoring male fertility in plants with the C source of sterility-inducing cytoplasm

Genotyping by sequencing (GBS) is an efficient method of genotyping in numerous plant species. One of the crucial steps toward the application of GBS markers in crop improvement is anchoring them on particular chromosomes. In rye (Secale cereale L.), chromosomal localization of GBS markers has not yet been reported. In this paper, the application of GBS markers generated by the DArTseq platform for extending the high-density map of rye is presented. Additionally, their application is used for the localization of the Rfc1 gene that restores male fertility in plants with the C source of sterility-inducing cytoplasm. The total number of markers anchored on the current version of the map is 19,081, of which 18,132 were obtained from the DArTseq platform. Numerous markers co-segregated within the studied mapping population, so, finally, only 3397 unique positions were located on the map of all seven rye chromosomes. The total length of the map is 1593 cM and the average distance between markers is 0.47 cM. In spite of the resolution of the map being not very high, it should be a useful tool for further studies of the Secale cereale genome because of the presence on this map of numerous GBS markers anchored for the first time on rye chromosomes. The Rfc1 gene was located on high-density maps of the long arm of the 4R chromosome obtained for two mapping populations. Genetic maps were composed of DArT, DArTseq, and PCR-based markers. Consistent mapping results were obtained and DArTs tightly linked to the Rfc1 gene were successfully applied for the development of six new PCR-based markers useful in marker-assisted selection.

Genetic dissection of maize plant architecture with an ultra-high density bin map based on recombinant inbred lines

Background

Plant architecture attributes, such as plant height, ear height, and internode number, have played an important role in the historical increases in grain yield, lodging resistance, and biomass in maize (Zea mays L.). Analyzing the genetic basis of variation in plant architecture using high density QTL mapping will be of benefit for the breeding of maize for many traits. However, the low density of molecular markers in existing genetic maps has limited the efficiency and accuracy of QTL mapping. Genotyping by sequencing (GBS) is an improved strategy for addressing a complex genome via next-generation sequencing technology. GBS has been a powerful tool for SNP discovery and high-density genetic map construction. The creation of ultra-high density genetic maps using large populations of advanced recombinant inbred lines (RILs) is an efficient way to identify QTL for complex agronomic traits.

Results

A set of 314 RILs derived from inbreds Ye478 and Qi319 were generated and subjected to GBS. A total of 137,699,000 reads with an average of 357,376 reads per individual RIL were generated, which is equivalent to approximately 0.07-fold coverage of the maize B73 RefGen_V3 genome for each individual RIL. A high-density genetic map was constructed using 4183 bin markers (100-Kb intervals with no recombination events). The total genetic distance covered by the linkage map was 1545.65 cM and the average distance between adjacent markers was 0.37 cM with a physical distance of about 0.51 Mb. Our results demonstrated a relatively high degree of collinearity between the genetic map and the B73 reference genome. The quality and accuracy of the bin map for QTL detection was verified by the mapping of a known gene, pericarp color 1 (P1), which controls the color of the cob, with a high LOD value of 80.78 on chromosome 1. Using this high-density bin map, 35 QTL affecting plant architecture, including 14 for plant height, 14 for ear height, and seven for internode number were detected across three environments. Interestingly, pQTL10, which influences all three of these traits, was stably detected in three environments on chromosome 10 within an interval of 14.6 Mb. Two MYB transcription factor genes, GRMZM2G325907 and GRMZM2G108892, which might regulate plant cell wall metabolism are the candidate genes for qPH10.

Conclusions

Here, an ultra-high density accurate linkage map for a set of maize RILs was constructed using a GBS strategy. This map will facilitate identification of genes and exploration of QTL for plant architecture in maize. It will also be helpful for further research into the mechanisms that control plant architecture while also providing a basis for marker-assisted selection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2555-z) contains supplementary material, which is available to authorized users.

Development of cleaved amplified polymorphic sequence markers and a CAPS-based genetic linkage map in watermelon (Citrullus lanatus [Thunb.] Matsum. and Nakai) constructed using whole-genome re-sequencing data

Cleaved amplified polymorphic sequence (CAPS) markers are useful tools for detecting single nucleotide polymorphisms (SNPs). This study detected and converted SNP sites into CAPS markers based on high-throughput re-sequencing data in watermelon, for linkage map construction and quantitative trait locus (QTL) analysis. Two inbred lines, Cream of Saskatchewan (COS) and LSW-177 had been re-sequenced and analyzed by Perl self-compiled script for CAPS marker development. 88.7% and 78.5% of the assembled sequences of the two parental materials could map to the reference watermelon genome, respectively. Comparative assembled genome data analysis provided 225,693 and 19,268 SNPs and indels between the two materials. 532 pairs of CAPS markers were designed with 16 restriction enzymes, among which 271 pairs of primers gave distinct bands of the expected length and polymorphic bands, via PCR and enzyme digestion, with a polymorphic rate of 50.94%. Using the new CAPS markers, an initial CAPS-based genetic linkage map was constructed with the F2 population, spanning 1836.51 cM with 11 linkage groups and 301 markers. 12 QTLs were detected related to fruit flesh color, length, width, shape index, and brix content. These newly CAPS markers will be a valuable resource for breeding programs and genetic studies of watermelon.

Updated sesame genome assembly and fine mapping of plant height and seed coat color QTLs using a new high-density genetic map

Background

Sesame is an important high-quality oil seed crop. The sesame genome was de novo sequenced and assembled in 2014 (version 1.0); however, the number of anchored pseudomolecules was higher than the chromosome number (2n = 2x = 26) due to the lack of a high-density genetic map with 13 linkage groups.

Results

We resequenced a permanent population consisting of 430 recombinant inbred lines and constructed a genetic map to improve the sesame genome assembly. We successfully anchored 327 scaffolds onto 13 pseudomolecules. The new genome assembly (version 2.0) included 97.5 % of the scaffolds greater than 150 kb in size present in assembly version 1.0 and increased the total pseudomolecule length from 233.7 to 258.4 Mb with 94.3 % of the genome assembled and 97.2 % of the predicted gene models anchored. Based on the new genome assembly, a bin map including 1,522 bins spanning 1090.99 cM was generated and used to identified 41 quantitative trait loci (QTLs) for sesame plant height and 9 for seed coat color. The plant height-related QTLs explained 3–24 % the phenotypic variation (mean value, 8 %), and 29 of them were detected in at least two field trials. Two major loci (qPH-8.2 and qPH-3.3) that contributed 23 and 18 % of the plant height were located in 350 and 928-kb spaces on Chr8 and Chr3, respectively. qPH-3.3, is predicted to be responsible for the semi-dwarf sesame plant phenotype and contains 102 candidate genes. This is the first report of a sesame semi-dwarf locus and provides an interesting opportunity for a plant architecture study of the sesame. For the sesame seed coat color, the QTLs of the color spaces L*, a*, and b* were detected with contribution rates of 3–46 %. qSCb-4.1 contributed approximately 39 % of the b* value and was located on Chr4 in a 199.9-kb space. A list of 32 candidate genes for the locus, including a predicted black seed coat-related gene, was determined by screening the newly anchored genome.

Conclusions

This study offers a high-density genetic map and an improved assembly of the sesame genome. The number of linkage groups and pseudomolecules in this assembly equals the number of sesame chromosomes for the first time. The map and updated genome assembly are expected to serve as a platform for future comparative genomics and genetic studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2316-4) contains supplementary material, which is available to authorized users.

Genome-Wide Differentiation of Various Melon Horticultural Groups for Use in GWAS for Fruit Firmness and Construction of a High Resolution Genetic Map

Melon (Cucumis melo L.) is a phenotypically diverse eudicot diploid (2n = 2x = 24) has climacteric and non-climacteric morphotypes and show wide variation for fruit firmness, an important trait for transportation and shelf life. We generated 13,789 SNP markers using genotyping-by-sequencing (GBS) and anchored them to chromosomes to understand genome-wide fixation indices (Fst) between various melon morphotypes and genomewide linkage disequilibrium (LD) decay. The F_ST between accessions of cantalupensis and inodorus was 0.23. The F_ST between cantalupensis and various agrestis accessions was in a range of 0.19-0.53 and between inodorus and agrestis accessions was in a range of 0.21-0.59 indicating sporadic to wide ranging introgression. The EM (Expectation Maximization) algorithm was used for estimation of 1436 haplotypes. Average genome-wide LD decay for the melon genome was noted to be 9.27 Kb. In the current research, we focused on the genome-wide divergence underlying diverse melon horticultural groups. A high-resolution genetic map with 7153 loci was constructed. Genome-wide segregation distortion and recombination rate across various chromosomes were characterized. Melon has climacteric and non-climacteric morphotypes and wide variation for fruit firmness, a very important trait for transportation and shelf life. Various levels of QTLs were identified with high to moderate stringency and linked to fruit firmness using both genome-wide association study (GWAS) and biparental mapping. Gene annotation revealed some of the SNPs are located in β-D-xylosidase, glyoxysomal malate synthase, chloroplastic anthranilate phosphoribosyltransferase, and histidine kinase, the genes that were previously characterized for fruit ripening and softening in other crops.