Chengzhang Liu

Water-yield relations and water-use efficiency of winter wheat in the North China Plain

Abstract. Limited precipitation restricts yield of winter wheat (Triticum aestivum L.) grown in the North China Plain. Water stress effects on yield can be avoided or minimized by application of irrigation. We examined the multiseasonal irrigation experiments in four locations of the piedmont and lowland in the region, and developed crop water-stress sensitivity index, relationship between seasonal evapotranspiration (ET) and yield, and crop water production functions. By relating relative yield to relative ET deficit, we found that the crop was more sensitive to water stress from stem elongation to heading and from heading to milking. For limited irrigation, irrigation is recommended during the stages sensitive to water stress. Grain yield was 258–322 g m−2 in the piedmont and 260–280 g m−2 in the lowland under rainfed conditions. The corresponding seasonal ET was 242–264 mm in the piedmont and 247–281 mm in the lowland. Irrigation significantly increased seasonal ET and therefore grain yield as a result of increased kernel numbers per m−2 and kernels per ear. On average, one irrigation increased grain yield by 21–43% and two to four irrigations by 60–100%. Grain yield was linearly related to seasonal ET with a slope of 1.15 kg m−3 in the lowland and 1.73 kg m−3 in the piedmont. Water-use efficiency was 0.98–1.22 kg m−3 for rainfed wheat and 1.20–1.40 kg m−3 for the wheat irrigated 2–4 times. Grain yield response to the amount of irrigation (IRR) was developed using a quadratic function and used to analyze different irrigation scenarios. To achieve the maximum grain yield, IRR was 240 mm in the piedmont and 290 mm in the lowland. When the maximum net profit was achieved, IRR was 195 mm and 250 mm in the piedmont and lowland, respectively. The yield response curve to IRR showed a plateau over a large range of IRR, indicating a great potential in saving IRR while maintaining reasonable high levels of grain yield.

SNP discovery in the transcriptome of white Pacific shrimp Litopenaeus vannamei by next generation sequencing.

The application of next generation sequencing technology has greatly facilitated high throughput single nucleotide polymorphism (SNP) discovery and genotyping in genetic research. In the present study, SNPs were discovered based on two transcriptomes of Litopenaeus vannamei (L. vannamei) generated from Illumina sequencing platform HiSeq 2000. One transcriptome of L. vannamei was obtained through sequencing on the RNA from larvae at mysis stage and its reference sequence was de novo assembled. The data from another transcriptome were downloaded from NCBI and the reads of the two transcriptomes were mapped separately to the assembled reference by BWA. SNP calling was performed using SAMtools. A total of 58,717 and 36,277 SNPs with high quality were predicted from the two transcriptomes, respectively. SNP calling was also performed using the reads of two transcriptomes together, and a total of 96,040 SNPs with high quality were predicted. Among these 96,040 SNPs, 5,242 and 29,129 were predicted as non-synonymous and synonymous SNPs respectively. Characterization analysis of the predicted SNPs in L. vannamei showed that the estimated SNP frequency was 0.21% (one SNP per 476 bp) and the estimated ratio for transition to transversion was 2.0. Fifty SNPs were randomly selected for validation by Sanger sequencing after PCR amplification and 76% of SNPs were confirmed, which indicated that the SNPs predicted in this study were reliable. These SNPs will be very useful for genetic study in L. vannamei, especially for the high density linkage map construction and genome-wide association studies.

Proteomic analysis of differentially expressed proteins in lymphoid organ of Fenneropenaeus chinensis response to Vibrio anguillarum stimulation

To gain an insight into the function of shrimp lymphoid organ at protein level, we analyzed the proteome of lymphoid organ in healthy Chinese shrimp Fenneropenaeus chinensis (F. chinensis) through two-dimensional gel electrophoresis (2-DE) based proteomic approach. A total of 95 spots representing 75 protein entries were identified by liquid chromatography tandem mass spectrometry (LC-MS/MS) with both online and in-house database. According to Gene Ontology (GO) annotation of biological process, the identified proteins were classified into 13 categories. Among them, approximately 36% of proteins related to cytoskeleton are noticeable. Then, a comparative proteomic approach was employed to investigate the differentially expressed proteins in lymphoid organ of Vibrio anguillarum-challenged F. chinensis. At 24 h post-injection (hpi), 17 differentially expressed protein spots were successfully identified, including 4 up-regulated protein spots (represent 4 proteins: cathepsin L, protein similar to squid CG16901-PC, protein kinase C and protein similar to T-complex Chaperonin 5 CG8439-PA), and 13 down-regulated protein spots (represent 9 proteins: actin, beta-actin, cytoplasmic actin CyII, alpha tubulin, beta tubulin, protein similar to proteasome delta, vacuolar ATP synthase subunit B, elongation factor 2, carboxypeptidase B). These data may help us to understand the function of lymphoid organ and the molecular immune mechanism of shrimp responsive to pathogen infection.

A BAC-Based Physical Map of Zhikong Scallop (Chlamys farreri Jones et Preston)

Zhikong scallop (Chlamys farreri) is one of the most economically important aquaculture species in China. Physical maps are crucial tools for genome sequencing, gene mapping and cloning, genetic improvement and selective breeding. In this study, we have developed a genome-wide, BAC-based physical map for the species. A total of 81,408 clones from two BAC libraries of the scallop were fingerprinted using an ABI 3130xl Genetic Analyzer and a fingerprinting kit developed in our laboratory. After data processing, 63,641 (∼5.8× genome coverage) fingerprints were validated and used in the physical map assembly. A total of 3,696 contigs were assembled for the physical map. Each contig contained an average of 10.0 clones, with an average physical size of 490 kb. The combined total physical size of all contigs was 1.81 Gb, equivalent to approximately 1.5 fold of the scallop haploid genome. A total of 10,587 BAC end sequences (BESs) and 167 markers were integrated into the physical map. We evaluated the physical map by overgo hybridization, BAC-FISH (fluorescence in situ hybridization), contig BAC pool screening and source BAC library screening. The results have provided evidence of the high reliability of the contig physical map. This is the first physical map in mollusc; therefore, it provides an important platform for advanced research of genomics and genetics, and mapping of genes and QTL of economical importance, thus facilitating the genetic improvement and selective breeding of the scallop and other marine molluscs.

Horizontally transferred genes in the genome of Pacific white shrimp, Litopenaeus vannamei.

BackgroundIn recent years, as the development of next-generation sequencing technology, a growing number of genes have been reported as being horizontally transferred from prokaryotes to eukaryotes, most of them involving arthropods. As a member of the phylum Arthropoda, the Pacific white shrimp Litopenaeus vannamei has to adapt to the complex water environments with various symbiotic or parasitic microorganisms, which provide a platform for horizontal gene transfer (HGT).ResultsIn this study, we analyzed the genome-wide HGT events in L. vannamei. Through homology search and phylogenetic analysis, followed by experimental PCR confirmation, 14 genes with HGT event were identified: 12 of them were transferred from bacteria and two from fungi. Structure analysis of these genes showed that the introns of the two fungi-originated genes were substituted by shrimp DNA fragment, two genes transferred from bacteria had shrimp specific introns inserted in them. Furthermore, around other three bacteria-originated genes, there were three large DNA segments inserted into the shrimp genome. One segment was a transposon that fully transferred, and the other two segments contained only coding regions of bacteria. Functional prediction of these 14 genes showed that 6 of them might be related to energy metabolism, and 4 others related to defense of the organism.ConclusionsHGT events from bacteria or fungi were happened in the genome of L. vannamei, and these horizontally transferred genes can be transcribed in shrimp. This is the first time to report the existence of horizontally transferred genes in shrimp. Importantly, most of these genes are exposed to a negative selection pressure and appeared to be functional.

The sea cucumber genome provides insights into morphological evolution and visceral regeneration

Apart from sharing common ancestry with chordates, sea cucumbers exhibit a unique morphology and exceptional regenerative capacity. Here we present the complete genome sequence of an economically important sea cucumber, A. japonicus, generated using Illumina and PacBio platforms, to achieve an assembly of approximately 805 Mb (contig N50 of 190 Kb and scaffold N50 of 486 Kb), with 30,350 protein-coding genes and high continuity. We used this resource to explore key genetic mechanisms behind the unique biological characters of sea cucumbers. Phylogenetic and comparative genomic analyses revealed the presence of marker genes associated with notochord and gill slits, suggesting that these chordate features were present in ancestral echinoderms. The unique shape and weak mineralization of the sea cucumber adult body were also preliminarily explained by the contraction of biomineralization genes. Genome, transcriptome, and proteome analyses of organ regrowth after induced evisceration provided insight into the molecular underpinnings of visceral regeneration, including a specific tandem-duplicated prostatic secretory protein of 94 amino acids (PSP94)-like gene family and a significantly expanded fibrinogen-related protein (FREP) gene family. This high-quality genome resource will provide a useful framework for future research into biological processes and evolution in deuterostomes, including remarkable regenerative abilities that could have medical applications. Moreover, the multiomics data will be of prime value for commercial sea cucumber breeding programs.

Virus-derived small RNAs in the penaeid shrimp Fenneropenaeus chinensis during acute infection of the DNA virus WSSV

Small interfering RNAs (siRNAs) and microRNAs (miRNAs) are two classes of small RNAs (sRNAs) that are critical for virus-host interplay via the RNA interference (RNAi) pathway. One virus-derived siRNA and numerous miRNAs has been reported for the double-stranded DNA virus white spot syndrome virus (WSSV), however, the expression profiles of these different types of sRNAs have not been assessed. Here, by sequencing the sRNAs and mRNAs of WSSV-infected Chinese shrimp (Fenneropenaeus chinensis), we found that the viral transcripts were universally targeted by WSSV-derived siRNAs, supporting a pivotal role for RNAi in the anti-viral immunity of shrimp. The genesis of WSSV-derived siRNAs was associated with long RNA structures. Moreover, by separating miRNAs from siRNAs, 12 WSSV miRNAs were identified. Investigation of conserved viral miRNA targets in different host species indicated the involvement of viral miRNAs in host immune responses. Collectively, our data provide new insights into the role of the RNAi pathway in the interplay between DNA viruses and crustaceans.

Genome Sequences of Marine Shrimp Exopalaemon carinicauda Holthuis Provide Insights into Genome Size Evolution of Caridea

Crustacea, particularly Decapoda, contains many economically important species, such as shrimps and crabs. Crustaceans exhibit enormous (nearly 500-fold) variability in genome size. However, limited genome resources are available for investigating these species. Exopalaemon carinicauda Holthuis, an economical caridean shrimp, is a potential ideal experimental animal for research on crustaceans. In this study, we performed low-coverage sequencing and de novo assembly of the E. carinicauda genome. The assembly covers more than 95% of coding regions. E. carinicauda possesses a large complex genome (5.73 Gb), with size twice higher than those of many decapod shrimps. As such, comparative genomic analyses were implied to investigate factors affecting genome size evolution of decapods. However, clues associated with genome duplication were not identified, and few horizontally transferred sequences were detected. Ultimately, the burst of transposable elements, especially retrotransposons, was determined as the major factor influencing genome expansion. A total of 2 Gb repeats were identified, and RTE-BovB, Jockey, Gypsy, and DIRS were the four major retrotransposons that significantly expanded. Both recent (Jockey and Gypsy) and ancestral (DIRS) originated retrotransposons responsible for the genome evolution. The E. carinicauda genome also exhibited potential for the genomic and experimental research of shrimps.

EST-derived SNP discovery and selective pressure analysis in Pacific white shrimp (Litopenaeus vannamei)

Pacific white shrimp has become a major aquaculture and fishery species worldwide. Although a large scale EST resource has been publicly available since 2008, the data have not yet been widely used for SNP discovery or transcriptome-wide assessment of selective pressure. In this study, a set of 155 411 expressed sequence tags (ESTs) from the NCBI database were computationally analyzed and 17 225 single nucleotide polymorphisms (SNPs) were predicted, including 9 546 transitions, 5 124 transversions and 2 481 indels. Among the 7 298 SNP substitutions located in functionally annotated contigs, 58.4% (4 262) are non-synonymous SNPs capable of introducing amino acid mutations. Two hundred and fifty nonsynonymous SNPs in genes associated with economic traits have been identified as candidates for markers in selective breeding. Diversity estimates among the synonymous nucleotides were on average 3.49 times greater than those in non-synonymous, suggesting negative selection. Distribution of non-synonymous to synonymous substitutions (Ka/Ks) ratio ranges from 0 to 4.01, (average 0.42, median 0.26), suggesting that the majority of the affected genes are under purifying selection. Enrichment analysis identified multiple gene ontology categories under positive or negative selection. Categories involved in innate immune response and male gamete generation are rich in positively selected genes, which is similar to reports in Drosophila and primates. This work is the first transcriptome-wide assessment of selective pressure in a Penaeid shrimp species. The functionally annotated SNPs provide a valuable resource of potential molecular markers for selective breeding.

BAC end sequencing of Pacific white shrimp Litopenaeus vannamei : a glimpse into the genome of Penaeid shrimp

Little is known about the genome of Pacific white shrimp (Litopenaeus vannamei). To address this, we conducted BAC (bacterial artificial chromosome) end sequencing of L. vannamei. We selected and sequenced 7 812 BAC clones from the BAC library LvHE from the two ends of the inserts by Sanger sequencing. After trimming and quality filtering, 11 279 BAC end sequences (BESs) including 4 609 pairedends BESs were obtained. The total length of the BESs was 4 340 753 bp, representing 0.18% of the L. vannamei haploid genome. The lengths of the BESs ranged from 100 bp to 660 bp with an average length of 385 bp. Analysis of the BESs indicated that the L. vannamei genome is AT-rich and that the primary repeats patterns were simple sequence repeats (SSRs) and low complexity sequences. Dinucleotide and hexanucleotide repeats were the most common SSR types in the BESs. The most abundant transposable element was gypsy, which may contribute to the generation of the large genome size of L. vannamei. We successfully annotated 4 519 BESs by BLAST searching, including genes involved in immunity and sex determination. Our results provide an important resource for functional gene studies, map construction and integration, and complete genome assembly for this species.