Advancements in genetic techniques and functional genomics for enhancing crop traits and agricultural sustainability

Genetic variability is essential for the development of new crop varieties with economically beneficial traits. The traits can be inherited from wild relatives or induced through mutagenesis. Novel genetic elements can then be identified and new gene functions can be predicted. In this study, forward and reverse genetics approaches were described, in addition to their applications in modern crop improvement programs and functional genomics. By using heritable phenotypes and linked genetic markers, forward genetics searches for genes by using traditional genetic mapping and allele frequency estimation. Despite recent advances in sequencing technology, omics and computation, genetic redundancy remains a major challenge in forward genetics. By analyzing close-related genes, we will be able to dissect their functional redundancy and predict possible traits and gene activity patterns. In addition to these predictions, sophisticated reverse gene editing tools can be used to verify them, including TILLING, targeted insertional mutagenesis, gene silencing, gene targeting and genome editing. By using gene knock-down, knock-up and knock-out strategies, these tools are able to detect genetic changes in cells. In addition, epigenome analysis and editing enable the development of novel traits in existing crop cultivars without affecting their genetic makeup by increasing epiallelic variants. Our understanding of gene functions and molecular dynamics of various biological phenomena has been revised by all of these findings. The study also identifies novel genetic targets in crop species to improve yields and stress tolerances through conventional and non-conventional methods. In this article, genetic techniques and functional genomics are specifically discussed and assessed for their potential in crop improvement.

Unveiling the molecular regulatory mechanisms underlying sucrose accumulation and oil reduction in peanut kernels through genetic mapping and transcriptome analysis

Sucrose content is a key factor for the flavor of edible peanut, which determines the sweet taste of fresh peanut and also attribute to pleasant flavor of roasted peanut. To explore the genetic mechanism of the sucrose content in peanut, an F2 population was created by crossing the sweet cultivar Zhonghuatian 1 (ZHT1) with Nanyangbaipi (NYBP). A genomic region spanning 28.26 kb on chromosome A06 was identified for the sucrose content through genetic mapping, elucidating 47.5% phenotypic variance explained. As the sucrose content had a significantly negative correlation with the oil content, this region was also found to be related to the oil content explaining 37.2% of phenotype variation. In this region, Arahy.42CAD1 was characterized as the most likely candidate gene through a comprehensive analysis. The nuclear localization of Arahy.42CAD1 suggests its potential involvement in the regulation of gene expression for sucrose and oil contents in peanut. Transcriptome analysis of the developing seeds in both parents revealed that genes involved in glycolysis and triacylglycerol biosynthesis pathways were not significantly down-regulated in ZHT1, indicating that the sucrose accumulation was not attributed to the suppression of triacylglycerol biosynthesis. Based on the WGCNA analysis, Arahy.42CAD1 was co-expressed with the genes involved in vesicle transport and oil body assembly, suggesting that the sucrose accumulation may be caused by disruptions in TAG transportation or storage mechanisms. These findings offer new insights into the molecular mechanisms governing sucrose accumulation in peanut, and also provide a potential gene target for enhancing peanut flavor.

Contrasting roles of cytochrome P450s in amitraz and chlorfenapyr resistance in the crop pest Tetranychus urticae

The molecular mechanisms of amitraz and chlorfenapyr resistance remain only poorly understood for major agricultural pests and vectors of human diseases. This study focusses on a multi-resistant field strain of the crop pest Tetranychus urticae, which could be readily selected in the laboratory to high levels of amitraz and chlorfenapyr resistance. Toxicity experiments using tralopyril, the active toxophore of chlorfenapyr, suggested decreased activation as a likely mechanism underlying resistance. Starting from the same parental strain, transcriptome profiling revealed that a cluster of detoxifying genes was upregulated after amitraz selection, but unexpectedly downregulated after chlorfenapyr selection. Further functional validation associated the upregulation of CYP392A16 with amitraz metabolism and the downregulation of CYP392D8 with reduced activation of chlorfenapyr to tralopyril. Genetic mapping (QTL analysis by BSA) was conducted in an attempt to unravel the genetic mechanisms of expression variation and resistance. This revealed that chlorfenapyr resistance was associated with a single QTL, while 3 QTLs were uncovered for amitraz resistance. Together with the observed contrasting gene expression patterns, we argue that transcriptional regulators most likely underly the distinct expression profiles associated with resistance, but these await further functional validation.

Opioid trail: Tracking contributions to opioid use disorder from host genetics to the gut microbiome

Opioid use disorder (OUD) is a worldwide public health crisis with few effective treatment options. Traditional genetics and neuroscience approaches have provided knowledge about biological mechanisms that contribute to OUD-related phenotypes, but the complexity and magnitude of effects in the brain and body remain poorly understood. The gut-brain axis has emerged as a promising target for future therapeutics for several psychiatric conditions, so characterizing the relationship between host genetics and the gut microbiome in the context of OUD will be essential for development of novel treatments. In this review, we describe evidence that interactions between host genetics, the gut microbiome, and immune signaling likely play a key role in mediating opioid-related phenotypes. Studies in humans and model organisms consistently demonstrated that genetic background is a major determinant of gut microbiome composition. Furthermore, the gut microbiome is susceptible to environmental influences such as opioid exposure. Additional work focused on gene by microbiome interactions will be necessary to gain improved understanding of their effects on OUD-related behaviors.

Efficient targeted recombination with CRISPR/Cas9 in hybrids of Caenorhabditis nematodes with suppressed recombination

BACKGROUND

Homology-based recombination (HR) is the cornerstone of genetic mapping. However, a lack of sufficient sequence homology or the presence of a genomic rearrangement prevents HR through crossing, which inhibits genetic mapping in relevant genomic regions. This is particularly true in species hybrids whose genomic sequences are highly divergent along with various genome arrangements, making the mapping of genetic loci, such as hybrid incompatibility (HI) loci, through crossing impractical. We previously mapped tens of HI loci between two nematodes, Caenorhabditis briggsae and C. nigoni, through the repeated backcrossing of GFP-linked C. briggsae fragments into C. nigoni. However, the median introgression size was over 7 Mb, indicating apparent HR suppression and preventing the subsequent cloning of the causative gene underlying a given HI phenotype. Therefore, a robust method that permits recombination independent of sequence homology is desperately desired.

RESULTS

Here, we report a method of highly efficient targeted recombination (TR) induced by CRISPR/Cas9 with dual guide RNAs (gRNAs), which circumvents the HR suppression in hybrids between the two species. We demonstrated that a single gRNA was able to induce efficient TR between highly homologous sequences only in the F1 hybrids but not in the hybrids that carry a GFP-linked C. briggsae fragment in an otherwise C. nigoni background. We achieved highly efficient TR, regardless of sequence homology or genetic background, when dual gRNAs were used that each specifically targeted one parental chromosome. We further showed that dual gRNAs were able to induce efficient TR within genomic regions that had undergone inversion, in which HR-based recombination was expected to be suppressed, supporting the idea that dual-gRNA-induced TR can be achieved through nonhomology-based end joining between two parental chromosomes.

CONCLUSIONS

Recombination suppression can be circumvented through CRISPR/Cas9 with dual gRNAs, regardless of sequence homology or the genetic background of the species hybrid. This method is expected to be applicable to other situations in which recombination is suppressed in interspecies or intrapopulation hybrids.

A review of the molecular mechanisms of acaricide resistance in mites and ticks

The Arachnida subclass of Acari comprises many harmful pests that threaten agriculture as well as animal health, including herbivorous spider mites, the bee parasite Varroa, the poultry mite Dermanyssus and several species of ticks. Especially in agriculture, acaricides are often used intensively to minimize the damage they inflict, promoting the development of resistance. Beneficial predatory mites used in biological control are also subjected to acaricide selection in the field. The development and use of new genetic and genomic tools such as genome and transcriptome sequencing, bulked segregant analysis (QTL mapping), and reverse genetics via RNAi or CRISPR/Cas9, have greatly increased our understanding of the molecular genetic mechanisms of resistance in Acari, especially in the spider mite Tetranychus urticae which emerged as a model species. These new techniques allowed to uncover and validate new resistance mutations in a larger range of species. In addition, they provided an impetus to start elucidating more challenging questions on mechanisms of gene regulation of detoxification associated with resistance.

Natural variation and genetic loci underlying resistance to grain shattering in standing crop of modern wheat

Modern wheat (Triticum aestivum L.) cultivars have a free-threshing habit, which allows for easy manual or mechanical threshing. However, when harvesting is delayed or extreme weather events occur at harvest time, grain shattering can cause severe loss of harvestable yield. In the past, grain size was considered a predisposing factor as large, plump kernels can lead to buckling and breaking of the outer glume, but the correlation between glume strength and shattering is not strong in modern wheat, and it is hypothesised that there may be other genetic mechanisms. Data from two bi-parent populations and a wheat diversity panel were analyzed to explore the underlying genetic basis for grain shattering observed in multiple field experiments through quantitative trait loci (QTL) analysis. Grain shattering had a significant and negative association with grain yield, irrespective of populations and environments. The correlation with plant height was positive in all populations, but correlations with phenology were population specific, being negative in the diversity panel and the Drysdale × Waagan population, and positive in the Crusader × RT812 population. In the wheat diversity panel, allelic variations at well-known major genes (Rht-B1, Rht-D1 and Ppd-D1) showed minimal association with grain shattering. Instead, the genome-wide analysis identified a single locus on chromosome 2DS, which explained 50% of the phenotypic variation, and mapping to ~ 10 Mb from Tenacious glume (Tg) gene. In the Drysdale × Waagan cross, however, the reduced height (Rht) genes showed major effects on grain shattering. At the Rht-B1 locus, the Rht-B1b allele was associated with 10.4 cm shorter plant height, and 18% decreased grain shattering, whereas Rht-D1b reduced plant height by 11.4 cm and reduced grain shattering by 20%. Ten QTL were detected in the Crusader × RT812, including a major locus detected on the long arm of chromosome 5A. All the QTL identified in this population were non-pleiotropic, as they were still significant even after removing the influence of plant height. In conclusion, these results indicated a complex genetic system for grain shattering in modern wheat, which varied with genetic background, involved pleiotropic as well as independent gene actions, and which might be different from shattering in wild wheat species caused by major domestication genes. The influence of Rht genes was confirmed, and this provides valuable information in breeding crops of the future. Further, the SNP marker close to Tg on chromosome 2DS should be considered for utility in marker-assisted selection.

Identification and Genomic Localization of Autosomal sdY Locus in a Population of Atlantic Salmon (Salmo salar)

The determination of sex in salmonid fishes is controlled by genetic mechanisms, with males being the heterogametic sex. The master sex-determining gene, the sexually dimorphic gene on the Y chromosome (sdY), is a conserved gene across various salmonid species. Nevertheless, variations in the genomic location of sdY have been observed both within and between species. Furthermore, different studies have reported discordances in the association between the sdY and the phenotypic gender. While some males seem to lack this locus, there have been reports of females carrying sdY. Although the exact reasons behind this discordance remain under investigation, some recent studies have proposed the existence of an autosomal, non-functional copy of sdY as a potential cause. In this study, we confirmed the presence of this autosomal sdY in the SalmoBreed strain of Atlantic salmon using a genotyping platform through a novel approach that allows for high-throughput screening of a large number of individuals. We further characterized the segregation profile of this locus across families and found the ratio of genetically assigned female-to-male progeny to be in accordance with the expected profile of a single autosomal sdY locus. Additionally, our mapping efforts localized this locus to chromosome 3 and suggested a putative copy on chromosome 6.

Identification and high-resolution mapping of a novel tiller number gene (tin6) by combining forward genetics screen and MutMap approach in bread wheat

Wheat (Triticum aestivum) is one of the most important food crops worldwide, providing up to 20% of the caloric intake per day. Developing high-yielding wheat cultivars with tolerance against abiotic and biotic stresses is important to keep up with the increasing human population. Tiller number is one of the major yield-related traits, directly affecting the number of grains produced per plant; however, only a small number of QTL and underlining genes have been identified for this important factor. Identification of novel genetic variation underlying contrasting traits and their precise genetic mapping in wheat is considered difficult due to the complexity and size of the genome; however, advancements in genomic resources have made efficient gene localization more possible. In this study, we report the characterization of a novel tillering number gene using a mutant identified in the forward genetic screen of an ethyl methane sulfonate (EMS)-treated population of cv. "Jagger." By crossing the low tillering mutant with the Jagger wild-type plant, we generated an F2 population and used the MutMap approach to identify a novel physical interval on 11 Mb on chromosome 2DS. Using an F₂ population of 442 gametes and polymorphic SNP markers, we were able to delineate the tin6 locus to a 2.1 Mb region containing 22 candidate genes.

Structural variation among assembled genomes facilitates development of rapid and low-cost NOR-linked markers and NOR-telomere junction mapping in Arabidopsis

KEY MESSAGE

Genome-wide structural variants we identified and new NOR-linked markers we developed would be useful for future genome-wide association studies (GWAS), and for new gene/trait mapping purposes. Bioinformatic alignment of the assembled genomes of Col-0 and Sha ecotypes of Arabidopsis thaliana revealed ~ 13,000 genome-wide structural variants involving simple insertions or deletions and repeat contractions or expansions. Using some of these structural variants, we developed new, rapid, and low-cost PCR-based molecular markers that are genetically linked to the nucleolus organizer regions (NORs). A. thaliana has two NORs, one each on chromosome 2 (NOR2) and chromosome 4 (NOR4). Both NORs are ~ 4 Mb each, and hundreds of 45S ribosomal RNA (rRNA) genes are tandemly arrayed at these loci. Using previously characterized recombinant inbred lines (RILs) derived from Sha x Col-0  crosses, we validated the utility of the newly developed NOR-linked markers in genetically mapping rRNA genes and the associated telomeres to either NOR2 or NOR4. Lastly, we sequenced Sha genome using Oxford Nanopore Technology (ONT) and used the data to obtain sequences of NOR-telomere junctions, and with the help of RILs, we mapped them as new genetic markers to their respective NORs (NOR2-TEL2N and NOR4-TEL4N). The structural variants obtained from this study would serve as valuable data for genome-wide association studies (GWAS), and to rapidly design more genome-wide genetic (molecular) markers for new gene/trait mapping purposes.

Mapping of a soybean rust resistance in PI 594756 at the Rpp1 locus

Asian soybean rust (ASR), caused by the fungus Phakopsora pachyrhizi, is the main disease affecting soybean in Brazil. This study aimed at investigating and mapping the resistance of the PI 594756 to P. pachyrhizi, by using Bulked Segregant Analysis (BSA). The PI 594756 and the susceptible PI 594891 were crossed and the resulting F2 and F2:3 populations (208 and 1770 plants, respectively) were tested against ASR. Also, these PIs and differential varieties were tested against a panel of monosporic isolates. Plants presenting tan lesions were classified as susceptible (S) while plants presenting reddish-brown (RB) lesions were classified as resistant. DNA bulks were genotyped with Infinium BeadChips and the genomic region identified was further analyzed in the F2 individuals with target GBS (tGBS). PI 594,56 presented a unique resistance profile compared to the differential varieties. The resistance was monogenic dominant; however, it was classified as incompletely dominant when quantitatively studied. Genetic and QTL mapping placed the PI 594756 gene between the genomic region located at 55,863,741 and 56,123,516 bp of chromosome 18. This position is slightly upstream mapping positions of Rpp1 (PI 200492) and Rpp1-b (PI 594538A). Finally, we performed a haplotype analysis in a whole genomic sequencing-SNP database composed of Brazilian historical germplasm and sources of Rpp genes. We found SNPs that successfully differentiated the new PI 594756 allele from Rpp1 and Rpp1-b sources. The haplotype identified can be used as a tool for marker-assisted selection (MAS).

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01358-4.

Genetic architecture of fresh-market tomato yield

BACKGROUND

The fresh-market tomato (Solanum lycopersicum) is bred for direct consumption and is selected for a high yield of large fruits. To understand the genetic variations (distinct types of DNA sequence polymorphism) that influence the yield, we collected the phenotypic variations in the yields of total fruit, extra-large-sized fruit, small-sized fruit, or red-colored fruit from 68 core inbred contemporary U.S. fresh-market tomatoes for three consecutive years and the genomic information in 8,289,741 single nucleotide polymorphism (SNP) positions from the whole-genome resequencing of these tomatoes.

RESULTS

Genome-wide association (GWA) mapping using the SNP data with or without SNP filtering steps using the regularization methods, validated with quantitative trait loci (QTL) linkage mapping, identified 18 significant association signals for traits evaluated. Among them, 10 of which were not located within genomic regions previously identified as being associated with fruit size/shape. When mapping-driven association signals [558 SNPs associated with 28 yield (component) traits] were used to calculate genomic estimated breeding values (GEBVs) of evaluated traits, the prediction accuracies of the extra-large-sized fruit and small-sized fruit yields were higher than those of the total and red-colored fruit yields, as we tested the generated breeding values in inbred tomatoes and F2 populations. Improved accuracy in GEBV calculation of evaluated traits was achieved by using 364 SNPs identified using the regularization methods.

CONCLUSIONS

Together, these results provide an understanding of the genetic variations underlying the heritable phenotypic variability in yield in contemporary tomato breeding and the information necessary for improving such economically important and complex quantitative trait through breeding.

Gamete Binning to Achieve Haplotype-Resolved Genome Assembly

Haplotype-resolved genome assemblies remain a challenge in practice. Here, we provide a step-by-step guide on gamete binning, a method to generate haplotype-resolved genome assemblies for diploid species. The protocol starts by phasing heterozygous variants to individual haplotypes of specific chromosomes using the genome information of individual haploid gametes of the focal individual. Using phased variants, the whole-genome sequencing reads from the diploid genome can be genotyped and assigned into groups, which represent the individual haplotypes of each of the chromosomes. Finally, haplotype-specific chromosomes can be assembled independently using standard assembly tools. First applications of gamete binning revealed a haplotyping accuracy over 99%, which outperformed sequence-only or Hi-C-based haplotype-resolved genome assemblies.Availability: github.com/schneebergerlab/GameteBinning_prac .

Genetic mapping of maize metabolites using high-throughput mass profiling

Plant metabolites are the basis of human nutrition and have biological relevance in ecology. Farmers selected plants with favorable characteristics since prehistoric times and improved the cultivars, but without knowledge of underlying mechanisms. Understanding the genetic basis of metabolite production can facilitate the successful breeding of plants with augmented nutritional value. To identify genetic factors related to the metabolic composition in maize, we generated mass profiles of 198 recombinant inbred lines (RILs) and their parents (B73 and Mo17) using direct-injection electrospray ionization mass spectrometry (DLI-ESI MS). Mass profiling allowed the correct clustering of samples according to genotype. We quantified 71 mass features from grains and 236 mass features from leaf extracts. For the corresponding ions, we identified tissue-specific metabolic 'Quantitative Trait Loci' (mQTLs) distributed across the maize genome. These genetic regions could regulate multiple metabolite biosynthesis pathways. Our findings demonstrate that DLI-ESI MS has sufficient analytical resolution to map mQTLs. These identified genetic loci will be helpful in metabolite-focused maize breeding. Mass profiling is a powerful tool for detecting mQTLs in maize and enables the high-throughput screening of loci responsible for metabolite biosynthesis.

In search of the relationship between the rye polyamine oxidase (PAO) gene and resistance to powdery mildew (PM)

Powdery mildew (PM), a common cereal disease in cultivated areas, including Europe and other temperate regions, is caused by the fungus Blumeria graminis. While PM is one of the most important wheat leaf diseases globally, rye is highly tolerant to PM. It has been reported that in barley infected with PM, polyamine oxidase (PAO) activity related to the production of hydrogen peroxide (H₂O₂) has increased, which may promote defense against biotrophic or hemibiotrophic pathogens. The current study aimed to assess the relationship between the segregation of the polymorphic marker for rye PAO (ScPAO) and the level of PM infection in plants. The genetic mapping in two interline populations shows that ScPAO is located on chromosome 7R. Further analysis comparing ScPAO location to mapped wheat (Triticum aestivum L.) PAO duplicates suggests the ScPAO homology with TaPAO6 or TaPAO7. A possible association of ScPAO from 7R with PM resistance is demonstrated in the recombinant inbred lines (RIL)-L population phenotyped for PM infection. Finally, three novel QTLs for PM resistance on the 7R chromosome of rye are detected.

Genome-wide association study and genetic mapping of BhWAX conferring mature fruit cuticular wax in wax gourd

BACKGROUND

Wax gourd [Benincasa hispida (Thunb) Cogn. (2n = 2x = 24)] is an economically important vegetable crop of genus Benincasa in the Cucurbitaceae family. Fruit is the main consumption organ of wax gourd. The mature fruit cuticular wax (MFCW) is an important trait in breeding programs, which is also of evolutionary significance in wax gourd. However, the genetic architecture of this valuable trait remains unrevealed.

RESULTS

In this study, genetic analysis revealed that the inheritance of MFCW was controlled by a single gene, with MFCW dominant over non-MFCW, and the gene was primarily named as BhWAX. Genome-wide association study (GWAS) highlighted a 1.1 Mb interval on chromosome 9 associated with MFCW in wax gourd germplasm resources. Traditional fine genetic mapping delimited BhWAX to a 0.5 Mb region containing 12 genes. Based on the gene annotation, expression analysis and co-segregation analysis, Bhi09G001428 that encodes a membrane bound O-acyltransferase (MBOAT) was proposed as the candidate gene for BhWAX. Moreover, it was demonstrated that the efficiency of a cleaved amplified polymorphic sequences (CAPS) marker in the determination of MFCW in wax gourd reached 80%.

CONCLUSIONS

In closing, the study identified the candidate gene controlling MFCW and provided an efficient molecular marker for the trait in wax gourd for the first time, which will not only be beneficial for functional validation of the gene and marker-assisted breeding of wax gourd, but also lay a foundation for analysis of its evolutionary meaning among cucurbits.

Heat Stress Tolerance 2 confers basal heat stress tolerance in allohexaploid wheat (Triticum aestivum L.)

Heat stress substantially reduces the yield potential of wheat (Triticum aestivum L.), one of the most widely cultivated staple crops, and greatly threatens global food security in the context of global warming. However, few studies have explored the heat stress tolerance (HST)-related genetic resources in wheat. Here, we identified and fine-mapped a wheat HST locus, TaHST2, which is indispensable for HST in both the vegetative and reproductive stages of the wheat life cycle. The studied pair of near isogenic lines (NILs) exhibited diverse morphologies under heat stress, based on which we mapped TaHST2 to a 485 kb interval on chromosome arm 4DS. Under heat stress, TaHST2 confers a superior conversion rate from soluble sugars to starch in wheat grains, resulting in faster grain filling and a higher yield potential. A further exploration of genetic resources indicated that TaHST2 underwent strong artificial selection during wheat domestication, suggesting it is an essential locus for basal HST in wheat. Our findings provide deeper insights into the genetic basis of wheat HST and might be useful for global efforts to breed heat-stress-tolerant cultivars.

A SNP mutation in the CsCLAVATA1 leads to pleiotropic variation in plant architecture and fruit morphogenesis in cucumber (Cucumis sativus L.)

Plant architectures is predominantly determined by branching pattern, internode elongation, phyllotaxis, shoot determinacy and reproductive organs. Domestication or improvement of this critical agronomic trait played an important role in the breakthrough of crop yield. Here, we identified a mutant with fasciated plant architecture, named fas, from an ethyl methanesulfonate (EMS) induced mutant population in cucumber. The mutant exhibited abnormal phyllotaxy, flattened main stem, increased number of floral organs, and significantly shorter and thicker fruits. However, the molecular mechanism conferring this pleiotropic effect remains unknown. Using a map-based cloning strategy, we isolated the gene CsaV3_3G045960, encoding a leucine-rich repeat receptor-like kinase, a putative direct homolog of the Arabidopsis CLAVATA1 protein referred to as CsCLV1. Endogenous hormone assays showed that IAA and GA₃ levels in fas stems and ovaries were significantly reduced. Conformably, RNA-seq analysis showed that CsCLV1 regulates cucumber stem and ovary development by coordinating hormones and transcription factors. Our results contribute to the understanding of the function of CsCLV1 throughout the growth cycle, provide new evidence that the CLV signaling system is functionally conserved in Cucurbitaceae.

inGAP-family: Accurate Detection of Meiotic Recombination Loci and Causal Mutations by Filtering Out Artificial Variants due to Genome Complexities

Accurately identifying DNA polymorphisms can bridge the gap between phenotypes and genotypes and is essential for molecular marker assisted genetic studies. Genome complexities, including large-scale structural variations, bring great challenges to bioinformatic analysis for obtaining high-confidence genomic variants, as sequence differences between non-allelic loci of two or more genomes can be misinterpreted as polymorphisms. It is important to correctly filter out artificial variants to avoid false genotyping or estimation of allele frequencies. Here, we present an efficient and effective framework, inGAP-family, to discover, filter, and visualize DNA polymorphisms and structural variants (SVs) from alignment of short reads. Applying this method to polymorphism detection on real datasets shows that elimination of artificial variants greatly facilitates the precise identification of meiotic recombination points as well as causal mutations in mutant genomes or quantitative trait loci. In addition, inGAP-family provides a user-friendly graphical interface for detecting polymorphisms and SVs, further evaluating predicted variants and identifying mutations related to genotypes. It is accessible at https://sourceforge.net/projects/ingap-family/.

Genetic mapping and physiological analysis of chlorophyll-deficient mutant in Brassica napus L

BACKGROUND

Leaf color mutants have reduced photosynthetic efficiency, which has severely negative impacts on crop growth and economic product yield. There are different chlorophyll mutants in Arabidopsis and crops that can be used for genetic control and molecular mechanism studies of chlorophyll biosynthesis, chloroplast development and photoefficiency. Chlorophyll mutants in Brassica napus are mostly used for mapping and location research but are rarely used for physiological research. The chlorophyll-deficient mutant in this experiment were both genetically mapped and physiologically analyzed.

RESULTS

In this study, yellow leaf mutant of Brassica napus L. mutated by ethyl methyl sulfone (EMS) had significantly lower chlorophyll a, b and carotenoid contents than the wild type, and the net photosynthetic efficiency, stomatal conductance and transpiration rate were all significantly reduced. The mutant had sparse chloroplast distribution and weak autofluorescence. The granule stacks were reduced, and the shape was extremely irregular, with more broken stromal lamella. Transcriptome data analysis enriched the differentially expressed genes mainly in phenylpropane and sugar metabolism. The mutant was mapped to a 2.72 Mb region on A01 by using BSA-Seq, and the region was validated by SSR markers.

CONCLUSIONS

The mutant chlorophyll content and photosynthetic efficiency were significantly reduced compared with those of the wild type. Abnormal chloroplasts and thylakoids less connected to the stroma lamella appeared in the mutant. This work on the mutant will facilitate the process of cloning the BnaA01.cd gene and provide more genetic and physiological information concerning chloroplast development in Brassica napus.

Genome-wide association mapping and genomic prediction for kernel color traits in intermediate wheatgrass (Thinopyrum intermedium)

BACKGROUND

Intermediate wheatgrass (IWG) is a novel perennial grain crop currently undergoing domestication. It offers important ecosystem benefits while producing grain suitable for human consumption. Several aspects of plant biology and genetic control are yet to be studied in this new crop. To understand trait behavior and genetic characterization of kernel color in IWG breeding germplasm from the University of Minnesota was evaluated for the CIELAB components (L*, a*, b*) and visual differences. Trait values were used in a genome-wide association scan to reveal genomic regions controlling IWG's kernel color. The usability of genomic prediction in predicting kernel color traits was also evaluated using a four-fold cross validation method.

RESULTS

A wide phenotypic variation was observed for all four kernel color traits with pairwise trait correlations ranging from - 0.85 to 0.27. Medium to high estimates of broad sense trait heritabilities were observed and ranged from 0.41 to 0.78. A genome-wide association scan with single SNP markers detected 20 significant marker-trait associations in 9 chromosomes and 23 associations in 10 chromosomes using multi-allelic haplotype blocks. Four of the 20 significant SNP markers and six of the 23 significant haplotype blocks were common between two or more traits. Evaluation of genomic prediction of kernel color traits revealed the visual score to have highest mean predictive ability (r² = 0.53); r² for the CIELAB traits ranged from 0.29-0.33. A search for candidate genes led to detection of seven IWG genes in strong alignment with MYB36 transcription factors from other cereal crops of the Triticeae tribe. Three of these seven IWG genes had moderate similarities with R-A1, R-B1, and R-D1, the three genes that control grain color in wheat.

CONCLUSIONS

We characterized the distribution of kernel color in IWG for the first time, which revealed a broad phenotypic diversity in an elite breeding germplasm. Identification of genetic loci controlling the trait and a proof-of-concept that genomic selection might be useful in selecting genotypes of interest could help accelerate the breeding of this novel crop towards specific end-use.

Missense mutation of a class B heat shock factor is responsible for the tomato bushy root-2 phenotype

The bushy root-2 (brt-2) tomato mutant has twisting roots, and slower plant development. Here we used whole genome resequencing and genetic mapping to show that brt-2 is caused by a serine to cysteine (S75C) substitution in the DNA binding domain (DBD) of a heat shock factor class B (HsfB) encoded by SolycHsfB4a. This gene is orthologous to the Arabidopsis SCHIZORIZA gene, also known as AtHsfB4. The brt-2 phenotype is very similar to Arabidopsis lines in which the function of AtHsfB4 is altered: a proliferation of lateral root cap and root meristematic tissues, and a tendency for lateral root cap cells to easily separate. The brt-2 S75C mutation is unusual because all other reported amino acid substitutions in the highly conserved DBD of eukaryotic heat shock factors are dominant negative mutations, but brt-2 is recessive. We further show through reciprocal grafting that brt-2 exerts its effects predominantly through the root genotype even through BRT-2 is expressed at similar levels in both root and shoot meristems. Since AtHsfB4 is induced by root knot nematodes (RKN), and loss-of-function mutants of this gene are resistant to RKNs, BRT-2 could be a target gene for RKN resistance, an important trait in tomato rootstock breeding.Gene & accession numbersSolycHsfB4a - Solyc04g078770.

Development and application of KASP markers associated with Restorer-of-fertility gene in Capsicum annuum L

Fertility restoration of cytoplasmic male sterility (CMS) in Capsicum annuum is controlled by multiple alleles of Restorer-of-fertility (Rf) genes. The isolation of additional Rf genes should therefore enrich the knowledge of CMS/Rf systems and accelerate their exploitation in hybrid seed production. In this study, the fertility restorer gene CaRfm of '0601 M', a non-pungent bell pepper, was genetically mapped to a 1.2-cM region flanked by KASP markers S761 and S183. CaRfm was then physically mapped to a 128.96-Kb interval predicted from 24 recombinants with two co-segregated markers, S423 and S424. CaPPR6 encoding a pentatricopeptide repeat (PPR) protein was suggested as the most likely candidate gene for the CaRfm locus on the basis of sequence alignment as well as genotyping of tightly linked markers. In addition, molecular markers S1597 and S1609, which are immediately adjacent to CaRfm at 15.7 and 57.8-Kb respectively, were developed and applied to marker-assisted selection. The results provided friendly markers for breeding pepper restorer lines and laid the foundation for elucidating the male fertility restoration mechanism.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01109-9.

Identification of QTLs for stagnant flooding tolerance in rice employing genotyping by sequencing of a RIL population derived from Swarna × Rashpanjor

In lowland rice ecosystems stagnant flooding or partial submergence has a significant negative impact on important yield attributing traits resulting in substantial grain yield reduction. Genetics of this stress is not yet studied intensively. Rashpanjor (IC 575321), a landrace from India, was identified and used as the tolerant donor for stagnant flooding and was crossed with high yielding variety Swarna to develop the RIL population for the present investigation. Yield and yield attributing traits of 180 F_2:8 lines in rainfed non-stressed and stressed (stagnant flooding with 45 ± 5 cm standing water) conditions were recorded in the wet season of 2018 and stress susceptibility and tolerance indices of yield component traits were deduced. Homo-polymorphic high-quality SNPs between two parents derived from genotyping by sequencing were employed and 17 putative QTLs for plant height, shoot elongation, panicle number, grain weight, panicle length in control and stagnant flooding conditions were identified. Tolerance and susceptibility indexes for these traits were detected in chromosomes 1, 3, 4, 5, 6, 10, 11, and 12 with PVE ranging from 6.53 to 57.89%. Two major QTLs clusters were found for stress susceptibility index of grain and panicle weight on chromosome 1 and plant height in non-stress condition and stress tolerance index of elongation ability on chromosome 3. Putative functional genes present either in associated non-synonymous SNPs or inside the QTL regions were also predicted. Some of them were directly associated with ethylene biosynthesis and encoding auxin responsive factors for better adaptation under stagnant flooding and also coded for different transcription factors viz. NAC domain-binding protein, WRKY gene family, and MYB class known for ROS scavenging and production of metabolites to enhance tolerance to stagnant flooding.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01107-x.

A proximal LAVA method for genome-wide association and prediction of traits with mixed inheritance patterns

BACKGROUND

The genetic basis of phenotypic traits is highly variable and usually divided into mono-, oligo- and polygenic inheritance classes. Relatively few traits are known to be monogenic or oligogeneic. The majority of traits are considered to have a polygenic background. To what extent there are mixtures between these classes is unknown. The rapid advancement of genomic techniques makes it possible to directly map large amounts of genomic markers (GWAS) and predict unknown phenotypes (GWP). Most of the multi-marker methods for GWAS and GWP falls into one of two regularization frameworks. The first framework is based on [Formula: see text]-norm regularization (e.g. the LASSO) and is suitable for mono- and oligogenic traits, whereas the second framework regularize with the [Formula: see text]-norm (e.g. ridge regression; RR) and thereby is favourable for polygenic traits. A general framework for mixed inheritance is lacking.

RESULTS

We have developed a proximal operator algorithm based on the recent LAVA regularization method that jointly performs [Formula: see text]- and [Formula: see text]-norm regularization. The algorithm is built on the alternating direction method of multipliers and proximal translation mapping (LAVA ADMM). When evaluated on the simulated QTLMAS2010 data, it is shown that the LAVA ADMM together with Bayesian optimization of the regularization parameters provides an efficient approach with lower test prediction mean-squared-error (65.89) than the LASSO (66.11), Ridge regression (83.41) and Elastic net (66.11). For the real pig data the test MSE of the LAVA ADMM is 0.850 compared to the LASSO, RR and EN with 0.875, 0.853 and 0.853, respectively.

CONCLUSIONS

This study presents the LAVA ADMM that is capable of joint modelling of monogenic major genetic effects and polygenic minor genetic effects which can be used for both genome-wide assoiciation and prediction purposes. The statistical evaluations based on both simulated and real pig data set shows that the LAVA ADMM has better prediction properies than the LASSO, RR and EN. Julia code for the LAVA ADMM is available at: https://github.com/patwa67/LAVAADMM .

Development of molecular markers lınked to QTL/genes controllıng Zn effıcıency

BACKGROUND

Zinc (Zn) deficiency is a widespread problem in reducing the yield and quality of crop plants worldwide. It is important to utilize molecular markers linked to Zn efficiency to develop high Zn efficient cultivars in pepper (Capsicum annuum L.).

METHODS AND RESULTS

In present study, genetic map was constructed using a F₂ populations derived from C. annuum L. (Alata 21A) × C. frutescens L. (PI 281420) cross. The quantitative trait locus (QTLs) for Zn efficiency were mapped using F_2:3 population. A genetic map with 929.6 cM in length and 12 linkage groups were obtained using 62 markers (31 sequence-related amplified polymorphism (SRAP), 19 simple sequence repeat (SSR) and 11 random amplified polymorphic DNA (RAPD) markers). The 41 linked QTLs related with nine (9) Zn efficiency characters were mapped on linkage groups. Results suggest that selecting plants for tolerance to Zn deficiency are expected to be rather responsive among segregating populations for breeding and developing Zn efficient genotypes in pepper.

CONCLUSIONS

The molecular markers are expected to aid selection and expedite breeding peppers resistant to Zn deficiency in soils low for available Zn contents.

Genetic mapping of host resistance to the Pyrenophora teres f. maculata isolate 13IM8.3

Spot form net blotch (SFNB), caused by the necrotrophic fungal pathogen Pyrenophora teres f. maculata (Ptm), is a foliar disease of barley that results in significant yield losses in major growing regions worldwide. Understanding the host-parasite interactions between pathogen virulence/avirulence genes and the corresponding host susceptibility/resistance genes is important for the deployment of genetic resistance against SFNB. Two recombinant inbred mapping populations were developed to characterize genetic resistance/susceptibility to the Ptm isolate 13IM8.3, which was collected from Idaho (ID). An Illumina Infinium array was used to produce a genome wide marker set. Quantitative trait loci (QTL) analysis identified ten significant resistance/susceptibility loci, with two of the QTL being common to both populations. One of the QTL on 5H appears to be novel, while the remaining loci have been reported previously. Single nucleotide polymorphisms (SNPs) closely linked to or delimiting the significant QTL have been converted to user-friendly markers. Loci and associated molecular markers identified in this study will be useful in genetic mapping and deployment of the genetic resistance to SFNB in barley.

QTL analysis of seed germination traits in tobacco (Nicotiana tabacum L.)

Genetic mapping of seed germination traits has been performed with many plant species. In tobacco, however, investigations are rare. In the present study, a bi-parental mapping population consisting of 118 doubled haploid lines and derived from a cross between 'Beinhart-1000' and 'Hicks' was investigated. Four germination-related traits, total germination (TG), normal germination (NG), time to reach 50% of total germination (T50), and the area under the curve after 200 h of germination (AUC) were considered by examining seeds either untreated or after a moderate controlled deterioration (CD). Quantitative trait loci were found for all traits distributed on 11 out of the 24 linkage groups. It was demonstrated that, as in many other species, germination-related traits are very complex and under polygenic control.

GBS-SNP and SSR based genetic mapping and QTL analysis for drought tolerance in upland cotton

A recombinant inbred line mapping population of intra-species upland cotton was generated from a cross between the drought-tolerant female parent (AS2) and the susceptible male parent (MCU13). A linkage map was constructed deploying 1,116 GBS-based SNPs and public domain-based 782 SSRs spanning a total genetic distance of 28,083.03 cM with an average chromosomal span length of 1,080.12 cM with inter-marker distance of 10.19 cM.A total of 19 quantitative trait loci (QTLs) were identified in nine chromosomes for field drought tolerance traits. Chromosomes 3 and 8 harbored important drought tolerant QTLs for chlorophyll stability index trait while for relative water content trait, three QTLs on chromosome 8 and one QTL each on chromosome 4, 12 were identified. One QTL on each chromosome 8, 5, and 7, and two QTLs on chromosome 15 linking to proline content were identified. For the nitrate reductase activity trait, two QTLs were identified on chromosome 3 and one on each chromosome 8, 13, and 26. To complement our QTL study, a meta-analysis was conducted along with the public domain database and resulted in a consensus map for chromosome 8. Under field drought stress, chromosome 8 harbored a drought tolerance QTL hotspot with two in-house QTLs for chlorophyll stability index (qCSI01, qCSI02) and three public domain QTLs (qLP.FDT_1, qLP.FDT_2, qCC.ST_3). Identified QTL hotspot on chromosome 8 could play a crucial role in exploring abiotic stress-associated genes/alleles for drought trait improvement.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01041-y.

How much epigenetics and quantitative trait loci (QTL) mapping tell us about parasitism maintenance and resistance/susceptibility to hosts

Interactions between the environment, parasites, vectors, and/or intermediate hosts are complex and involve several factors that define the success or failure of an infection. Among these interactions that can affect infections by a parasite, it is possible to highlight the genetic and epigenetic mechanisms in hosts and parasites. The interaction between genetics, epigenetics, infection, and the host's internal and external environment is decisive and dictates the outcome of a parasitic infection and the resistance, susceptibility, and transmission of this parasite. Epigenetic changes become important mediators in the regulation of gene expression, allowing the evasion of the parasite to immune host barriers, its transmission to new hosts, and the end of its development cycle. Epigenetics is a new frontier in the understanding of the interaction mechanisms between parasite and host that, along with information from the gene regions associated with complex phenotypic variations, the Quantitative Trait Loci, brings new possibilities to investigate more modern and efficient approaches to the treatment, control, and eradication of parasitic diseases. In this brief review, a general overview of the use of epigenetic information and mapping of Quantitative Trait Loci was summarized, both in genes of parasites and hosts, for understanding the mechanisms of resistance and/or susceptibility in parasitic relationships; also, the main search platforms were quantitatively compared, aiming to facilitate access data produced over a period of twenty years.

A genome-wide genetic linkage map and reference quality genome sequence for a new race in the wheat pathogen Pyrenophora tritici-repentis

Pyrenophora tritici-repentis is an ascomycete fungus that causes tan spot of wheat. The disease has a worldwide distribution and can cause significant yield and quality losses in wheat production. The fungal pathogen is homothallic in nature, which means it can undergo sexual reproduction by selfing to produce pseudothecia on wheat stubble for seasonal survival. Since homothallism precludes the development of bi-parental fungal populations, no genetic linkage map has been developed for P. tritici-repentis for mapping and map-based cloning of fungal virulence genes. In this work, we created two heterothallic strains by deleting one of the mating type genes in each of two parental isolates 86-124 (race 2) and AR CrossB10 (a new race) and developed a bi-parental fungal population between them. The draft genome sequences of the two parental isolates were aligned to the Pt-1C-BFP reference sequence to mine single nucleotide polymorphisms (SNPs). A total of 225 SNP markers were developed for genotyping the entire population. Additionally, 75 simple sequence repeat, and two gene markers were also developed and used in the genotyping. The resulting linkage map consisted of 13 linkage groups spanning 5,075.83 cM in genetic distance. Because the parental isolate AR CrossB10 is a new race and produces Ptr ToxC, it was sequenced using long-read sequencing platforms and de novo assembled into contigs. The majority of the contigs were further anchored into chromosomes with the aid of the linkage maps. The whole genome comparison of AR CrossB10 to the reference genome of M4 revealed a few chromosomal rearrangements. The genetic linkage map and the new AR CrossB10 genome sequence are valuable tools for gene cloning in P. tritici-repentis.

Genetic analysis of the barley variegation mutant, grandpa1.a

BACKGROUND

Providing the photosynthesis factory for plants, chloroplasts are critical for crop biomass and economic yield. However, chloroplast development is a complicated process, coordinated by the cross-communication between the nucleus and plastids, and the underlying biogenesis mechanism has not been fully revealed. Variegation mutants have provided ideal models to identify genes or factors involved in chloroplast development. Well-developed chloroplasts are present in the green tissue areas, while the white areas contain undifferentiated plastids that are deficient in chlorophyll. Unlike albino plants, variegation mutants survive to maturity and enable investigation into the signaling pathways underlying chloroplast biogenesis. The allelic variegated mutants in barley, grandpa 1 (gpa1), have long been identified but have not been genetically characterized.

RESULTS

We characterized and genetically analyzed the grandpa1.a (gpa1.a) mutant. The chloroplast ultrastructure was evaluated using transmission electron microscopy (TEM), and it was confirmed that chloroplast biogenesis was disrupted in the white sections of gpa1.a. To determine the precise position of Gpa1, a high-resolution genetic map was constructed. Segregating individuals were genotyped with the barley 50 k iSelect SNP Array, and the linked SNPs were converted to PCR-based markers for genetic mapping. The Gpa1 gene was mapped to chromosome 2H within a gene cluster functionally related to photosynthesis or chloroplast differentiation. In the variegated gpa1.a mutant, we identified a large deletion in this gene cluster that eliminates a putative plastid terminal oxidase (PTOX).

CONCLUSIONS

Here we characterized and genetically mapped the gpa1.a mutation causing a variegation phenotype in barley. The PTOX-encoding gene in the delimited region is a promising candidate for Gpa1. Therefore, the present study provides a foundation for the cloning of Gpa1, which will elevate our understanding of the molecular mechanisms underlying chloroplast biogenesis, particularly in monocot plants.

Mapping chromosomal regions associated with anther indehiscence with exerted stigmas in CRI-48 and Jasmine 85 cross of rice (Oryza sativa L)

Anther indehiscence in certain wide crosses combines male sterility with stigma exertion, a phenomenon that is desirable for hybrid rice seed production. This study sought to identify chromosomal region(s) that combine anther indehiscence with exerted stigmas. A mapping population consisting of 189 BC1F1 plants was derived from a cross between CRI-48 and Jasmine 85 and backcrossing the resulting F1 to Jasmine 85. Contrary to the three complementary genes mode of inheritance reported earlier, a single locus (AI6-1) was mapped on chromosome 6 at 27.4 cM for anther indehiscence with exerted stigmas through a mixed model-based composite interval mapping (MCIM). This locus was flanked by two single nucleotide polymorphism (SNP) markers, K_ID6002884 and K_ID6003341 within a range of 23.1-28.9 cM. The allele at the locus was contributed by the CRI-48 parent which has Oryza glaberrima ancestry. This locus is suggested to control anther indehiscence and stigma exertion through pleiotropic gene action or cluster of genes.

Differentiation in the genetic basis of stem trichome development between cultivated tetraploid cotton species

BACKGROUND

Cotton stem trichomes and seed fibers are each single celled structures formed by protrusions of epidermal cells, and were found sharing the overlapping molecular mechanism. Compared with fibers, cotton stem trichomes are more easily observed, but the molecular mechanisms underlying their development are still poorly understood.

RESULTS

In this study, Gossypium hirsutum (Gh) and G. barbadense (Gb) were found to differ greatly in percentages of varieties/accessions with glabrous stems and in trichome density, length, and number per trichopore. Gh varieties normally had long singular and clustered trichomes, while Gb varieties had short clustered trichomes. Genetic mapping using five F₂ populations from crosses between glabrous varieties and those with different types of stem trichomes revealed that much variation among stem trichome phenotypes could be accounted for by different combinations of genes/alleles on Chr. 06 and Chr. 24. The twenty- six F₁ generations from crosses between varieties with different types of trichomes had varied phenotypes, further suggesting that the trichomes of tetraploid cotton were controlled by different genes/alleles. Compared to modern varieties, a greater proportion of Gh wild accessions were glabrous or had shorter and denser trichomes; whereas a smaller proportion of Gb primitive accessions had glabrous stems. A close correlation between fuzz fiber number and stem trichome density was observed in both Gh and Gb primitive accessions and modern varieties.

CONCLUSION

Based on these findings, we hypothesize that stem trichomes evolved in parallel with seed fibers during the domestication of cultivated tetraploid cotton. In addition, the current results illustrated that stem trichome can be used as a morphological index of fiber quality in cotton conventional breeding.

Multiparental Population in Crops: Methods of Development and Dissection of Genetic Traits

Multiparental populations are located midway between association mapping that relies on germplasm collections and classic linkage analysis, based upon biparental populations. They provide several key advantages such as the possibility to include a higher number of alleles and increased level of recombination with respect to biparental populations, and more equilibrated allelic frequencies than association mapping panels. Moreover, in these populations new allele's combinations arise from recombination that may reveal transgressive phenotypes and make them a useful pre-breeding material. Here we describe the strategies for working with multiparental populations, focusing on nested association mapping populations (NAM) and multiparent advanced generation intercross populations (MAGIC). We provide details from the selection of founders, population development, and characterization to the statistical methods for genetic mapping and quantitative trait detection.

Identification of Solanum Immune Receptors by Bulked Segregant RNA-Seq and High-Throughput Recombinant Screening

The identification, understanding, and deployment of immune receptors are crucial to achieve high-level and durable resistance for crops against pathogens. In potato, many R genes have been identified using map-based cloning strategies. However, this is a challenging and laborious task that involves the development of a high number of molecular markers for the initial mapping, and the screening of thousands of plants for fine mapping. Bulked segregant RNA-Seq (BSR-Seq) has proven to be an efficient technique for the mapping of resistance genes. The RNA from two bulks of plants with contrasting phenotypes is sequenced and analyzed to identify single-nucleotide polymorphism (SNPs) markers linked to the target gene. Subsequently, the SNP markers that are identified can be used to delimit the mapping interval. Additionally, we designed an in vitro recombinant screening strategy that is advantageous for analyzing a large number of plants, in terms of time, space, and cost. Tips and detailed protocols, including BSR-Seq, bioinformatic analysis, and recombinant screening, are provided in this chapter.

Genetic and Physical Localization of the Gene Controlling Leaf Pigmentation Pattern in Medicago truncatula

In Medicago truncatula, some ecotypes form a black or purple stain in the middle of adaxial leaf surface due to accumulation of anthocyanins. However, this morphological marker is missing in some other ecotypes, although anthocyanin biosynthesis pathway is not disrupted. Genetic analysis indicated that the lack of the leaf spot of anthocyanins accumulation is a dominant trait, which is controlled by a single gene, LPP1. Genetic mapping indicated that the LPP1 gene was delimited to a 280 kb-region on Chromosome 7. A total of 8 protein-coding genes were identified in the LPP1 locus through gene annotation and sequence analysis. Of those, two genes, putatively encoding MYB-transcriptional suppressors, were selected as candidates for functional validation.

Genetic Mapping in Autohexaploid Sweet Potato with Low-Coverage NGS-Based Genotyping Data

Next-generation sequencing (NGS)-based genotyping methods can generate numerous genetic markers in a single experiment and have contributed to plant genetic mapping. However, for high precision genetic analysis, the complicated genetic segregation mode in polyploid organisms requires high-coverage NGS data and elaborate analytical algorithms. In the present study, we propose a simple strategy for the genetic mapping of polyploids using low-coverage NGS data. The validity of the strategy was investigated using simulated data. Previous studies indicated that accurate allele dosage estimation from low-coverage NGS data (read depth < 40) is difficult. Therefore, we used allele dosage probabilities calculated from read counts in association analyses to detect loci associated with phenotypic variations. The allele dosage probabilities showed significant detection power, although higher allele dosage estimation accuracy resulted in higher detection power. On the contrary, differences in the segregation patterns between the marker and causal genes resulted in a drastic decrease in detection power even if the marker and casual genes were in complete linkage and the allele dosage estimation was accurate. These results indicated that the use of a larger number of markers is advantageous, even if the accuracy of allele dosage estimation is low. Finally, we applied the strategy for the genetic mapping of autohexaploid sweet potato (Ipomoea batatas) populations to detect loci associated with agronomic traits. Our strategy could constitute a cost-effective approach for preliminary experiments done performed to large-scale studies.

Construction of a high-density genetic linkage map and identification of gene controlling resistance to cucumber mosaic virus in Luffa cylindrica (L.) Roem. based on specific length amplified fragment sequencing

Luffa cylindrica L. is a cash crop which has important health, medicinal and industrial value, but no high saturation genetic map has been constructed owing to a lack of efficient markers. Furthermore, no genes were reportedly responsible for CMV resistance in Luffa spp. Specific length amplified fragment sequencing (SLAF-seq) is a valuable tool for large-scale discovery of markers and genetic mapping. The present study reported the construction of a high-density genetic map and the mapping of CMV resistant genes by using an F₂ population of 130 individuals and their two inbred line parents. A total of 271.01 Mb pair-end reads were generated. 100,077 high-quality SLAFs were detected, and 7404 of them were polymorphic. Finally, 3701 of the polymorphic markers were selected for genetic map construction, and 13 linkage groups were generated. The map spanned 1518.56 cM with an average distance of 0.41 cM between adjacent markers. Based on the newly constructed high-density map, one gene located on chromosome 1 (100.072-100.457 cM) was identified to regulate CMV resistance in L. cylindrica. A gag-polypeptide of LTR copia-type retrotransposon was predicted as the candidate gene responsible for CMV resistance in L. cylindrica. The high-density genetic map and the CMV resistant gene mapped and predicted in this study will be useful in future research.

Identification of the vernalization gene VRN-B1 responsible for heading date variation by QTL mapping using a RIL population in wheat

BACKGROUND

Heading time is one of the most important agronomic traits in wheat, as it largely affects both adaptation to different agro-ecological conditions and yield potential. Identification of genes underlying the regulation of wheat heading and the development of diagnostic markers could facilitate our understanding of genetic control of this process.

RESULTS

In this study, we developed 400 recombinant inbred lines (RILs) by crossing a γ-ray-induced early heading mutant (eh1) with the late heading cultivar, Lunxuan987. Bulked Segregant Analysis (BSA) of both RNA and DNA pools consisting of various RILs detected a quantitative trait loci (QTL) for heading date located on chromosomes 5B, and further genetic linkage analysis limited the QTL to a 3.31 cM region. We then identified a large deletion in the first intron of the vernalization gene VRN-B1 in eh1, and showed it was associated with the heading phenotype in the RIL population. However, it is not the mutation loci that resulted in early heading phonotype in the mutant compared to that of wildtype. RNA-seq analysis suggested that Vrn-B3 and several newly discovered genes, including beta-amylase 1 (BMY1) and anther-specific protein (RTS), were highly expressed in both the mutant and early heading pool with the dominant Vrn-B1 genotype compared to that of Lunxuan987 and late heading pool. Enrichment analysis of differentially expressed genes (DEGs) identified several key pathways previously reported to be associated with flowering, including fatty acid elongation, starch and sucrose metabolism, and flavonoid biosynthesis.

CONCLUSION

The development of new markers for Vrn-B1 in this study supplies an alternative solution for marker-assisted breeding to optimize heading time in wheat and the DEGs analysis provides basic information for VRN-B1 regulation study.

Genetic localization of the SPC gene controlling pod coiling direction in Medicago truncatula

BACKGROUND

Handedness in plants introduced by helical growth of organs is frequently observed, and it has fascinated plant scientists for decades. However, the genetic control of natural handedness has not been revealed. In the model legume Medicago truncatula, pods can be coiled in a clockwise or anti-clockwise manner, providing a model for genetic analysis of plant handedness.

OBJECTIVE

We aimed to localize the Sense of Pod Coiling (SPC) gene controlling pod coiling direction in M. truncatula.

METHODS

Linkage analysis was used with a biparental population for fine mapping of the SPC gene. The genome sequence of M. truncatula Mt4.0 was used for marker identification and physical mapping. Single nucleotide polymorphisms (SNPs) between the parental lines were converted to CAPS (cleaved amplified polymorphic sequences) markers. Genetic map was constructed using the software JoinMap version 3.0. Gene predication and annotation provided by the M. truncatula genome database (http://www.medicagogenome.org) was confirmed with the programs of FGENESH and Pfam 32.0, respectively. Quantitative reverse transcription PCR (qRT-PCR) was used to analyze the relative expression levels of candidate genes.

RESULTS

The genetic analysis indicated that the anti-clockwise coiling is dominant to clockwise and is controlled by the single gene, SPC. The SPC gene was delimited to a 250 kb-region on Chromosome 7. Total of 15 protein-coding genes were identified in the SPC locus through gene annotation and sequence analysis. Of those, two genes, potentially encoding a receptor-like kinase and a vacuolar cation/proton exchanger respectively, were selected as candidates for the SPC gene.

CONCLUSIONS

The result presented here lay a foundation for gene cloning of SPC, which will help us to understand the molecular mechanisms underlying helical growth in plant organs.

Identification of two recessive etiolation genes (py1, py2) in pakchoi (Brassica rapa L. ssp. chinensis)

BACKGROUND

Leaf color is a major agronomic trait, which has a strong influence on crop yields. Isolating leaf color mutants can represent valuable materials for research in chlorophyll (Chl) biosynthesis and metabolism regulation.

RESULTS

In this study, we identified a stably inherited yellow leaf mutant derived from 'Huaguan' pakchoi variety via isolated microspore culture and designated as pylm. This mutant displayed yellow leaves after germination. Its etiolated phenotype was nonlethal and stable during the whole growth period. Its growth was weak and its hypocotyls were markedly elongated. Genetic analysis revealed that two recessive nuclear genes, named py1 and py2, are responsible for the etiolation phenotype. Bulked segregant RNA sequencing (BSR-Seq) showed that py1 and py2 were mapped on chromosomes A09 and A07, respectively. The genes were single Mendelian factors in F_3:4 populations based on a 3:1 phenotypic segregation ratio. The py1 was localized to a 258.3-kb interval on a 34-gene genome. The differentially expressed gene BraA09004189 was detected in the py1 mapping region and regulated heme catabolism. One single-nucleotide polymorphism (SNP) of BraA09004189 occurred in pylm. A candidate gene-specific SNP marker in 1520 F_3:4 yellow-colored individuals co-segregated with py1. For py2, 1860 recessive homozygous F_3:4 individuals were investigated and localized py2 to a 4.4-kb interval. Of the five genes in this region, BraA07001774 was predicted as a candidate for py2. It encoded an embryo defective 1187 and a phosphotransferase related to chlorophyll deficiency and hypocotyl elongation. One SNP of BraA07001774 occurred in pylm. It caused a single amino acid mutation from Asp to Asn. According to quantitative real-time polymerase chain reaction (qRT-PCR), BraA07001774 was downregulated in pylm.

CONCLUSIONS

Our study identified a Chl deficiency mutant pylm in pakchoi. Two recessive nuclear genes named py1 and py2 had a significant effect on etiolation. Candidate genes regulating etiolation were identified as BraA09004189 and BraA07001774, respectively. These findings will elucidate chlorophyll metabolism and the molecular mechanisms of the gene interactions controlling pakchoi etiolation.

Using Heterogeneous Stocks for Fine-Mapping Genetically Complex Traits

In this chapter we will review both the rationale and experimental design for using Heterogeneous Stock (HS) populations for fine-mapping of complex traits in mice and rats. We define an HS as an outbred population derived from an intercross between two or more inbred strains. HS have been used to perform genome-wide association studies (GWAS) for multiple behavioral, physiological, and gene expression traits. GWAS using HS require four key steps, which we review: selection of an appropriate HS population, phenotyping, genotyping, and statistical analysis. We provide advice on the selection of an HS, comment on key issues related to phenotyping, discuss genotyping methods relevant to these populations, and describe statistical genetic analyses that are applicable to genetic analyses that use HS.

Enhancing Upland cotton for drought resilience, productivity, and fiber quality: comparative evaluation and genetic dissection

To provision the world sustainably, modern society must increase overall crop production, while conserving and preserving natural resources. Producing more with diminishing water resources is an especially daunting endeavor. Toward the goal of genetically improving drought resilience of cultivated Upland cotton (Gossypium hirsutum L.), this study addresses the genetics of differential yield components referred to as productivity and fiber quality traits under regular-water versus low-water (LW) field conditions. We used ten traits to assess water stress deficit, which included six productivity and four fiber quality traits on two recombinant inbred line (RIL) populations from reciprocally crossed cultivars, Phytogen 72 and Stoneville 474. To facilitate genetic inferences, we genotyped RILs with the CottonSNP63K array, assembled high-density linkage maps of over 7000 SNPs and then analyzed quantitative trait variations. Analysis of variance revealed significant differences for all traits (p < 0.05) in these RIL populations. Although the LW irrigation regime significantly reduced all traits, except lint percent, the RILs exhibited a broad phenotypic spectrum of heritable differences across the water regimes. Transgressive segregation occurred among the RILs, suggesting the possibility of genetic gain through phenotypic selection for drought resilience and perhaps through marker-based selection. Analyses revealed more than 150 quantitative trait loci (QTLs) associated with productivity and fiber quality traits (p < 0.005) on different genomic regions of the cotton genome. The multiple-QTL models analysis with LOD > 3.0 detected 21 QTLs associated with productivity and 22 QTLs associated with fiber quality. For fiber traits, strong clustering and QTL associations occurred in c08 and its homolog c24 as well as c10, c14, and c21. Using contemporary genome sequence assemblies and bioinformatically related information, the identification of genomic regions associated with responses to plant stress/drought elevates the possibility of using marker-assisted and omics-based selection to enhance breeding for drought resilient cultivars and identifying candidate genes and networks. RILs with different responses to drought indicated that it is possible to maintain high fiber quality under LW conditions or reduce the of LW impact on quality. The heritable variation among elite bi-parental RILs for productivity and quality under field drought conditions, and their association of QTLs, and thus specific genomic regions, indicate opportunities for breeding-based gains in water resource conservation, i.e., enhancing cotton's agricultural sustainability.

Preferential insertion of a Ty1 LTR-retrotransposon into the A sub-genome's HD1 gene significantly correlated with the reduction in stem trichomes of tetraploid cotton

Stem trichomes and seed fibers originate from epidermal cells and partially share a regulatory pathway at the molecular level. In Gossypium barbadense, two insertions of a Ty1 long-terminal repeat-retrotransposon [transposable element TE1 and TE2] in a homeodomain-leucine zipper gene (HD1) result in glabrous stems. The primers used to identify the TE insertions in G. barbadense were applied to screen for the same events in 81 modern G. hirsutum varieties and 31 wild races. Three wild races were found carrying the same TEs as G. barbadense. However, the TE insertions in two of these wild races occurred at different sites (4th exon), therefore, named TE3, while the TE in the other wild race occurred at the same site as TE2. An RNA sequencing and qRT-PCR analysis indicated that the loss of HD1 function was caused by the TE insertion. Genetic mapping revealed a strong association between glabrous stems and TE3 insertions, confirming that HD1 is a critical gene for stem trichome initiation in G. hirsutum, as in G. barbadense. Using the long-terminal repeat sequence as a query to search against the Texas Marker-1 reference genome sequence, we found that the TE occurred after tetraploid cotton formation and evolved at different rates in G. hirsutum and G. barbadense. Interestingly, at least three independent insertion events of the same retrotransposon occurred preferentially in the A sub-genome's HD1 gene, but not in the D sub-genome of G. hirsutum or G. barbadense, suggesting that an unknown TE insertion mechanism and resultant gene function changes may have hastened cotton speciation.

A chromosome-scale genome assembly and dense genetic map for Xenopus tropicalis

The Western clawed frog Xenopus tropicalis is a diploid model system for both frog genetics and developmental biology, complementary to the paleotetraploid X. laevis. Here we report a chromosome-scale assembly of the X. tropicalis genome, improving the previously published draft genome assembly through the use of new assembly algorithms, additional sequence data, and the addition of a dense genetic map. The improved genome enables the mapping of specific traits (e.g., the sex locus or Mendelian mutants) and the characterization of chromosome-scale synteny with other tetrapods. We also report an improved annotation of the genome that integrates deep transcriptome sequence from diverse tissues and stages. The exon-intron structures of these genes are highly conserved relative to both X. laevis and human, as are chromosomal linkages ("synteny") and local gene order. A network analysis of developmental gene expression will aid future studies.

Mapping and characterization of a quantitative trait locus resistance to the brown planthopper in the rice variety IR64

Background

Rice planthoppers (main brown planthopper, Nilaparvata lugens Stål; BPH) was one of substantial threats to Asia rice production as its serious destruction and difficulties in control under field conditions. Notably, host-plant resistance was proved to be one of the effective ways to manage the pest. And stronger virulence will probably emergence when continuous use of insecticides. Therefore, more resistance genes with different resistance mechanisms were needed to be detected and then applied in the rice breeding practice.

Results

Resistance genes in the rice variety IR64 were evaluated considering the seedling bulk test and seedling survival rate. As a result, a locus with a large LOD score of 7.23 was found between markers RM302 and YM35 on chromosome 1. The locus explained 36.9% of phenotypic variation and was tentatively denominated Bph37. Moreover, Bph1 was detected to be harbored by the markers RM28366 and RM463, and had the largest LOD score of 2.08, explaining 7.7% of phenotypic variance in the same mapping population. Finally, the preliminary-near-isogenic-lines (pre-NILs) carrying Bph37 exhibited significant tolerance to the insects. But no antibiotic or antixenotic effects were observed in the resistant plants when infested with the insects.

Conclusions

We mapped one major BPH resistance gene Bph37 in consideration of seedling survival rate and the resistance lines showed tolerance to BPH. The detected gene should be beneficial for understanding the resistance mechanism of rice to BPH and for insect-resistance rice breeding programs.

Development and characterisation of interspecific hybrid lines with genome-wide introgressions from Triticum timopheevii in a hexaploid wheat background

BACKGROUND

Triticum timopheevii (2n = 4x = 28; A^tA^tGG), is an important source for new genetic variation for wheat improvement with genes for potential disease resistance and salt tolerance. By generating a range of interspecific hybrid lines, T. timopheevii can contribute to wheat's narrow gene-pool and be practically utilised in wheat breeding programmes. Previous studies that have generated such introgression lines between wheat and its wild relatives have been unable to use high-throughput methods to detect the presence of wild relative segments in such lines.

RESULTS

A whole genome introgression approach, exploiting homoeologous recombination in the absence of the Ph1 locus, has resulted in the transfer of different chromosome segments from both the A^t and G genomes of T. timopheevii into wheat. These introgressions have been detected and characterised using single nucleotide polymorphism (SNP) markers present on a high-throughput Axiom® Genotyping Array. The analysis of these interspecific hybrid lines has resulted in the detection of 276 putative unique introgressions from T. timopheevii, thereby allowing the generation of a genetic map of T. timopheevii containing 1582 SNP markers, spread across 14 linkage groups representing each of the seven chromosomes of the A^t and G genomes of T. timopheevii. The genotyping of the hybrid lines was validated through fluorescence in situ hybridisation (FISH). Comparative analysis of the genetic map of T. timopheevii and the physical map of the hexaploid wheat genome showed that synteny between the two species is highly conserved at the macro-level and confirmed the presence of inter- and intra-genomic translocations within the A^t and G genomes of T. timopheevii that have been previously only detected through cytological techniques.

CONCLUSIONS

In this work, we report a set of SNP markers present on a high-throughput genotyping array, able to detect the presence of T. timopheevii in a hexaploid wheat background making it a potentially valuable tool for marker assisted selection (MAS) in wheat pre-breeding programs. These valuable resources of high-density molecular markers and wheat-T. timopheevii hybrid lines will greatly enhance the work being undertaken for wheat improvement through wild relative introgressions.

Genetic analysis of drought response of wheat following either chemical desiccation or the use of a rain-out shelter

Simulating drought stress during the breeding process has been proposed as a way to select varieties under naturally non-stressful conditions. The aim of the study was to characterise the genetic basis of the response of 111 spring wheat (Triticum aestivum L.) varieties and landraces to chemical desiccation and to rain-out shelter drought. The effect of the rain-out shelter was a 15% reduction in plant height, spike length and thousand seed weight (TSW); in contrast, the desiccant treatment induced a 15% reduction in seed number, a 35-72% loss in TSW and a reduction in subsequent germination of 12%. A genome-wide association analysis revealed 263 significant marker-trait associations (MTAs), of which 246 involved days to anthesis, plant height, spike length, number of spikelets, seed number, TSW and germination from the non-treated plants. Only four and five MTAs involved TSW from plants grown under the rain-out shelter and the chemical desiccation, respectively, and harboured the Sugar-Dependent6 gene. Seven MTAs involved seed number for chemical desiccated plants. Both, chemical desiccation and rain-out shelter drought identified same tolerant genotypes. Concluding, both approaches are suitable to simulate different drought scenarios. However, there was a strong environmental impact for chemical desiccation which may increase the complexity of this tolerance mechanism.

Population genomic and evolutionary modelling analyses reveal a single major QTL for ivermectin drug resistance in the pathogenic nematode, Haemonchus contortus

BACKGROUND

Infections with helminths cause an enormous disease burden in billions of animals and plants worldwide. Large scale use of anthelmintics has driven the evolution of resistance in a number of species that infect livestock and companion animals, and there are growing concerns regarding the reduced efficacy in some human-infective helminths. Understanding the mechanisms by which resistance evolves is the focus of increasing interest; robust genetic analysis of helminths is challenging, and although many candidate genes have been proposed, the genetic basis of resistance remains poorly resolved.

RESULTS

Here, we present a genome-wide analysis of two genetic crosses between ivermectin resistant and sensitive isolates of the parasitic nematode Haemonchus contortus, an economically important gastrointestinal parasite of small ruminants and a model for anthelmintic research. Whole genome sequencing of parental populations, and key stages throughout the crosses, identified extensive genomic diversity that differentiates populations, but after backcrossing and selection, a single genomic quantitative trait locus (QTL) localised on chromosome V was revealed to be associated with ivermectin resistance. This QTL was common between the two geographically and genetically divergent resistant populations and did not include any leading candidate genes, suggestive of a previously uncharacterised mechanism and/or driver of resistance. Despite limited resolution due to low recombination in this region, population genetic analyses and novel evolutionary models supported strong selection at this QTL, driven by at least partial dominance of the resistant allele, and that large resistance-associated haplotype blocks were enriched in response to selection.

CONCLUSIONS

We have described the genetic architecture and mode of ivermectin selection, revealing a major genomic locus associated with ivermectin resistance, the most conclusive evidence to date in any parasitic nematode. This study highlights a novel genome-wide approach to the analysis of a genetic cross in non-model organisms with extreme genetic diversity, and the importance of a high-quality reference genome in interpreting the signals of selection so identified.