来自易变山羊草抗禾谷孢囊线虫基因的克隆与序列分析

禾谷孢囊线虫（Heterodera avenae）是严重危害禾谷类作物的病原线虫之一,它广泛分布于澳大利亚、欧洲、北美、印度和中国等世界主要小麦产区,使作物严重减产,造成巨大的经济损失。目前最有效的防治措施之一是将外源抗性基因导入栽培小麦(Triticum aestivum L.)，培育抗禾谷孢囊线虫的新品种。但迄今为止抗禾谷孢囊线虫基因克隆研究的相关报道却很少。 本实验根据此前从抗禾谷孢囊线虫材料E-10扩增得到的与来自节节麦（Aegilops tauschii）的抗禾谷孢囊线虫基因Cre3高度同源的序列Rccn4,设计出三条嵌套引物,采用SON-PCR(single oligonucleotide nested PCR)方法,从E-10基因组DNA中得到一个长为1264 bp的扩增产物(命名为Rccn-L),测序比对结果显示,这一序列将Rccn4的3’端延伸了1209 bp,与抗禾谷孢囊线虫Cre3基因核苷酸同源性为86﹪，核苷酸编码区长1026 bp，含一个不完整的开放阅读框,一个终止密码子，没有起始密码子和内含子结构,编码一个342个氨基酸残基的蛋白质。该蛋白质等电点为5.19，分子量为38112.6Da。从序列的第113位开始到第332位是NBS-LRR类抗病性基因LRR区,呈现XXLXXLXXL重复。LRR编码区内亮氨酸残基的含量达17﹪，与抗禾谷孢囊线虫Cre3基因LRR编码区的核苷酸和氨基酸同源性分别为89﹪和78﹪。本实验首次将SON-PCR成功地运用于植物基因克隆，为植物基因克隆提供了又一有效方法。 此外,还根据Cre3基因及其他的NBS-LRR类植物抗性基因的NBS和LRR区保守序列设计了两对特异性引物，从禾谷孢囊线虫抗性材料易变山羊草基因组DNA中扩增到两个相应的目标条带。测序分析结果表明，它们的长度分别为532bp和1175bp，构成了一个有32bp的共同序列的NBS-LRR编码区。其序列总长为1675bp（命名为RCCN），含有一个不完整的开放阅读框，没有起始密码子、终止密码子和内含子结构。其中编码序列为1673bp，可编码一个557个氨基酸的蛋白质，等电点（pI）为5.39，分子量为63537.5Da。与Cre3的核苷酸和氨基酸同源性分别为87.8﹪和77﹪。RCCN氨基酸序列中含有已知抗病基因NBS区域的几个保守模体：kinase2区的ILDD、kinase3的（ⅰ）ESKILVTTRSK，（ⅱ）KGSPLAARTVGG，（ⅲ）RRCFAYCS及EGF。RCCN NBS区与Cre3 NBS区的核苷酸和氨基酸的同源性分别为96.4﹪和94﹪。从氨基酸序列的274位到548位为LRR保守区，呈现不规则的aXXLXXLXXL（其中a代表I，V，L，F或M）重复，其中亮氨酸的含量为15.6﹪。该区域与Cre3的LRR区的核苷酸和氨基酸同源性分别为80.8﹪和74﹪。推测该序列可能为一个抗禾谷孢囊线虫的新基因。 本文对抗禾谷孢囊线虫基因的克隆研究，为进一步克隆基因全序列，探索其结构与功能，和研究该基因表达与调控提供了关键信息。同时也为通过基因工程途径将抗性基因向优良小麦品种高效、定向转移，最终培育出小麦抗禾谷孢囊线虫新品种奠定了基础。 Cereal cyst nematode (CCN) is a damaging pathogen of broad acre cereal crops in Australia, Europe, North America, India and China. It affects wheat, barley, oat and triticale and causes yield loss of up to 80%. At present, Transferring resistance genes against CCN into wheat cultivars and breeding varieties are considered one of the most effective methods for controlling the CCN. However, there are very limited reports concerning the cloning studies of resistance genes against the cereal cyst nematode. According to the sequence of Rccn4 which had high similarity to the nucleotide binding site (NBS) coding region of cereal cyst nematode resistance gene, Cre3， We designed three 3’ nested primers. Using single oligonucleotide nested PCR (SON-PCR) we successfully amplified one band, Rccn-L, of 1264bp from E-10 which is the wheat-Ae.variabilis translocation line containing the cereal cyst nematode resistance gene of Ae.variabilis. We found that this band of interesting is the 3’ flanking sequence of 1209bp in size of Rccn4. The coding region was 1026bp, which contained an incomplete open reading frame and a terminator codon, without initiation codon and intron, encoding a peptide of 342 amino acid residues, and shared 86﹪nucleotide sequence identity with Cre3. This peptide had a conserved LRR domain, containing the imperfect repeats,XXLXXLXXL, which contains 17﹪ leucine residues and shares, respectively, 89﹪ nucleotide sequence and 78﹪ amino acid sequence identity with the LRR sequence of Cre3 locus. This research firstly used SON-PCR in the research of plant genome successfully, which indicated that SON-PCR is another method of cloning plant gene. At the same time, According to the conversed motif of NBS and LRR region of cereal cyst nematode resistance gene Cre3 from wild wheat (Triticum tauschlii L.) and the known NBS-LRR group resistance genes, we designed two pairs of specific primers for NBS and LRR region respectively. One band of approximately 530bp was amplified using the specific primers for conversed NBS region and one band of approximately 1200bp was amplified with the specific primers for conversed LRR region. After sequencing, we found that these two sequences included 32bp common nucleotide sequence and have 1675 bp in total, which was registered as RCCN in the Genbank. RCCN contained a NBS-LRR domain and an incomplete open reading frame without initiation codon, terminator codon and inxon. Its exon encodes a peptide of 557 amino acid residues. The molecular weight of the protein from the amino acid was 63.537 KDa. The amino acid sequence of RCCN contained conserved motif: ILDD, ESKILVTTRSK, KGSPLAARTVGG, RRCFAYCS, EGF，LRR. RCCN shares 87.8﹪ nucleotide sequence and 77﹪ amino acid sequence identity with cereal cyst nematode gene Cre3. It might be a novel cereal cyst nematode resistance gene. These research results of cloning the resistance genes against cereal cyst nematode bring a great promise for transferring resistance genes into wheat cultivars and breeding new wheat varieties against cereal cyst nematode by gene engineering. And these results also lay the hard foundation for the expressing researches of these genes.

小麦抗禾谷孢囊线虫（ Heterodera avenae）相关基因的研究

禾谷孢囊线虫严重影响禾谷类作物的产量，在小麦中由禾谷孢囊线虫引起的产量损失可达30-100%。尤其在澳大利亚、欧洲、印度和中东危害严重，目前禾谷孢囊线虫已成为危害我国作物的主要病源。控制禾谷孢囊线虫的方法主要有：作物轮作、杀线虫剂、寄主抗性等等，其中基因工程方法培育抗线虫小麦品种被认为是最经济有效的方法。分离抗禾谷类孢囊线虫基因对揭示抗性基因结构与功能及其表达调控具有重要意义。 尽管小麦是重要的粮食作物，在小麦中已发现的抗禾谷孢囊线虫的基因很少，而比其近缘属如节节麦、易变山羊草、偏凸山羊草中含有丰富的抗源。目前已鉴定出禾谷孢囊线虫抗性位点Cre，并发现了9个禾谷孢囊线虫抗性基因(Cre1,2, 3, 4, 5, 6, 7, 8, and R) ，其中只有Cre1和Cre8直接从普通小麦中获得。从节节麦中获得的Cre3基因能最有效的控制线虫数量，其次是Cre1和Cre8。这些基因的克隆对于了解禾谷孢囊线虫抗性机制及进一步的育种应用都是非常关键的。然而，目前为止仅有Cre3基因通过图位克隆的方法从节节麦中被分离得到。该基因已被克隆得到的多数线虫抗性基因一样均属于核苷酸结合位点区（NBS）－亮氨酸重复序列区（LRR）基因家族。目前，已有很多抗性基因被分离，这些已知的NBS-LRR类抗性基因的保守序列为应用PCR的方法克隆新的抗性基因提供了可能。 因此本课题的目的是采用保守区同源克隆、3′RACE 和5′RACE 等方法从抗禾谷孢囊线虫小麦-易变山羊草小片段易位系E10 中克隆小麦抗禾谷孢囊线虫基因全序列，进而通过半定量PCR 和荧光定量PCR 研究该基因的表达模式。同时通过mRNA 差别显示技术和任意引物PCR（RAP-PCR）技术分离克隆植物禾谷孢囊线虫抗性基因及其相关基因，为阐明植物抗病性分子机制以及改良作物抗病性和作物育种提供基础，为通过分子标记辅助育种和基因工程方法实现高效、定向转移抗病基因到优良小麦品种奠定了重要的理论和物质基础。主要研究结果： 1． 本实验根据此前从抗禾谷孢囊线虫材料E-10 扩增得到的与来自节节麦的抗禾谷孢囊线虫Cre3 基因及其他的NBS-LRR 类抗性基因的NBS 和LRR 保守区序列设计了两对特异性引物，从E10 中扩增到532bp 和1175bp 的两个目标条带，它们有一个32bp 的共同序列，连接构成总长为1675bp 的NBS-LRR 编码区（命名为RCCN）。根据RCCN设计引物，利用NBS-LRR区序列设计引物，通过5′RACE 和3′RACE 技术采用3′-Full RACE Core Set(TaKaRa)和5'-Full RACE Kit (TaKaRa)试剂盒，反转录后通过嵌套引物GSP1 和GSP2 分别进行两轮基因特异性扩增，分别将NBS_LRR 区向5′端和3′端延伸了1173bp 和449bp，并包含了起始密码子和终止密码子。根据拼接的得到的序列重新设计引物扩增进行全基因扩增的结果与上面获得的一致。拼接后得到全长2775 bp 的基因序列（记作CreZ, GenBank 号：EU327996）。CreZ 基因包括完整的开放阅读框，全长2775 bp，编码924个氨基酸。序列分析表明它与已知的禾谷孢囊线虫抗性基因Cre3的一致性很高，并且它与已经报到的NBS-LRR 类疾病抗性基因有着相同的保守结构域。推测CreZ基因可能是一个新的NBS-LRR 类禾谷孢囊线虫抗性基因，该基因的获得为通过基因工程途径培育抗禾谷孢囊线虫小麦新品种奠定了基础，并为抗禾谷孢囊线虫基因的调控表达研究提供了参考。 2． 通过半定量PCR和SYBR Green荧光定量PCR技术对CreZ基因的相对表达模式进行了研究。以α-tubulin 2作为参照，采用半定量PCR 分析CreZ 基因在不同接种时期1d, 5d, 10, 15d 的E-10的根和叶的的表达情况。在内参扩增一致的条件下，CreZ 在E-10的根部随着侵染时间的增加表达量有明显的增加，在没有侵染的E-10的根部其表达量没有明显变化，而在叶中没有检测表达，说明该基因只在抗性材料的根部表达。SYBR Green定量PCR分析接种前后E10根部基因CreZ基因的表达水平为检测CreZ基因的表达建立了一套灵敏、可靠的SYBRGreen I 荧光定量PCR 检测方法。接种禾谷孢囊线虫后E10根内CreZ基因的相对表达水平显著高于接种前。随接种时间的延长持续增加，最终CreZ基因的相对表达量达到未接种的对照植株的10.95倍。小麦禾谷孢囊线虫抗性基因CreZ的表达量与胁迫呈正相关，表明其与小麦的的禾谷孢囊线虫抗性密切相关，推测CreZ基因可能是一个新的禾谷孢囊线虫候选抗性基因。 3． 针对小麦基因组庞大、重复序列较多，禾谷孢囊线虫抗性基因及其相关基因的片断难以有效克隆的问题，通过mRNA 差别显示技术及RAP-PCR 技术分离克隆植物禾谷孢囊线虫抗性及其相关基因。试验最终得到154 条差异表达条带，将回收得到的差异条带的二次PCR 扩增产物经纯化后点到带正电的尼龙膜上，进行反向Northern 杂交筛选，最终筛选得到102 个阳性差异点。将其中81 个进行测序，并将序列提交到Genbank 中的dbEST 数据库，分别获得登录号(FE192210 -FE192265，FE193048- FE193074 )。序列比对分析发现，其中26 个序列与已知功能的基因序列同源；有28 条EST 序列在已有核酸数据库中未找到同源已知基因和EST，属新的ESTs 序列；另外27 个EST 序列与已知核酸数据库中的ESTs 具有一定相似性，但功能未知。其所得ESTs 序列补充了Genbank ESTs 数据库，为今后进一步开展抗禾谷类孢囊线虫基因研究工作打下了基础。结合本试验功能基因的相关信息，对小麦接种禾谷孢囊线虫后产生的抗性机制进行了探讨。接种禾谷孢囊线虫后植物在mRNA 水平上的应答是相当复杂的，同时植物的抗病机制是一个复杂的过程，涉及到多个代谢途径的相互作用。 The cereal cyst nematode (CCN), Heterodera avenae Woll, causes severe yieldreductions in cereal crops. The losses caused by CCN can be up to 30-100% in somewheat fields. At present, cereal cyst nematode has become the major disease sourcein China and it also damaged heavily in Australia, Europe, India and Middle East.The damage caused by CCN can be mitigated through several methods, includingcrop rotation, nematicide application, cultural practice, host resistance, and others.Of these methods, incorporating resistance genes into wheat cultivars and breedingresistant lines is considered to be the most cost-effective control measure forreducing nematode populations. Although wheat is an economically important crop around the world, far fewergenes resistant to CCN were found in wheat than were detected in its relatives, suchas Aegilops taucchi, Aegilops variabilis and Aegilops ventricosa. Cloning these genesis essential for understanding the mechanism of this resistance and for furtherapplication in breeding. Because of the huge genome and high repeat sequencescontent, the efficient methods to clone genes from cereal crops, are still lacking. A resistance locus, Cre, has been identified and 9 genes resistant to CCN (designatedCre1, 2, 3, 4, 5, 6, 7, 8, and R) have been described, in which Cre1 and Cre8 werederived directly from common wheat. The Cre3 locus, which was derived from Ae.tauschii, has the greatest impact on reducing the number of female cysts, followed byCre1 and Cre8. Cloning these genes is essential for understanding the mechanism ofthis resistance and for further application in breeding. However, to this point, only Cre3, a NBS-LRR disease resistance gene, has been obtained through mappingcloning in Ae. tauschii. The majority of nematode resistance genes cloned so far belong to a super familywhich contains highly conserved nucleotide-binding sites (NBS) and leucine-richrepeat (LRR) domains. To date, many NBS-LRR resistance genes have been isolated.The conserved sequences of these recognized NBS-LRR resistance genes provide thepossibility to isolate novel resistance genes using a PCR-based strategy. The aim of the present study was to clone the resistance gene of CCN fromWheat/Aegilops variabilis small fragment chromosome translocation line E10 whichis resistant to CCN and investigate the espression profiles of this gene withsemi-quantitative PCR and real-time PCR. Another purpose of this study is cloningthe relational resistance gene for CCN by mRNA differential display PCR andRAP-PCR. These works will offer a foundation for disease defence of crop andbreeding and directional transferring resistance gene into wheat with geneengineering. Primary results as following: 1．According to the conversed motif of NBS and LRR region of cereal cystnematode resistance gene Cre3 from wild wheat (Triticum tauschlii) and the knownNBS-LRR group resistance genes, we designed two pairs of specific primers for NBSand LRR region respectively. One band of approximately 530bp was amplified usingthe specific primers for conversed NBS region and one band of approximately 1175bpwas amplified with the specific primers for conversed LRR region. After sequencing,we found that these two sequences included 32bp common nucleotide having 1675bpin total, which was registered as RCCN in the Genbank. Based on the conservedregions of known resistance genes, a NBS-LRR type CCN resistance gene analog wasisolated from the CCN resistant line E-10 of the wheat near isogenic lines (NILs), by5′RACE and 3′ RACE.designated as CreZ (GenBank accession number: EU327996) .It contained a comlete ORF of 2775 bp and encoded 924 amino acids. Sequencecomparison indicated that it shared 92% nucleotide and 87% amino acid identitieswith those of the known CCN-resistance gene Cre3 and it had the same characteristic of the conserved motifs as other established NBS-LRR disease resistance genes. 2． Usingα-tubulin 2 as exoteric reference, semi-quantitative PCR and real-timePCR analysis were conducted. The expression profiling of CreZ indicated that it wasspecifically expressed in the roots of resistant plants and its relative expression levelincreased sharply when the plants were inoculated with cereal cyst nematodes. therelative expression level of the 15days-infected E10 is the 10.95 times as that ofuninfected E10,ultimately. It was inferred that the CreZ gene be a novel potentialresistance gene to CCN. 3．We cloned the relational resistance gene for CCN by mRNA differentialdisplay PCR and arbitrarily primed PCR fingerprinting of RNA from wheat whichpossess huge and high repeat sequence content genomes. Total 154 differentialexpression bands were separated and second amplified by PCR. The products werenylon membrane. The 102 positive clones were filtrated by reverse northern dot blotand 81 of those were sent to sequence. The EST sequences were submitted toGenbank (Genbank accession: FE192210 - FE192265, FE193048 - FE193074). Thesequences alignment analysis indicated 26 of them were identical with known genes;28 were not found identical sequence in nucleic acid database; another 27 ests wereidentical with some known ests, but their functions were not clear. These ESTsenriched Genbank ESTs database and offered foundation for further research ofresistance gene of CCN.

Atomic force microscope investigation of large-circle DNA molecules

A circular bacterial artificial chromosome of 148.9 kbp on human chromosome 3 has been extended and fixed on bare mica substrates using a developed fluid capillary flow method in evaporating liquid drops. Extended circular DNA molecules were imaged with an atomic force microscope (AFM) under ambient conditions. The measured total lengths of the whole DNA molecules were in agreement with sequencing analysis data with an error range of +/-3.6%. This work is important groundwork for probing single nucleotide polymorphisms in the human genome, mapping genomic DNA, manipulating biomolecular nanotechnology, and studying the interaction of DNA-protein complexes investigated by AFM.

The complete mitochondrial genome of the large yellow croaker, Larimichthys crocea (Perciformes, Sciaenidae): Unusual features of its control region and the phylogenetic position of the Sciaenidae

To understand the systematic status of Larimichthys crocea in the Percoidei, we determined the complete mitochondrial (mt) genome sequence using 454 sequencing-by-synthesis technology. The complete mt genome is 16,466 bp in length including the typical structure of 22 tRNAs, 2 rRNAs, 13 protein-coding genes and the noncoding control region (CR). Further sequencing for the complete CR was performed using the primers Cyt b-F and 12S-R on six L crocea individuals and two L polyactis individuals. Interestingly, all seven CR sequences from L crocea were identical while the three sequences from L polyactis were distinct (including one from GenBank). Although the conserved blocks such as TAS and CSB-1, -2, and -3 are readily identifiable in the control regions of the two species, the typical central conserved blocks CSB-D, -E, and -F could not be detected, while they are found in Cynoscion acoupa of Sciaenidae and other Percoidei species. Phylogenetic analysis shows that L crocea is a relatively recently emerged species in Sciaenidae and this family is closely related to family Pomacanthidae within the Percoidei. L crocea, as the first species of Sciaenidae with complete mitochondrial genome available, will provide important information on the molecular evolution of the group. Moreover, the genus-specific pair of primers designed in this study for amplifying the complete mt control region will be very useful in studies on the population genetics and conservation biology of Larimichthys. (c) 2008 Elsevier B.V. All rights reserved.

The complete mitochondrial genome sequence of the cutlassfish Trichiurus japonicus (Perciformes: Trichiuridae) : Genome characterization and phylogenetic considerations

Mitochondrial genome sequence and structure analysis has become a powerful tool for studying molecular evolution and phylogenetic relationships. To understand the systematic status of Trichiurus japonicus in suborder Scombroidei, we determined the complete mitochondrial genome (mitogenome) sequence using the long-polymerase chain reaction (long-PCR) and shotgun sequencing method. The entire mitogenome is 16,796 by in length and has three unusual features, including (1) the absence of tRNA(Pro) gene, (2) the possibly nonfunctional light-strand replication origin (O-L) showing a shorter loop in secondary structure and no conserved motif (5'-GCCGG-3'), (3) two sets of the tandem repeats at the 5' and 3' ends of the control region. The three features seem common for Trichiurus mitogenomes, as we have confirmed them in other three T. japonicus individuals and in T nanhaiensis. Phylogenetic analysis does not support the monophyly of Trichiuridae, which is against the morphological result. T. japonicus is most closely related to those species of family Scombridae; they in turn have a sister relationship with Perciformes members including suborders Acanthuroidei, Caproidei, Notothenioidei, Zoarcoidei, Trachinoidei, and some species of Labroidei, based on the current dataset of complete mitogenome. T japonicus together with T. brevis, T lepturus and Aphanopus carbo form a clade distinct from Lepidopus caudatus in terms of the complete Cyt b sequences. T. japonicus mitogenome, as the first discovered complete mitogenome of Trichiuridae, should provide important information on both genomics and phylogenetics of Trichiuridae. (C) 2009 Elsevier B.V. All rights reserved.

Analysis of tandem repeats in the genome of Chinese shrimp Fenneropenaeus chinensis

Through random sequencing, we found a total of 884000 base-pairs (bp) of random genomic sequences in the genome of Chinese shrimp (Fenneropenaeus chinensis). Using bio-soft Tandem Repeat Finder (TRF) software, 2159 tandem repeats were found, in which there were 1714 microsatellites and 445 minisatellites, accounting for 79.4% and 20.6% of repeat sequences, respectively. The cumulative length of repeat sequences was found to be 116685 bp, accounting for 13.2% of the total DNA sequence; the cumulative length of microsatellites occupied 9.78% of the total DNA sequence, and that of minisatellites occupied 3.42%. In decreasing order, the 20 most abundant repeat sequence classes were as follows: AT (557), AC (471), AG (274), AAT (92), A (56), AAG (28), ATC (27), ATAG (27), AGG (18), ACT (15), C (11), AAC (11), ACAT (11), CAGA (10), AGAA (9), AGGG (7), CAAA (7), CGCA (6), ATAA (6), AGAGAA (6). Dinucleotide repeats, not only in the aspect of the number, but also in cumulative length, were the preponderant repeat type. There were few classes and low copy numbers of repeat units of the pentanucleotide repeat type, which included only three classes: AGAGA, GAGGC and AAAGA. The classes and copy numbers of heptanucleotide, eleven-nucleotide and thirteen-nucleotide primer-number-composed repeats were distinctly less than that of repeat types beside them.

Complete mitochondrial genome of the Asian paddle crab Charybdis japonica (Crustacea: Decapoda: Portunidae): gene rearrangement of the marine brachyurans and phylogenetic considerations of the decapods

Given the commercial and ecological importance of the Asian paddle crab, Charybdis japonica, there is a clearly need for genetic and molecular research on this species. Here, we present the complete mitochondrial genome sequence of C. japonica, determined by the long-polymerase chain reaction and primer walking sequencing method. The entire genome is 15,738 bp in length, encoding a standard set of 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes, plus the putative control region, which is typical for metazoans. The total A+T content of the genome is 69.2%, lower than the other brachyuran crabs except for Callinectes sapidus. The gene order is identical to the published marine brachyurans and differs from the ancestral pancrustacean order by only the position of the tRNA (His) gene. Phylogenetic analyses using the concatenated nucleotide and amino acid sequences of 13 protein-coding genes strongly support the monophyly of Dendrobranchiata and Pleocyemata, which is consistent with the previous taxonomic classification. However, the systematic status of Charybdis within subfamily Thalamitinae of family Portunidae is not supported. C. japonica, as the first species of Charybdis with complete mitochondrial genome available, will provide important information on both genomics and molecular ecology of the group.

The ecological effects of exotic disease resistance genes introgressed into British gooseberries

Warren, J. and James, P. (2006). The ecological effects of exotic disease resistance genes introgressed into British gooseberries. Oecologia 147(1),69-75. RAE2008

Transmission of single HIV-1 genomes and dynamics of early immune escape revealed by ultra-deep sequencing.

We used ultra-deep sequencing to obtain tens of thousands of HIV-1 sequences from regions targeted by CD8+ T lymphocytes from longitudinal samples from three acutely infected subjects, and modeled viral evolution during the critical first weeks of infection. Previous studies suggested that a single virus established productive infection, but these conclusions were tempered because of limited sampling; now, we have greatly increased our confidence in this observation through modeling the observed earliest sample diversity based on vastly more extensive sampling. Conventional sequencing of HIV-1 from acute/early infection has shown different patterns of escape at different epitopes; we investigated the earliest escapes in exquisite detail. Over 3-6 weeks, ultradeep sequencing revealed that the virus explored an extraordinary array of potential escape routes in the process of evading the earliest CD8 T-lymphocyte responses--using 454 sequencing, we identified over 50 variant forms of each targeted epitope during early immune escape, while only 2-7 variants were detected in the same samples via conventional sequencing. In contrast to the diversity seen within epitopes, non-epitope regions, including the Envelope V3 region, which was sequenced as a control in each subject, displayed very low levels of variation. In early infection, in the regions sequenced, the consensus forms did not have a fitness advantage large enough to trigger reversion to consensus amino acids in the absence of immune pressure. In one subject, a genetic bottleneck was observed, with extensive diversity at the second time point narrowing to two dominant escape forms by the third time point, all within two months of infection. Traces of immune escape were observed in the earliest samples, suggesting that immune pressure is present and effective earlier than previously reported; quantifying the loss rate of the founder virus suggests a direct role for CD8 T-lymphocyte responses in viral containment after peak viremia. Dramatic shifts in the frequencies of epitope variants during the first weeks of infection revealed a complex interplay between viral fitness and immune escape.

Sequencing three crocodilian genomes to illuminate the evolution of archosaurs and amniotes.

The International Crocodilian Genomes Working Group (ICGWG) will sequence and assemble the American alligator (Alligator mississippiensis), saltwater crocodile (Crocodylus porosus) and Indian gharial (Gavialis gangeticus) genomes. The status of these projects and our planned analyses are described.

Conversion of natural forests to managed forest plantations decreases tree resistance to prolonged droughts

© 2015 Published by Elsevier B.V.Throughout the southern US, past forest management practices have replaced large areas of native forests with loblolly pine plantations and have resulted in changes in forest response to extreme weather conditions. However, uncertainty remains about the response of planted versus natural species to drought across the geographical range of these forests. Taking advantage of a cluster of unmanaged stands (85-130year-old hardwoods) and managed plantations (17-20year-old loblolly pine) in coastal and Piedmont areas of North Carolina, tree water use, cavitation resistance, whole-tree hydraulic (Ktree) and stomatal (Gs) conductances were measured in four sites covering representative forests growing in the region. We also used a hydraulic model to predict the resilience of those sites to extreme soil drying. Our objectives were to determine: (1) if Ktree and stomatal regulation in response to atmospheric and soil droughts differ between species and sites; (2) how ecosystem type, through tree water use, resistance to cavitation and rooting profiles, affects the water uptake limit that can be reached under drought; and (3) the influence of stand species composition on critical transpiration that sets a functional water uptake limit under drought conditions. The results show that across sites, water stress affected the coordination between Ktree and Gs. As soil water content dropped below 20% relative extractable water, Ktree declined faster and thus explained the decrease in Gs and in its sensitivity to vapor pressure deficit. Compared to branches, the capability of roots to resist high xylem tension has a great impact on tree-level water use and ultimately had important implications for pine plantations resistance to future summer droughts. Model simulations revealed that the decline in Ktree due to xylem cavitation aggravated the effects of soil drying on tree transpiration. The critical transpiration rate (Ecrit), which corresponds to the maximum rate at which transpiration begins to level off to prevent irreversible hydraulic failure, was higher in managed forest plantations than in their unmanaged counterparts. However, even with this higher Ecrit, the pine plantations operated very close to their critical leaf water potentials (i.e. to their permissible water potentials without total hydraulic failure), suggesting that intensively managed plantations are more drought-sensitive and can withstand less severe drought than natural forests.

Computational Inference of Genome-Wide Protein-DNA Interactions Using High-Throughput Genomic Data

Transcriptional regulation has been studied intensively in recent decades. One important aspect of this regulation is the interaction between regulatory proteins, such as transcription factors (TF) and nucleosomes, and the genome. Different high-throughput techniques have been invented to map these interactions genome-wide, including ChIP-based methods (ChIP-chip, ChIP-seq, etc.), nuclease digestion methods (DNase-seq, MNase-seq, etc.), and others. However, a single experimental technique often only provides partial and noisy information about the whole picture of protein-DNA interactions. Therefore, the overarching goal of this dissertation is to provide computational developments for jointly modeling different experimental datasets to achieve a holistic inference on the protein-DNA interaction landscape.

We first present a computational framework that can incorporate the protein binding information in MNase-seq data into a thermodynamic model of protein-DNA interaction. We use a correlation-based objective function to model the MNase-seq data and a Markov chain Monte Carlo method to maximize the function. Our results show that the inferred protein-DNA interaction landscape is concordant with the MNase-seq data and provides a mechanistic explanation for the experimentally collected MNase-seq fragments. Our framework is flexible and can easily incorporate other data sources. To demonstrate this flexibility, we use prior distributions to integrate experimentally measured protein concentrations.

We also study the ability of DNase-seq data to position nucleosomes. Traditionally, DNase-seq has only been widely used to identify DNase hypersensitive sites, which tend to be open chromatin regulatory regions devoid of nucleosomes. We reveal for the first time that DNase-seq datasets also contain substantial information about nucleosome translational positioning, and that existing DNase-seq data can be used to infer nucleosome positions with high accuracy. We develop a Bayes-factor-based nucleosome scoring method to position nucleosomes using DNase-seq data. Our approach utilizes several effective strategies to extract nucleosome positioning signals from the noisy DNase-seq data, including jointly modeling data points across the nucleosome body and explicitly modeling the quadratic and oscillatory DNase I digestion pattern on nucleosomes. We show that our DNase-seq-based nucleosome map is highly consistent with previous high-resolution maps. We also show that the oscillatory DNase I digestion pattern is useful in revealing the nucleosome rotational context around TF binding sites.

Finally, we present a state-space model (SSM) for jointly modeling different kinds of genomic data to provide an accurate view of the protein-DNA interaction landscape. We also provide an efficient expectation-maximization algorithm to learn model parameters from data. We first show in simulation studies that the SSM can effectively recover underlying true protein binding configurations. We then apply the SSM to model real genomic data (both DNase-seq and MNase-seq data). Through incrementally increasing the types of genomic data in the SSM, we show that different data types can contribute complementary information for the inference of protein binding landscape and that the most accurate inference comes from modeling all available datasets.

This dissertation provides a foundation for future research by taking a step toward the genome-wide inference of protein-DNA interaction landscape through data integration.

Genome- and transcriptome-assisted development of nuclear insertion/deletion markers for Calanus species (Copepoda: Calanoida) identification

Copepods of the genus Calanus are key zooplankton species in temperate to arctic marine ecosystems. Despite their ecological importance, species identification remains challenging. Furthermore, the recent report of hybrids among Calanus species highlights the need for diagnostic nuclear markers to efficiently identify parental species and hybrids. Using next-generation sequencing analysis of both the genome and transcriptome from two sibling species, Calanus finmarchicus and Calanus glacialis, we developed a panel of 12 nuclear insertion/deletion markers. All the markers showed species-specific amplicon length. Furthermore, most of the markers were successfully amplified in other Calanus species, allowing the molecular identification of Calanus helgolandicus, Calanus hyperboreus and Calanus marshallae.

Next Generation Sequencing Reveals the Hidden Diversity of Zooplankton Assemblages

Zooplankton play an important role in our oceans, in biogeochemical cycling and providing a food source for commercially important fish larvae. However, difficulties in correctly identifying zooplankton hinder our understanding of their roles in marine ecosystem functioning, and can prevent detection of long term changes in their community structure. The advent of massively parallel next generation sequencing technology allows DNA sequence data to be recovered directly from whole community samples. Here we assess the ability of such sequencing to quantify richness and diversity of a mixed zooplankton assemblage from a productive time series site in the Western English Channel. Methodology/Principle Findings Plankton net hauls (200 µm) were taken at the Western Channel Observatory station L4 in September 2010 and January 2011. These samples were analysed by microscopy and metagenetic analysis of the 18S nuclear small subunit ribosomal RNA gene using the 454 pyrosequencing platform. Following quality control a total of 419,041 sequences were obtained for all samples. The sequences clustered into 205 operational taxonomic units using a 97% similarity cut-off. Allocation of taxonomy by comparison with the National Centre for Biotechnology Information database identified 135 OTUs to species level, 11 to genus level and 1 to order, <2.5% of sequences were classified as unknowns. By comparison a skilled microscopic analyst was able to routinely enumerate only 58 taxonomic groups. Conclusions Metagenetics reveals a previously hidden taxonomic richness, especially for Copepoda and hard-to-identify meroplankton such as Bivalvia, Gastropoda and Polychaeta. It also reveals rare species and parasites. We conclude that Next Generation Sequencing of 18S amplicons is a powerful tool for elucidating the true diversity and species richness of zooplankton communities. While this approach allows for broad diversity assessments of plankton it may become increasingly attractive in future if sequence reference libraries of accurately identified individuals are better populated.

Components of resistance to Fusarium head blight in wheat detected in a seed germination assay with Microdochium majus and the relationship to FHB disease development and mycotoxin accumulation from Fusarium graminearum infection.

Components of partial disease resistance (PDR) to fusarium head blight (FHB), detected in a seed-germination assay, were compared with whole-plant FHB resistance of 30 USA soft red winter wheat entries in the 2002 Uniform Southern FHB Nursery. Highly significant (P <0·001) differences between cultivars in the in vitro seed-germination assay inoculated with Microdochium majus were correlated to FHB disease incidence (r = -0·41; P <0·05), severity (r = -0·47; P <0·01), FHB index (r = -0·46; P <0·01), damaged kernels (r = -0·52; P <0·01), grain deoxynivalenol (DON) concentration (r = -0·40; P <0·05) and incidence/severity/kernel-damage index (ISK) (r = -0·45; P <0·01) caused by Fusarium graminearum. Multiple linear regression analysis explained a greater percentage of variation in FHB resistance using the seed-germination assay and the previously reported detached-leaf assay PDR components as explanatory factors. Shorter incubation periods, longer latent periods, shorter lesion lengths in the detached-leaf assay and higher germination rates in the seed-germination assay were related to greater FHB resistance across all disease variables, collectively explaining 62% of variation for incidence, 49% for severity, 56% for F. graminearum-damaged kernels (FDK), 39% for DON and 59% for ISK index. Incubation period was most strongly related to disease incidence and the early stages of infection, while resistance detected in the seed germination assay and latent period were more strongly related to FHB disease severity. Resistance detected using the seed-germination assay was notable as it related to greater decline in the level of FDK and a smaller reduction in DON than would have been expected from the reduction in FHB disease assessed by visual symptoms.