Genomic selection is a promising molecular breeding strategy enhancing genetic gain per unit time. The objectives of our study were to (1) explore the prediction accuracy of genomic selection for plant height and yield per plant in soybean [Glycine max (L.) Merr.], (2) discuss the relationship between prediction accuracy and numbers of markers, and (3) evaluate the effect of marker preselection based on different methods on the prediction accuracy. Our study is based on a population of 235 soybean varieties which were evaluated for plant height and yield per plant at multiple locations and genotyped by 5361 single nucleotide polymorphism markers. We applied ridge regression best linear unbiased prediction coupled with fivefold cross-validations and evaluated three strategies of marker preselection. For plant height, marker density and marker preselection procedure impacted prediction accuracy only marginally. In contrast, for grain yield, prediction accuracy based on markers selected with a haplotype block analyses-based approach increased by approximately 4 % compared with random or equidistant marker sampling. Thus, applying marker preselection based on haplotype blocks is an interesting option for a cost-efficient implementation of genomic selection for grain yield in soybean breeding.

Anthocyanins protect plants from biotic and abiotic stressors and provide great health benefits to consumers. In this study, we cloned four genes (Red Lettuce Leaves 1 to 4: RLL1 to RLL4) that contribute to colour variations in lettuce. The RLL1 gene encodes a bHLH transcription factor, and a 5-bp deletion in some cultivars abolishes its function to activate the anthocyanin biosynthesis pathway. The RLL2 gene encodes an R2R3-MYB transcription factor, which was derived from a duplication followed by mutations in its promoter region. The RLL3 gene encodes an R2-MYB transcription factor, which down-regulates anthocyanin biosynthesis through competing with RLL2 for interaction with RLL1; a mis-sense mutation compromises the capacity of RLL3 to bind RLL1. The RLL4 gene encodes a WD-40 transcription factor, homologous to the RUP genes suppressing the UV-B signal transduction pathway in Arabidopsis; a mis-sense mutation in rll4 attenuates its suppressing function, leading to a high concentration of anthocyanins. Sequence analysis of the RLL1-RLL4 genes from wild and cultivated lettuce showed that their function-changing mutations occurred after domestication. The mutations in rll1 disrupt anthocyanin biosynthesis, while the mutations in RLL2, rll3 and rll4 activate anthocyanin biosynthesis, showing disruptive selection for leaf colour during domestication of lettuce. The characterization of multiple polymorphic genes in this study provides the necessary molecular resources for the rational breeding of lettuce cultivars with distinct levels of red pigments and green cultivars with high levels of health-promoting flavonoids.

Hybridization and resulting introgression are important processes shaping the tree of life and appear to be far more common than previously thought. However, how the genome evolution was shaped by various genetic and evolutionary forces after hybridization remains unresolved. Here we used whole-genome resequencing data of 227 individuals from multiple widespread Populus species to characterize their contemporary patterns of hybridization and to quantify genomic signatures of past introgression. We observe a high frequency of contemporary hybridization and confirm that multiple previously ambiguous species are in fact F1 hybrids. Seven species were identified, which experienced different demographic histories that resulted in strikingly varied efficacy of selection and burdens of deleterious mutations. Frequent past introgression has been found to be a pervasive feature throughout the speciation of these Populus species. The retained introgressed regions, more generally, tend to contain reduced genetic load and to be located in regions of high recombination. We also find that in pairs of species with substantial differences in effective population size, introgressed regions are inferred to have undergone selective sweeps at greater than expected frequencies in the species with lower effective population size, suggesting that introgression likely have higher potential to provide beneficial variation for species with small populations. Our results, therefore, illustrate that demography and recombination have interplayed with both positive and negative selection in determining the genomic evolution after hybridization.

Worldwide distributed coccidiosis is caused by infection of both Eimeria species and Cystoisospora in the host intestine and causes huge economic losses to the livestock industry, especially the poultry industry. The control of such diseases relies mainly on chemoprophylaxis with anticoccidials, which has led to a very common drug resistance in this field. However, the genetic mechanisms underlying resistance to many anticoccidial drugs remain unknown. In this study, strains of E. tenella resistant to 250 mg/kg monensin were generated and characterized. Forward genetic approaches based on pooled genome sequencing, including experimental evolution and linkage group selection, were used to locate candidate targets responsible for resistance to monensin and diclazuril in E. tenella. A total of 16 nonsynonymous mutants in protein-coding genes were identified in monensin-resistant strains, and two genomic regions with strong selection signals were also detected in diclazuril-resistant strains. Our study reveals the genetic characterization of the experimental evolution and linkage group selection in Eimeria species, and also provides important information that contributes to the understanding of the molecular mechanism of drug resistance in coccidia.

As a diverse species, watermelon [Citrullus lanatus (Thunb.) Matsum. &Nakai var. lanatus] has different kinds of fruit sizes, shapes, flesh colors and skin colors. Skin color is among the major objectives for breeding. Yellow skin is an important trait in watermelon, but the underlying genetic mechanism is unknown. In this study, we identified a locus for yellow skin through BSA-seq and GWAS. A segregation analysis in F2 and BC1 populations derived from a cross of two inbred lines '94E1'(yellow skin) and 'Qingfeng'(green skin) suggested that skin color is a qualitative trait. BSA-seq mapping confirmed the locus in the F2 population, which was detected on chromosome 4 by GWAS among 330 varieties. Several major markers, namely, 15 CAPS markers, 6 SSR markers and 2 SNP markers, were designed to delimit the region to 59.8 kb region on chromosome 4. Utilizing the two populations consisting of 10 yellow and 10 green skin watermelons, we found a tightly linked functional SNP marker for the yellow skin phenotype. The application of this marker as a selection tool in breeding programs will help to improve the breeder's ability to make selections at early stages of growth, thus accelerating the breeding program.

Exploitation of crop heterosis is crucial for increasing global agriculture production. However, the quantitative genomic analysis of heterosis was lacking, and there is currently no effective prediction tool to optimize cross-combinations. Here 2,839 rice hybrid cultivars and 9,839 segregation individuals were resequenced and phenotyped. Our findings demonstrated that indica-indica hybrid-improving breeding was a process that broadened genetic resources, pyramided breeding-favorable alleles through combinatorial selection and collaboratively improved both parents by eliminating the inferior alleles at negative dominant loci. Furthermore, we revealed that widespread genetic complementarity contributed to indica-japonica intersubspecific heterosis in yield traits, with dominance effect loci making a greater contribution to phenotypic variance than overdominance effect loci. On the basis of the comprehensive dataset, a genomic model applicable to diverse rice varieties was developed and optimized to predict the performance of hybrid combinations. Our data offer a valuable resource for advancing the understanding and facilitating the utilization of heterosis in rice.

Raf-1 kinase inhibitor protein (RKIP) plays a critical role in tumor development by regulating cell functions such as invasion, apoptosis and differentiation. Down-regulation of RKIP expression has been implicated in the development and progression of renal cell carcinoma (RCC). Herein, we hypothesized that genetic polymorphisms in RKIP might be associated with susceptibility and progression of RCC.

Foxtail millet [Setaria italica (L.) P. Beauv.], a crop of historical importance in China, has been adopted as a model crop for studying C-4 photosynthesis, stress biology and biofuel traits. Construction of a high density genetic map and identification of stable quantitative trait loci (QTL) lay the foundation for marker-assisted selection for agronomic traits and yield improvement.

No abstract available

LncARSR (lncRNA Activated in RCC with Sunitinib Resistance, ENST00000424980) is a newly identified lncRNA to promote the sunitinib resistance of renal cell carcinoma (RCC), which may contribute to tumorigenesis and progression. This study aimed to explore the association of lncARSR tagSNPs with the risk and prognosis of RCC.

Aiming to investigate whether genetic risk factors (GRFs) for fracture and bone mineral density (BMD) identified from people of European descent can help improve the prediction of osteoporotic fracture (OF) risk and BMD in Chinese populations, we built assessment models for femoral neck (FN)-fracture prediction and BMD value prediction using 700 elderly Chinese Han subjects and 1,620 unrelated Chinese Han subjects, respectively. 17 fracture-associated genes and 82 FN-BMD associated genes identified in people of European descent were used to build a logistic regression model with clinical risk factors (CRFs) for FN-fracture prediction in Chinese. Meanwhile 107 BMD-associated genes from people of European descent were used to build a multiple linear regression model with CRFs for BMD prediction in Chinese. A Lasso algorithm was employed for informative SNP selection to construct the genetic risk score (GRS) with ten-fold cross-validation. The results showed that, adding fracture GRF and FN-BMD GRF to the model with CRFs, the area under the receiver operating characteristic curve (AUC) decrease from 0.653 to 0.587 and 0.588, respectively, for FN fracture prediction. 62.3% and 61.8% of the risk variation were explained by the Model with CRFs and fracture GRF and by the Model with CRFs and FN-BMD GRF, respectively, as compared to 65.5% in the Model with CRFs only. The net reclassification improvement (NRI) index in the reclassification analysis is 0.56% (P = 0.57) and 1.13% (P = 0.29), respectively. There is no significant difference either between the performance of the model with CRFs and that of the model with both CRFs and GRF for BMD prediction. We concluded that, in the current study, GRF of fracture identified in people of European descent does not contributes to improve the fracture prediction in Chinese; and GRF of BMD from people of European descent cannot help improve the accuracy of the fracture prediction in Chinese perhaps partially because GRF of BMD from people of European descent may not contribute to BMD prediction in Chinese. This study highlights the limited utility of the current genetics studies largely focused on people of European descent for disease or risk factor prediction in other ethnic groups, and calls for more and larger scale studies focused on other ethnic groups.

Six novel fruit weight QTLs were identified in tomato using multiple bi-parental populations developed from ancestral accessions. Beneficial alleles at these loci arose in semi-domesticated subpopulations and were likely left behind. This study paves the way to introgress these alleles into breeding programs. The size and weight of edible organs have been strongly selected during crop domestication. Concurrently, human have also focused on nutritional and cultural characteristics of fruits and vegetables, at times countering selective pressures on beneficial size and weight alleles. Therefore, it is likely that novel improvement alleles for organ weight still segregate in ancestral germplasm. To date, five domestication and diversification genes affecting tomato fruit weight have been identified, yet the genetic basis for increases in weight has not been fully accounted for. We found that fruit weight increased gradually during domestication and diversification, and semi-domesticated subpopulations featured high phenotypic and nucleotide diversity. Columella and septum fruit tissues were proportionally increased, suggesting targeted selection. We developed twenty-one F2 populations with parents fixed for the known fruit weight genes, corresponding to putative key transitions from wild to fully domesticated tomatoes. These parents also showed differences in fruit weight attributes as well as the developmental timing of size increase. A subset of populations was targeted for QTL-seq, leading to the identification of six uncloned fruit weight QTLs. Three QTLs, located on chromosomes 1, 2 and 3, were subsequently validated by progeny testing. By exploring the segregation of the known fruit weight genes and the identified QTLs, we estimated that most beneficial alleles in the newly identified loci arose in semi-domesticated subpopulations from South America and were not likely transmitted to fully domesticated landraces. Therefore, these alleles could be incorporated into breeding programs using the germplasm and genetic resources identified in this study.

Ginseng, which contains ginsenosides as bioactive compounds, has been regarded as an important traditional medicine for several millennia. However, the genetic background of ginseng remains poorly understood, partly because of the plant's large and complex genome composition. We report the entire genome sequence of Panax ginseng using next-generation sequencing. The 3.5-Gb nucleotide sequence contains more than 60% repeats and encodes 42 006 predicted genes. Twenty-two transcriptome datasets and mass spectrometry images of ginseng roots were adopted to precisely quantify the functional genes. Thirty-one genes were identified to be involved in the mevalonic acid pathway. Eight of these genes were annotated as 3-hydroxy-3-methylglutaryl-CoA reductases, which displayed diverse structures and expression characteristics. A total of 225 UDP-glycosyltransferases (UGTs) were identified, and these UGTs accounted for one of the largest gene families of ginseng. Tandem repeats contributed to the duplication and divergence of UGTs. Molecular modeling of UGTs in the 71st, 74th, and 94th families revealed a regiospecific conserved motif located at the N-terminus. Molecular docking predicted that this motif captures ginsenoside precursors. The ginseng genome represents a valuable resource for understanding and improving the breeding, cultivation, and synthesis biology of this key herb.

Leaf heading is an important and economically valuable horticultural trait in many vegetables. The formation of a leafy head is a specialized leaf morphogenesis characterized by the emergence of the enlarged incurving leaves. However, the transcriptional regulation mechanisms underlying the transition to leaf heading remain unclear. We carried out large-scale time-series transcriptome assays covering the major vegetative growth phases of two headingBrassica crops, Chinese cabbage and cabbage, with the non-heading morphotype Taicai as the control. A regulatory transition stage that initiated the heading process is identified, accompanied by a developmental switch from rosette leaf to heading leaf in Chinese cabbages. This transition did not exist in the non-heading control. Moreover, we reveal that the heading transition stage is also conserved in the cabbage clade. Chinese cabbage acquired through domestication a leafy head independently from the origins of heading in other cabbages; phylogenetics supports that the ancestor of all cabbages is non-heading. The launch of the transition stage is closely associated with the ambient temperature. In addition, examination of the biological activities in the transition stage identified the ethylene pathway as particularly active, and we hypothesize that this pathway was targeted for selection for domestication to form the heading trait specifically in Chinese cabbage. In conclusion, our findings on the transcriptome transition that initiated the leaf heading in Chinese cabbage and cabbage provide a new perspective for future studies of leafy head crops.

Lean body mass (LBM) is a complex trait for human health. To identify genomic loci underlying LBM, we performed a gene-based genome-wide association study of lean mass index (LMI) in 1000 unrelated Caucasian subjects, and replicated in 2283 unrelated Caucasians subjects. Gene-based association analyses highlighted the significant associations of three genes UQCR, TCF3 and MBD3 in one single locus 19p13.3 (discovery p = 6.10 × 10-5, 1.65 × 10-4 and 1.10 × 10-4; replication p = 2.21 × 10-3, 1.84 × 10-3 and 6.95 × 10-3; combined p = 2.26 × 10-6, 4.86 × 10-6 and 1.15 × 10-5, respectively). These results, together with the known functional relevance of the three genes to LMI, suggested that the 19p13.3 region containing UQCR, TCF3 and MBD3 genes was a novel locus underlying lean mass variation.

The genetic improvement of disease resistance in poultry continues to be a challenge. To identify candidate genes and loci responsible for these traits, genome-wide association studies using the chicken 60k high density single nucleotide polymorphism (SNP) array for six immune traits, total serum immunoglobulin Y (IgY) level, numbers of, and the ratio of heterophils and lymphocytes, and antibody responses against Avian Influenza Virus (AIV) and Sheep Red Blood Cell (SRBC), were performed. RT-qPCR was used to quantify the relative expression of the identified candidate genes. Nine significantly associated SNPs (P < 2.81E-06) and 30 SNPs reaching the suggestively significant level (P < 5.62E-05) were identified. Five of the 10 SNPs that were suggestively associated with the antibody response to SRBC were located within or close to previously reported QTL regions. Fifteen SNPs reached a suggestive significance level for AIV antibody titer and seven were found on the sex chromosome Z. Seven suggestive markers involving five different SNPs were identified for the numbers of heterophils and lymphocytes, and the heterophil/lymphocyte ratio. Nine significant SNPs, all on chromosome 16, were significantly associated with serum total IgY concentration, and the five most significant were located within a narrow region spanning 6.4kb to 253.4kb (P = 1.20E-14 to 5.33E-08). After testing expression of five candidate genes (IL4I1, CD1b, GNB2L1, TRIM27 and ZNF692) located in this region, changes in IL4I1, CD1b transcripts were consistent with the concentrations of IgY, while abundances of TRIM27 and ZNF692 showed reciprocal changes to those of IgY concentrations. This study has revealed 39 SNPs associated with six immune traits (total serum IgY level, numbers of, and the ratio of heterophils and lymphocytes, and antibody responses against AIV and SRBC) in Beijing-You chickens. The narrow region spanning 247kb on chromosome 16 is an important QTL for serum total IgY concentration. Five candidate genes related to IgY level validated here are novel and may play critical roles in the modulation of immune responses. Potentially useful candidate SNPs for marker-assisted selection for disease resistance are identified. It is highly likely that these candidate genes play roles in various aspects of the immune response in chickens.

Plasminogen activator inhibitor -2 (PAI-2) is an important molecular that plays a crucial role in vascular homeostasis and constitutes a critical response mechanism to cardiovascular injury, such as atherosclerosis, coronary artery disease (CAD).

The suppressor of the cytokine signaling (SOCS) family of proteins play an essential role in inhibiting cytokine receptor signaling by regulating immune signal pathways. Although SOCS gene functions have been examined extensively, no comprehensive study has been performed on this gene family's molecular evolution in reptiles. In this study, we identified eight canonical SOCS genes using recently-published reptilian genomes. We used phylogenetic analysis to determine that the SOCS genes had highly conserved evolutionary dynamics that we classified into two types. We identified positive SOCS4 selection signals in whole reptile lineages and SOCS2 selection signals in the crocodilian lineage. Selective pressure analyses using the branch model and Z-test revealed that these genes were under different negative selection pressures compared to reptile lineages. We also concluded that the nature of selection pressure varies across different reptile lineages on SOCS3, and the crocodilian lineage has experienced rapid evolution. Our results may provide a theoretical foundation for further analyses of reptilian SOCS genes' functional and molecular mechanisms, as well as their roles in reptile growth and development.

The hedgehog signaling pathway plays a vital role in human and animal patterning and cell proliferation during the developmental process. The hedgehog gene family of vertebrate species includes three genes, Shh, Dhh, and Ihh, which possess different functions and expression patterns. Despite the importance of hedgehog genes, genomic evidence of this gene family in reptiles is lacking. In this study, the available genomes of a number of representative reptile species were explored by utilizing adaptive evolutionary analysis methods to characterize the evolutionary patterns of the hedgehog gene family. Altogether, 33 sonic hedgehog (Shh), 25 desert hedgehog (Dhh), and 20 Indian hedgehog (Ihh) genes were obtained from reptiles, and six avian and five mammalian sequences were added to the analysis. The phylogenetic maximum likelihood (ML) tree of the Shh, Dhh, and Ihh genes revealed a similar topology, which is approximately consistent with the traditional taxonomic group. No shared positive selection site was identified by the PAML site model or the three methods in the Data Monkey Server. Branch model and Clade model C analyses revealed that the Dhh and Ihh genes experienced different evolutionary forces in reptiles and other vertebrates, while the Shh gene was not significantly different in terms of selection pressure. The different evolutionary rates of the Dhh and Ihh genes suggest that these genes may be potential contributors to the discrepant sperm and body development of different clades. The different adaptive evolutionary history of the Shh, Dhh, and Ihh genes among reptiles may be due to their different functions in regulating cellular events of development from the embryonic stages to adulthood. Overall, this study has provided meaningful information regarding the evolution of the hedgehog gene family in reptiles and a theoretical foundation for further analyses on the functional and molecular mechanisms that have shaped the reptilian hedgehog genes.

The inhibitor of differentiation/DNA-binding (ID) is a member of the helix-loop-helix (HLH) transcription factor family, and plays a role in tumorigenesis, invasiveness and angiogenesis. The aims were to investigate the expression patterns and prognostic values of individual ID family members in lung cancer, and the potential functional roles. The expression levels of ID family were assessed using the Oncomine online database and GEPIA database. Furthermore, the prognostic value of ID family members was evaluated using the Kaplan-Meier plotter database. The genetic mutations of ID family members were investigated using the cBioPortal database. Moreover, enrichment analysis was performed using STRING database and Funrich software. It was found that all the ID family members were significantly down-regulated in lung cancer. Prognostic results indicated that low mRNA expression levels of ID1 or increased mRNA expression levels of ID2/3/4 were associated with improved overall survival, first progression and post progression survival. Additionally, genetic mutations of ID family members were identified in lung cancer, and it was suggested that amplification and deep deletion were the main mutation types. Furthermore, functional enrichment analysis results suggested that ID1/2/4 were significantly enriched in 'regulation of nucleobase, nucleoside, nucleotide and nucleic acid metabolism' for biological process, 'transcription factor activity' for molecular function and 'HLH domain' for protein domain. However, it was found that ID3 was not enriched in the above functions. The aberrant expression of ID family members may affect the occurrence and prognosis of lung cancer, and may be related to cell metabolism and transcriptional regulation.