Genome- and Transcriptome-wide Association Studies to Discover Candidate Genes for Diverse Root Phenotypes in Cultivated Rice.

Root system architecture plays a crucial role in nutrient and water absorption during rice production. Genetic improvement of the rice root system requires elucidating its genetic control. Genome-wide association studies (GWASs) have identified genomic regions responsible for rice root phenotypes. However, candidate gene prioritization around the peak region often suffers from low statistical power and resolution. Transcriptomics enables other statistical mappings, such as transcriptome-wide association study (TWAS) and expression GWAS (eGWAS), which improve candidate gene identification by leveraging the natural variation of the expression profiles. To explore the genes responsible for root phenotypes, we conducted GWAS, TWAS, and eGWAS for 12 root phenotypes in 57 rice accessions using 427,751 single nucleotide polymorphisms (SNPs) and the expression profiles of 16,901 genes expressed in the roots. The GWAS identified three significant peaks, of which the most significant peak responsible for seven root phenotypes (crown root length, crown root surface area, number of crown root tips, lateral root length, lateral root surface area, lateral root volume, and number of lateral root tips) was detected at 6,199,732 bp on chromosome 8. In the most significant GWAS peak region, OsENT1 was prioritized as the most plausible candidate gene because its expression profile was strongly negatively correlated with the seven root phenotypes. In addition to OsENT1, OsEXPA31, OsSPL14, OsDEP1, and OsDEC1 were identified as candidate genes responsible for root phenotypes using TWAS. Furthermore, a cis-eGWAS peak SNP was detected for OsDjA6, which showed the eighth strongest association with lateral root volume in the TWAS. The cis-eGWAS peak SNP for OsDjA6 was in strong linkage disequilibrium (LD) with a GWAS peak SNP on the same chromosome for lateral root volume and in perfect LD with another SNP variant in a putative cis-element at the 518 bp upstream of the gene. These candidate genes provide new insights into the molecular breeding of root system architecture.

Independent control of organ number and distribution pattern in rice panicle.

As the determinants of yield products, rice panicle traits are important targets for breeding. Despite their importance in grain filling and subsequent yield productivity, knowledge on the organ distribution pattern in rice panicles is limited owing to the lack of objective evaluation methods. In this study, we developed a method for quantifying rice panicle organ distribution patterns. To validate our method for practical application in biology, we integrated this method into a quantitative trait locus (QTL) analysis and identified QTLs for panicle organ distribution patterns in rice. Interestingly, Grain number 1 (Gn1), a major QTL of organ number, was not identified as a QTL for distribution pattern, indicating that the number and distribution of panicle organs are independently controlled. This study provides insight into rice panicle organ distribution patterns that will help improve breeding targeting rice panicle architecture.

Genetic basis of maturity time is independent from that of flowering time and contributes to ecotype differentiation in common buckwheat (Fagopyrum esculentum Moench).

BACKGROUND: Common buckwheat is considered a quantitative short-day plant and is classified into the autumn (highly photoperiod sensitive), summer (weakly photoperiod sensitive), and intermediate ecotype. Understanding ecotype differentiation is essential for adaptive expansion and maximizing yield. The genetic analysis for ecotype has focused on photoperiod-dependent flowering time, whereas post-flowering traits such as seed set and maturity time might also regulate ecotype differentiation. RESULTS: A field experiment revealed that ecotype differentiation is mainly defined by the timing of seed set and maturation, whereas flowering time is less relevant. Thus, we focused on maturity time as a trait that defines the ecotype. To detect QTLs for maturity time, we developed two F2 populations derived from early × late-maturing accessions and intermediate × late-maturing accessions. Using genotyping by random amplicon sequencing-direct analysis, we generated a high-density linkage map. QTL analysis detected two major QTLs for maturity time, one in each F2 population. We also detected QTLs for flowering time at loci different from maturity time QTLs, which suggests that different genetic mechanisms regulate flowering and maturity. Association analysis showed that both QTLs for maturity time were significantly associated with variations in the trait across years. CONCLUSIONS: Maturity time appeared to be more suitable for explaining ecotype differentiation than flowering time, and different genetic mechanisms would regulate the timing of flowering and maturation. The QTLs and QTL-linked markers for maturity time detected here may be useful to extend the cultivation area and to fine-tune the growth period to maximize yield in buckwheat.

Improving Quantitative Traits in Self-Pollinated Crops Using Simulation-Based Selection With Minimal Crossing.

Genomic selection and marker-assisted recurrent selection have been applied to improve quantitative traits in many cross-pollinated crops. However, such selection is not feasible in self-pollinated crops owing to laborious crossing procedures. In this study, we developed a simulation-based selection strategy that makes use of a trait prediction model based on genomic information to predict the phenotype of the progeny for all possible crossing combinations. These predictions are then used to select the best cross combinations for the selection of the given trait. In our simulated experiment, using a biparental initial population with a heritability set to 0.3, 0.6, or 1.0 and the number of quantitative trait loci set to 30 or 100, the genetic gain of the proposed strategy was higher or equal to that of conventional recurrent selection method in the early selection cycles, although the number of cross combinations of the proposed strategy was considerably reduced in each cycle. Moreover, this strategy was demonstrated to increase or decrease seed protein content in soybean recombinant inbred lines using SNP markers. Information on 29 genomic regions associated with seed protein content was used to construct the prediction model and conduct simulation. After two selection cycles, the selected progeny had significantly higher or lower seed protein contents than those from the initial population. These results suggest that our strategy is effective in obtaining superior progeny over a short period with minimal crossing and has the potential to efficiently improve the target quantitative traits in self-pollinated crops.

Analysis of the Segregation Distortion of FcRAN1 Genotypes Based on Whole-Genome Resequencing of Fig (Ficus carica L.) Breeding Parents.

The common fig (Ficus carica L.) has a gynodioecious breeding system, and its sex phenotype is an important trait for breeding because only female plant fruits are edible. During breeding to select for female plants, we analyzed the FcRAN1 genotype, which is strongly associated with the sex phenotype. In 12 F1 populations derived from 13 cross combinations, the FcRAN1 genotype segregation ratio was 1:1, whereas the M119-226 × H238-107 hybridization resulted in an extremely male-biased segregation ratio (178:7 = male:female). This finding suggests that the segregation distortion was caused by some genetic factor(s). A whole-genome resequencing of breeding parents (paternal and maternal lines) identified 9,061 high-impact SNPs in the parents. A genome-wide linkage analysis exploring the gene(s) responsible for the distortion revealed 194 high-impact SNPs specific to Caprifig6085 (i.e., seed parent ancestor) and 215 high-impact SNPs specific to H238-107 (i.e., pollen parent) in 201 annotated genes. A comparison between the annotated genes and the genes required for normal embryo or gametophyte development and function identified several candidate genes possibly responsible for the segregation distortion. This is the first report describing segregation distortion in F. carica.

Choosing the optimal population for a genome-wide association study: A simulation of whole-genome sequences from rice.

A genome-wide association study (GWAS) needs to have a suitable population. The factors that affect a GWAS (e.g. population structure, sample size, and sequence analysis and field testing costs) need to be considered. Mixed populations containing subpopulations of different genetic backgrounds may be suitable populations. We conducted simulation experiments to see if a population with high genetic diversity, such as a diversity panel, should be added to a target population, especially when the target population harbors small genetic diversity. The target population was 112 accessions of Oryza sativa L. subsp. japonica, mainly developed in Japan. We combined the target population with three populations that had higher genetic diversity. These were 100 indica accessions, 100 japonica accessions, and 100 accessions with various genetic backgrounds. The results showed that the GWAS's power with a mixed population was generally higher than with a separate population. Also, the optimal GWAS populations varied depending on the fixation index (FST ) of the quantitative trait nucleotides (QTNs) and the polymorphism of QTNs in each population. When a QTN was polymorphic in a target population, a target population combined with a higher diversity population improved the QTN's detection power. By investigating FST and the expected heterozygosity (He ) as factors influencing the detection power, we showed that single nucleotide polymorphisms with high FST or low He are less likely to be detected by GWAS with mixed populations. Sequenced or genotyped germplasm collections can improve the GWAS's detection power by using a subset of the collections with a target population.

Coupling day length data and genomic prediction tools for predicting time-related traits under complex scenarios.

Genomic selection (GS) has proven to be an efficient tool for predicting crop-rank performance of untested genotypes; however, when the traits have intermediate optima (phenology stages), this implementation might not be the most convenient. GS might deliver high-rank correlations but incurring in serious bias. Days to heading (DTH) is a crucial development stage in rice for regional adaptability with a significant impact on yield potential. The objective of this research consisted in develop a novel method that accurately predicts time-related traits such as DTH in unobserved environments. For this, we propose an implementation that incorporates day length information (DL) in the prediction process for two relevant scenarios: CV0, predicting tested genotypes in unobserved environments (C method); and CV00, predicting untested genotypes in unobserved environments (CB method). The use of DL has advantages over weather data since it can be determined in advance just by knowing the location and planting date. The proposed methods showed that DL information significantly helps to improve the predictive ability of DTH in unobserved environments. Under CV0, the C method returned a root-mean-square error (RMSE) of 3.9 days, a Pearson correlation (PC) of 0.98 and the differences between the predicted and observed environmental means (EMD) ranged between -4.95 and 4.67 days. For CV00, the CB method returned an RMSE of 7.3 days, a PC of 0.93 and the EMD ranged between -6.4 and 4.1 days while the conventional GS implementation produced an RMSE of 18.1 days, a PC of 0.41 and the EMD ranged between -31.5 and 28.7 days.

IonBreeders: bioinformatics plugins toward genomics-assisted breeding.

Polymorphism information generated by next-generation sequencing (NGS) technologies has enabled applications of genome-wide markers assisted breeding. However, handling such large-scale data remains a challenge for experimental researchers and breeders, calling for the urgent development of a flexible and straightforward analysis tool for NGS data. We developed "IonBreeders" as bioinformatics plugins that implement general analysis steps from genotyping to genomic prediction. IonBreeders comprises three plugins, "ABH", "IMPUTATION", and "GENOMIC PREDICTION", for format conversion of genotyping data, preprocessing and imputation of genotyping data, and genomic prediction, respectively. "ABH" converts genotyping data derived from NGS into the ABH format, which is acceptable for our further plugins and with other breeding software tools, R/qtl, MapMaker, and AntMap. "IMPUTATION" filters out non-informative markers and imputes missing marker genotypes. In "GENOMIC PREDICTION", users can use four statistical methods based on their target trait, quantitative trait locus effect, and number of markers, and construct a prediction model for genomic selection. IonBreeders is operated in Torrent Suite, but can also handle genotype data in standard formats, e.g., Variant Call Format (VCF), by format conversion using free software or our provided scripts.

Genomics-assisted breeding in minor and pseudo-cereals.

Minor and pseudo-cereals, which can grow with lower input and often produce specific nutrients compared to major cereal crops, are attracting worldwide attention. Since these crops generally have a large genetic diversity in a breeding population, rapid genetic improvement can be possible by the application of genomics-assisted breeding methods. In this review, we discuss studies related to biparental quantitative trait locus (QTL) mapping, genome-wide association study, and genomic selection for minor and pseudo-cereals. Especially, we focus on the current progress in a pseudo-cereal, buckwheat. Prospects for the practical utilization of genomics-assisted breeding in minor and pseudo-cereals are discussed including the issues to overcome especially for these crops.

Simulation-based optimization of genomic selection scheme for accelerating genetic gain while preventing inbreeding depression in onion breeding.

Genomic selection (GS) is being increasingly employed in plant breeding programs to accelerate genetic gain of economically important traits. However, its efficiency differs greatly across species, due to differences in reproduction and breeding strategies. Onion (Allium cepa L.) is an out-crossing crop but can be easily self-pollinated. High inbreeding depression occurs, and contamination of self-pollinated seeds is unavoidable in onion breeding. Taking this into consideration, 10-year breeding programs with and without GS were simulated. In addition to general GS, we proposed GS schemes to prevent inbreeding depression by avoiding co-selection of close relatives and combining the shortening of generation time and updating of the prediction model. The results showed that general GS with shortening of generation time yielded the highest genetic gain among the selection schemes in early years. However, inbreeding increased rapidly, reaching very high levels in later years. The proposed GS combining shortening of generation time with updating of the prediction model was superior to the others in later years, as it yielded relatively high genetic gain while maintaining significantly low levels of inbreeding. These results suggested that GS can be beneficial in onion breeding, and an optimal scheme should be selected depending on the selection period.

Whole genome sequencing data of Oryza sativa L. cv. Takanari.

Oryza sativa L. 'Takanari' is one of the most productive indica cultivars [1,2]. Reciprocal chromosome segment substitution lines (CSSLs) derived from a cross between 'Koshihikari' and 'Takanari' are useful tools for the detection and precise mapping of target quantitative trait loci (QTL) in 'Takanari'. Although the available Os-Nipponbare-Reference-IRGSP-1.0 reference genome is available and useful for evaluating genetic diversity among japonica cultivars, it is not always useful for evaluating genetic diversity harbored by indica cultivars such as 'Takanari'. To reveal sequence variants in 'Takanari' and to exploit these variants in rice breeding programs, the whole genome of 'Takanari' was sequenced using a combination of Illumina HiSeq X Ten (20,983,495 reads and %GC 43) and PacBio (2,847,220 high-quality subreads). NGS data obtained have been deposited in the DNA Data Bank of Japan (DDBJ) under accession number DRA007557.

Description of grain weight distribution leading to genomic selection for grain-filling characteristics in rice.

Grain-filling ability is one of the factors that controls grain yield in rice (Oryza sativa L.). We developed a method for describing grain weight distribution, which is the probability density function of single grain weight in a panicle, using 128 Japanese rice varieties. With this method, we quantitively analyzed genotypic differences in grain-filling ability and used the grain weight distribution parameters for genomic prediction subject to genetic improvement in grain yield in rice. The novel description method could represent the observed grain weight distribution with five genotype-specific parameters of a mixture of two gamma distributions. The estimated genotype-specific parameters representing the proportion of filled grains had applicability to explain the grain filling ability of genotypes comparable to that of sink-filling rate and the conventionally measured proportion of filled grains, which suggested the efficiency and flexibility of grain weight distribution parameters to handle several genotypes. We revealed that perfectly filled grains have to be prioritized over partially filled grains for the optimum allocation of the source of yield in a panicle, from the analysis for obtaining an ideal shape of grain weight distribution. We conducted genomic prediction of grain weight distribution considering five genotype-specific parameters of the distribution as phenotypes relating to grain filling ability. The proportion of filled grains, average weight of filled grains, and variance of filled grain weight, which were considered to control grain yield to a certain degree, were predicted with accuracies of 0.30, 0.28, and 0.53, respectively. The proposed description method of grain weight distribution facilitated not only the investigation of the optimum allocation of nutrients in a panicle for realizing high grain-filling ability, but also allowed genomic selection of grain weight distribution.

Influence of parasitic chytrids on the quantity and quality of algal dissolved organic matter (AOM).

Algae-derived dissolved organic matter (AOM) is an important nutrient source for heterotrophic bacteria, while AOM such as humic substances pose significant challenges during water treatment processing. We hypothesized that the parasitic infection of algae could change the composition and concentration of AOM. This study investigated the quality and quantity of DOM and bacterial abundance in diatom (Synedra) cultures, with and without parasitic fungi (chytrids). The quality of DOM was analyzed using three-dimensional excitation-emission matrix combined with parallel factor analysis (EEM-PARAFAC) and was compared to changes in algal and bacterial cell numbers. Bacterial abundance was higher and dissolved organic carbon concentrations were lower in the diatom cultures infected with parasitic fungi. Among the DOM compounds, the concentrations of tryptophan-like material derived from algae were significantly lower and the concentrations of humic substance-like material were higher in the infected treatment. The parasitic fungi may have consumed tryptophan-like material and stimulated the release of humic substances. These results provide the first evidence that fungal infection may modulate algal-bacterial interactions, which are associated with changes in the nature of AOM.

Potential of Genomic Selection in Mass Selection Breeding of an Allogamous Crop: An Empirical Study to Increase Yield of Common Buckwheat.

To evaluate the potential of genomic selection (GS), a selection experiment with GS and phenotypic selection (PS) was performed in an allogamous crop, common buckwheat (Fagopyrum esculentum Moench). To indirectly select for seed yield per unit area, which cannot be measured on a single-plant basis, a selection index was constructed from seven agro-morphological traits measurable on a single plant basis. Over 3 years, we performed two GS and one PS cycles per year for improvement in the selection index. In GS, a prediction model was updated every year on the basis of genotypes of 14,598-50,000 markers and phenotypes. Plants grown from seeds derived from a series of generations of GS and PS populations were evaluated for the traits in the selection index and other yield-related traits. GS resulted in a 20.9% increase and PS in a 15.0% increase in the selection index in comparison with the initial population. Although the level of linkage disequilibrium in the breeding population was low, the target trait was improved with GS. Traits with higher weights in the selection index were improved more than those with lower weights, especially when prediction accuracy was high. No trait changed in an unintended direction in either GS or PS. The accuracy of genomic prediction models built in the first cycle decreased in the later cycles because the genetic bottleneck through the selection cycles changed linkage disequilibrium patterns in the breeding population. The present study emphasizes the importance of updating models in GS and demonstrates the potential of GS in mass selection of allogamous crop species, and provided a pilot example of successful application of GS to plant breeding.

Impact of Mislabeling on Genomic Selection in Cassava Breeding.

In plant breeding, humans occasionally make mistakes. Genomic selection is particularly prone to human error because it involves more steps than conventional phenotypic selection. The impact of human mistakes should be determined to evaluate the cost effectiveness of controlling human error in plant breeding. We used simulation to evaluate the impact of mislabeling, where marker scores from one plant are associated with the performance records of another plant in cassava (Manihot esculenta Crantz) breeding. Results showed that, although selection with mislabeling reduced genetic gains, scenarios including six levels of mislabeling (from 5 to 50%) persisted in achieving gain because mislabeling decreased the genetic variance lost from the population. Breeding populations with higher rates of mislabeling experienced lower selection intensity, resulting in higher genetic variance, which partially compensated for the mislabeling. For low mislabeling rates (10% or less), the increased genetic variance observed under mislabeling led to improved accuracy of the prediction model in later selection cycles. Large-scale mislabeling should therefore be prevented, but the value of preventing small-scale mislabeling depends on the effort already being invested in preventing the loss of genetic variance during the course of selection. In a program, such as the one we simulated, that makes no effort to avoid loss of genetic variance, small-scale mislabeling has a less negative effect than expected. We assume that negative effects would be greater if best practices to avoid genetic variance loss were already implemented.

Island-Model Genomic Selection for Long-Term Genetic Improvement of Autogamous Crops.

Acceleration of genetic improvement of autogamous crops such as wheat and rice is necessary to increase cereal production in response to the global food crisis. Population and pedigree methods of breeding, which are based on inbred line selection, are used commonly in the genetic improvement of autogamous crops. These methods, however, produce a few novel combinations of genes in a breeding population. Recurrent selection promotes recombination among genes and produces novel combinations of genes in a breeding population, but it requires inaccurate single-plant evaluation for selection. Genomic selection (GS), which can predict genetic potential of individuals based on their marker genotype, might have high reliability of single-plant evaluation and might be effective in recurrent selection. To evaluate the efficiency of recurrent selection with GS, we conducted simulations using real marker genotype data of rice cultivars. Additionally, we introduced the concept of an "island model" inspired by evolutionary algorithms that might be useful to maintain genetic variation through the breeding process. We conducted GS simulations using real marker genotype data of rice cultivars to evaluate the efficiency of recurrent selection and the island model in an autogamous species. Results demonstrated the importance of producing novel combinations of genes through recurrent selection. An initial population derived from admixture of multiple bi-parental crosses showed larger genetic gains than a population derived from a single bi-parental cross in whole cycles, suggesting the importance of genetic variation in an initial population. The island-model GS better maintained genetic improvement in later generations than the other GS methods, suggesting that the island-model GS can utilize genetic variation in breeding and can retain alleles with small effects in the breeding population. The island-model GS will become a new breeding method that enhances the potential of genomic selection in autogamous crops, especially bringing long-term improvement.

Rapid genotyping with DNA micro-arrays for high-density linkage mapping and QTL mapping in common buckwheat (Fagopyrum esculentum Moench).

For genetic studies and genomics-assisted breeding, particularly of minor crops, a genotyping system that does not require a priori genomic information is preferable. Here, we demonstrated the potential of a novel array-based genotyping system for the rapid construction of high-density linkage map and quantitative trait loci (QTL) mapping. By using the system, we successfully constructed an accurate, high-density linkage map for common buckwheat (Fagopyrum esculentum Moench); the map was composed of 756 loci and included 8,884 markers. The number of linkage groups converged to eight, which is the basic number of chromosomes in common buckwheat. The sizes of the linkage groups of the P1 and P2 maps were 773.8 and 800.4 cM, respectively. The average interval between adjacent loci was 2.13 cM. The linkage map constructed here will be useful for the analysis of other common buckwheat populations. We also performed QTL mapping for main stem length and detected four QTL. It took 37 days to process 178 samples from DNA extraction to genotyping, indicating the system enables genotyping of genome-wide markers for a few hundred buckwheat plants before the plants mature. The novel system will be useful for genomics-assisted breeding in minor crops without a priori genomic information.

Codon pairs of the HIV-1 vif gene correlate with CD4+ T cell count.

BACKGROUND: The human APOBEC3G (A3G) protein activity is associated with innate immunity against HIV-1 by inducing high rates of guanosines to adenosines (G-to-A) mutations (viz., hypermutation) in the viral DNA. If hypermutation is not enough to disrupt the reading frames of viral genes, it may likely increase the HIV-1 diversity. To counteract host innate immunity HIV-1 encodes the Vif protein that binds A3G protein and form complexes to be degraded by cellular proteolysis. METHODS: Here we studied the pattern of substitutions in the vif gene and its association with clinical status of HIV-1 infected individuals. To perform the study, unique vif gene sequences were generated from 400 antiretroviral-naïve individuals. RESULTS: The codon pairs: 78-154, 85-154, 101-157, 105-157, and 105-176 of vif gene were associated with CD4+ T cell count lower than 500 cells per mm(3). Some of these codons were located in the (81)LGQGVSIEW(89) region and within the BC-Box. We also identified codons under positive selection clustered in the N-terminal region of Vif protein, between (21)WKSLVK(26) and (40)YRHHY(44) regions (i.e., 31, 33, 37, 39), within the BC-Box (i.e., 155, 159) and the Cullin5-Box (i.e., 168) of vif gene. All these regions are involved in the Vif-induced degradation of A3G/F complexes and the N-terminal of Vif protein binds to viral and cellular RNA. CONCLUSIONS: Adaptive evolution of vif gene was mostly to optimize viral RNA binding and A3G/F recognition. Additionally, since there is not a fully resolved structure of the Vif protein, codon pairs associated with CD4+ T cell count may elucidate key regions that interact with host cell factors. Here we identified and discriminated codons under positive selection and codons under functional constraint in the vif gene of HIV-1.