Citrus is a large genus that includes several major cultivated species, including C. sinensis (sweet orange), Citrus reticulata (tangerine and mandarin), Citrus limon (lemon), Citrus grandis (pummelo) and Citrus paradisi (grapefruit). In 2009, the global citrus acreage was 9 million hectares and citrus production was 122.3 million tons (FAO statistics, see URLs), which is the top ranked among all the fruit crops. Among the 10.9 million tons (valued at $9.3 billion) of citrus products traded in 2009, sweet orange accounted for approximately 60% of citrus production for both fresh fruit and processed juice consumption (FAO statistics, see URLs). Moreover, citrus fruits and juice are the prime human source of vitamin C, an important component of human nutrition.Citrus fruits also have some unique botanical features, such as nucellar embryony (nucellus cells can develop into apomictic embryos that are genetically identical to mother plant). Consequently, somatic embryos grow much more vigorously than the zygotic embryos in seeds such that seedlings are essentially clones of the maternal parent. Such citrus-unique characteristics have hindered the study of citrus genetics and breeding improvement 1,2 . Complete genome sequences would provide valuable genetic resources for improving citrus crops.Citrus is believed to be native to southeast Asia 3-5 , and cultivation of fruit crops occurred at least 4,000 years ago 3,6 . The genetic origin of the sweet orange is not clear, although there are some speculations that sweet orange might be derived from interspecific hybridization of some primitive citrus species 7,8 . Citrus is also in the order Sapindales, a sister order to the Brassicales in the Malvidae, making it valuable for comparative genomics studies with the model plant Arabidopsis.We aimed to sequence the genome of Valencia sweet orange (C. sinensis cv. Valencia), one of the most important sweet orange varieties cultivated worldwide and grown primarily for orange juice production. Normal sweet oranges are diploids, with nine pairs of chromosomes and an estimated genome size of ~367 Mb 9 . To reduce the complexity of the sequenced genome, we obtained a doublehaploid (dihaploid) line derived from the anther culture of Valencia sweet orange 10 . We first generated whole-genome shotgun pairedend-tag sequence reads from the dihaploid genomic DNA and built a de novo assembly as the citrus reference genome; we then produced shotgun sequencing reads from the parental diploid DNA and mapped the sequences to the haploid reference genome to obtain the complete genome information for Valencia sweet orange. In addition, we conducted comprehensive transcriptome sequencing analyses for four representative tissues using shotgun RNA sequencing (RNA-Seq) to capture all transcribed sequences and paired-end-tag RNA sequencing (RNA-PET) to demarcate the 5′ and 3′ ends of all transcripts. On the basis of the DNA and RNA sequencing data, we characterized the orange genome for its gene content, heterozygosity and evolutionary features. ...

The emergence of apomixis-the transition from sexual to asexual reproduction-is a prominent feature of modern citrus. Here we de novo sequenced and comprehensively studied the genomes of four representative citrus species. Additionally, we sequenced 100 accessions of primitive, wild and cultivated citrus. Comparative population analysis suggested that genomic regions harboring energy- and reproduction-associated genes are probably under selection in cultivated citrus. We also narrowed the genetic locus responsible for citrus polyembryony, a form of apomixis, to an 80-kb region containing 11 candidate genes. One of these, CitRWP, is expressed at higher levels in ovules of polyembryonic cultivars. We found a miniature inverted-repeat transposable element insertion in the promoter region of CitRWP that cosegregated with polyembryony. This study provides new insights into citrus apomixis and constitutes a promising resource for the mining of agriculturally important genes.

Mandarin (Citrus reticulata) is one of the most important citrus crops worldwide. Its domestication is believed to have occurred in South China, which has been one of the centers of mandarin cultivation for four millennia. We collected natural wild populations of mandarin around the Nanling region and cultivated landraces in the vicinity. We found that the citric acid level was dramatically reduced in cultivated mandarins. To understand genetic basis of mandarin domestication, we de novo assembled a draft genome of wild mandarin and analyzed a set of 104 citrus genomes. We found that the Mangshan mandarin is a primitive type and that two independent domestication events have occurred, resulting in two groups of cultivated mandarins (MD1 and MD2) in the North and South Nanling Mountains, respectively. Two bottlenecks and two expansions of effective population size were identified for the MD1 group of cultivated mandarins. However, in the MD2 group there was a long and continuous decrease in the population size. MD1 and MD2 mandarins showed different patterns of interspecific introgression from cultivated pummelo species. We identified a region of high divergence in an aconitate hydratase (ACO) gene involved in the regulation of citrate content, which was possibly under selection during the domestication of mandarin. This study provides concrete genetic evidence for the geographical origin of extant wild mandarin populations and sheds light on the domestication and evolutionary history of mandarin.

HighlightThe rice SLG gene, functioning as homomers, plays essential roles in regulating grain size and leaf angle via modulation of brassinosteroid homeostasis.

Current evidences show that copy number variations (CNVs) are linked to complex phenotypic traits. Leptin receptor (LEPR) gene plays a critical role in energy homeostasis and fat development and re-sequencing of the cattle genome revealed the CNV region (herein referred to as "I3 DNA") within the LEPR intron 3. In the present study, we qualified copy numbers of I3 DNA within LEPR gene in four cattle breeds (Qinchuan, Nanyang, Jinnan and Xianan) by quantitative PCR, and explored their impacts on LEPR gene expression and phenotypic traits in Qinchuan and Nanyang cattle. The results showed that more individuals in Nanyang are with loss of the I3 DNA copy number than that in the others. Additionally, I3 DNA CNVs exhibited a significant negative correlation with LEPR gene expression (P < 0.05). Association analysis showed that gain/normal copy number types performed better traits of body weight, body height and body length than the loss type in Nanyang. To the best of our knowledge, this is the first evidence of the association between LEPR CNVs and cattle traits, and this may help deep understanding of the function of CNVs which may be promising markers for beef cattle breeding and genetics.

Tung tree (Vernicia fordii) is an economically important woody oil plant that produces tung oil rich in eleostearic acid. Here, we report a high-quality chromosome-scale genome sequence of tung tree. The genome sequence was assembled by combining Illumina short reads, Pacific Biosciences single-molecule real-time long reads, and Hi-C sequencing data. The size of tung tree genome is 1.12 Gb, with 28,422 predicted genes and over 73% repeat sequences. The V. fordii underwent an ancient genome triplication event shared by core eudicots but no further whole-genome duplication in the subsequent ca. 34.55 million years of evolutionary history of the tung tree lineage. Insertion time analysis revealed that repeat-driven genome expansion might have arisen as a result of long-standing long terminal repeat retrotransposon bursts and lack of efficient DNA deletion mechanisms. The genome harbors 88 resistance genes encoding nucleotide-binding sites; 17 of these genes may be involved in early-infection stage of Fusarium wilt resistance. Further, 651 oil-related genes were identified, 88 of which are predicted to be directly involved in tung oil biosynthesis. Relatively few phosphoenolpyruvate carboxykinase genes, and synergistic effects between transcription factors and oil biosynthesis-related genes might contribute to the high oil content of tung seed. The tung tree genome constitutes a valuable resource for understanding genome evolution, as well as for molecular breeding and genetic improvements for oil production.

Naturally, resistant crop germplasms are important resources for managing the issues of agricultural product safety and environment deterioration. We found a spontaneous mutant of ‘Newhall’ navel orange (Citrus sinensis Osbeck) (MT) with broad-spectrum protections against fungal pathogens in the orchard, postharvest-storage, and artificial inoculation conditions. To understand the defense mechanism of MT fruit, we constructed a genome-scale metabolic network that integrated metabolome and transcriptome datasets. The coordinated transcriptomic and metabolic data were enriched in two sub-networks, showing the decrease in very long chain fatty acid (by 41.53%) and cuticular wax synthesis (by 81.34%), and increase in the synthesis of jasmonic acid (JA) (by 95.23%) and JA-induced metabolites such as 5-dimethylnobietin (by 28.37%) in MT. Furthermore, cytological and biochemical analyses confirmed that the response to fungal infection in MT was independent of wax deficiency and was correlated with the levels of jasmonates, and the expression of plant defensin gene PDF1.2. Results of exogenous application of MeJA and JA inhibitors such as propyl gallate proved that JA-mediated defense contributes to the strong tolerance against pathogens in MT. Our results indicated that jasmonate biosynthesis and signaling are stimulated by the fatty acid redirection of MT, and participate in the tolerance of pathogenic fungi.

BackgroundGrain size, which is determined by grain length, grain width, and grain thickness, is an important determinant for grain yield in rice. Identification and characterization of new genes that are associated with grain size will be helpful for the improvement of grain yield in rice.ResultsWe characterized the grain size mutant, larger grain size 1 (lgs1), derived from rice activation-tagged T-DNA insertion lines. Histological analysis showed that increased cell numbers in the longitudinal direction of spikelet hulls was responsible for the grain mutant phenotype in lgs1. Quantitative real-time PCR (qRT-PCR) analysis further showed that the expression levels of genes associated with the cell cycle in the young panicles of the lgs1 were higher than those in the wild type (WT), which might result in the increased cell numbers in lgs1 spikelet hulls. Insertion site analysis together with transgenic experiments confirmed that the lgs1 phenotype was caused by enhanced expression of truncated OsbHLH107, corresponding to the nucleotide (nt) 331–846 region (i.e., the transcriptional activation region of OsbHLH107) of the OsbHLH107 coding sequence (CDS). OsbHLH107 is a nucleus-localized bHLH transcription factor, which can form a homodimer with itself. Phylogenetic analysis showed that OsbHLH107 belonged to the same subfamily as OsPILs. OsPIL13 (OsPIL1) and OsPIL16 (APG) were reported to regulate grain size in rice. By transgenic experiments, we found that OsPIL11 could also regulate grain size.ConclusionWe concluded that OsbHLH107 and its homologs are important regulators of grain size development and might be useful for grain yield improvement in rice.Electronic supplementary materialThe online version of this article (10.1186/s12284-018-0237-y) contains supplementary material, which is available to authorized users.