中国青藏高原青稞抗穗发芽资源评价与α-淀粉酶基因的克隆及表达

穗发芽（PHS，preharvest sprouting）是影响禾本科作物生产的重要的灾害之一。收获时期如遇潮湿天气容易导致穗发芽发生。发生穗发芽的种子内部水解酶（主要是α-淀粉酶）活性急剧升高，胚乳贮藏物质开始降解，造成作物产量和品质严重降低。因此，选育低穗发芽风险的品种是当前作物育种工作中面临的重要任务。 青稞(Hordeum vulgare ssp. vulgare)主要分布于青藏高原，自古以来就是青藏高原人民的主要粮食。近年来，由于青稞丰富的营养成分和特有的保健品质、在燃料工业中的潜力以及在啤酒酿造工业中的利用前景，在发达国家日趋受到重视，掀起综合研究利用的热潮。我国拥有占全世界2/3 以上的青稞资源，具有发展青稞产业的得天独厚的条件。然而，由于青稞收获期间恰逢青藏高原雨季来临，常有穗发芽灾害发生，使青稞生产损失巨大。目前对青稞穗发芽研究很少，适用于育种的穗发芽抗性材料相对缺乏，不能很好的满足青稞穗发芽抗性育种的需要。本研究以青藏高原青稞为材料，对其穗发芽抗性的评价指标和体系进行构建，同时筛选青稞抗穗发芽品种并对其抗性进行评价，还利用分子生物学手段对青稞穗发芽抗性的分子机理进行了初步探讨。主要研究结果如下： 1. 本试验以来自于我国青藏高原地区的青稞为材料，对休眠性测定的温度范围进行探讨，并对各种穗发芽抗性测定方法的对青稞的适用性进行评测。通过探讨温度对13 个不同基因型的青稞籽粒发芽和休眠性表达的影响，对筛选青稞抗穗发芽资源的温度条件进行探索，并初步分析了其休眠性表达的机理。在10，15，20，25，30℃的黑暗条件下，选用新收获的13 个青稞品种为材料进行籽粒发芽实验，以发芽指数（GI）评价其休眠性。结果发现，不同品种对温度敏感性不同，其中温度不敏感品种，在各温度条件下均表现很低的休眠性；而温度敏感品种，其休眠性表达受低温抑制，受高温诱导。15℃至25℃是进行青稞休眠性鉴定的较适宜的温度范围。通过对供试材料发芽后的α-淀粉酶活性，发现温度对青稞种子的休眠性表达的影响至少在一定程度上表现在对α-淀粉酶活性的调控上。随后，对分别在马尔康和成都进行种植的34 份青稞穗发芽指数（SI），穗发芽率（SR），籽粒发芽指数（GI）和α-淀粉酶活性（AA）进行了测定和分析，发现它们均受基因型×栽培地点的极显著影响，且四个参数之间具有一定相关性。GI 参数由于其变异系数较低，在不同栽培地点稳定性好，且操作简便，是较可靠和理想的穗发芽评价参数。SI 参数可作为辅助，区别籽粒休眠性相似的材料（基因型）或全面评价材料（基因型）的穗发芽抗性特征。AA 参数稳定性较差，并且检测方法复杂，因此不建议在育种及大量材料筛选和评价时使用。此外，青稞穗发芽抗性受环境影响较大，评价时应考虑到尽可能多的抗性影响因素及其在不同栽培条件下的变异。 2. 对来自青藏高原的青稞穗发芽抗性特征及其与其它农艺性状间的关系进行研究。通过测定穗发芽指数（SI）、籽粒发芽指数（GI）和α-淀粉酶活性（AA），表明113 份青稞材料的穗发芽抗性具有显著差异。SI、GI 和AA 参数的变幅分别为1.00~8.86、0.01~0.97 和0.00~2.76，其均值分别为4.72、0.63 和1.22。根据SI 参数，六个基因型，包括‘XQ9-5’，‘XQ33-9’，‘XQ37-5’，‘XQ42-9’，‘XQ45-7’和‘JCL’被鉴定为抗性品种。综合SI、GI 和AA 参数，可以发现青稞的穗发芽抗性机制包含颖壳等穗部结构的抗性和种子自身的抗性（即种子休眠性），且供试材料中未发现较强的胚休眠品种，除‘XQ45-7’外，所有品种在发芽第四天均能检测出α-淀粉酶活性。穗部结构和种子休眠的抗性机制因基因型不同而不同，在穗发芽抗性中可单独作用或共同作用。农家品种和西藏群体分别比栽培品种和四川群体的穗发芽抗性强，而在不同籽粒颜色的青稞中未发现明显差异。相关性检验发现，青稞的穗发芽抗性，主要是种子休眠性，与百粒重、开花期、成熟期、穗长、芒长和剑叶长呈显著负相关关系，与株高相关性不显著。农艺性状可以作为穗发芽抗性材料选育中的辅助指标。本试验为青稞穗发芽抗性育种研究提供了必要的理论基础和可供使用的亲本材料。 3. α-淀粉酶是由多基因家族编码的蛋白质，在植物种子萌发时高度表达，与植物种子的萌发能力密切相关。在大麦种子发芽时，高等电点α-淀粉酶的活性远大于低等电点的α-淀粉酶。为了研究不同穗发芽抗性青稞品种中编码高等电点α-淀粉酶Amy1 基因结构与抗性间的关系，我们以筛选得到的抗性品种‘XQ32-5’（TR1）、‘XQ37-5’（TR2）、‘XQ45-7’（TR3），易感品种‘97-15’（TS1）、‘9657’（TS2）以及强休眠大麦品种‘SAMSON’（SAM）为材料，对其Amy1 基因的编码区序列进行克隆和结构分析，并对它们推导的氨基酸序列进行比较。结果显示，青稞Amy1 基因具有三个外显子、两个内含子，编码区中有13 个核苷酸变异位点，均位于2、3 号外显子，2 个变异位点位于2 号外显子。SAM 和TS1 分别在2 号外显子相应位置有5 个相同的碱基(GAACT)的插入片段。相应α-淀粉酶氨基酸序列推导发现，所有核苷酸变异中有8 个导致相应氨基酸残基的改变，其余位点为同义突变。青稞Amy1 基因编码区序列品种间相似度高达99%以上，部分序列变异可能与其穗发芽抗性有关。随后，我们又通过SYBR Green 荧光定量技术对该基因在不同发芽时间（1d~7d）的相对表达水平进行了差异性检测。结果发现，7 天内不能检测到SAM 的Amy1 基因表达，5 个青稞品种间的Amy1 基因的相对表达量均随着发芽时间延长而上升，但上升方式有所不同。弱抗品种该基因表达更早，转录本增加速率更大，且在4~5 天可达到平台期。发芽7 天中，抗性品种总转录水平明显低于易感品种。本研究结果表明，青稞Amy1 基因的转录水平是与其穗发芽抗性高度相关。 我国青藏高原青稞，尤其是农家品种的穗发芽抗性具有丰富的变异，蕴藏着穗发芽抗性育种的宝贵资源。本研究为青稞穗发芽抗性育种建立了合理抗性评价体系，筛选出可供育种使用的特殊材料，阐明了农艺性状可辅助穗发芽抗性育种，同时还对穗发芽抗性与α-淀粉酶基因的结构和表达关系进行分析，为青稞穗发芽抗性资源筛选奠定了基础。 Preharvest sprouting (PHS) is a serious problem in crop production. It often takes place when encountering damp, cold conditions at harvest time and results in the decrease of grain quality and great loss of yield by triggering the synthesis of endosperm degrading enzymes (mostly the α-amylase). Therefore, PHS is regarded as an important criterion for crop breeding. In order to minimize the risk of PHS, resistant genotypes are highly required. Hulless barley (Hordeum vulgare ssp. vulgare) is the staple food crop in Qinghai-Tibetan Plateau from of old, where is one of the origin and genetic diversity centers of hulless barley. Recently, interest in hulless barley has been sparked throughout the world due to the demonstrations of its great potential in health food industry and fuel alcohol production. Indeed, hulless barley can also be utilized to produce good quality malt if the appropriate malting conditions are used. In China, overcast and rainy conditions often occur at maturity of hulless barley and cause an adverse on its production and application. PHS resistant genotypes, therefore, are highly required for the hulless barley breeding programs. However, few investigations have been made so far on this issue. The objectives of this study were: 1) to assessment of methods used in testing preharvest sprouting resistance in hulless barley; 2) to evaluate the variability and characteristics of PHS resistance of hulless barley from Qinghai-Tibet Plateau in China; 3) to select potential parents for PHS resistance breeding; 4) to primarily study on the molecular mechanism of PHS resistance of hulless barley. Our results are as followed: 1. We investigated the temperature effects on seed germination and seed dormancy expression of hulless barley, discussed appropriate temperature range for screening of PHS resistant varieties, and analyzed the mechanism of seed dormancy expression of hulless barley. The dormancy level of 13 hulless barley were evaluated by GI (germination index) values calculating by seed germination tests at temperature of 10，15，20，25，30℃ in darkness. There were great differences in temperature sensitivity among these accessions. The insensitive accessions showed low dormancy at any temperature while the dormancy expression of sensitive accessions could be restrained by low temperature and induced by high temperature. The temperature range of 15℃ to 25℃ was workable for estimating of dormancy level of hulless barley according to our data. Analysis of α-amylase activity showed that the temperature effects on seed germination and the expression of seed dormancy be achieved probable via regulating of α-amylase activity. Furthermore, we evaluated the differences in sprouting index (SI), sprouting rate (SR), germination index (GI) and α-amylase activity (AA) between Maerkang and Chengdu among 34 accessions of hulless barley from Qinghai-Tibetan Plateau in China. These PHS sprouting parameters were significantly affected by accession×location, and they had correlation between each other. GI was the most reliable parameter because of its low CV value, good repeatability and simple operation. SI could assist in differentiating between accessions of similar dormancy or overall evaluation of the resistance. AA was bad in repeatability and had relatively complex testing method, therefore, not appropriate for breeding and evaluation and screening of PHS resistant materials. Besides, since PHS resistance of hulless barley was greatly influenced by its growth environment, possibly much influencing factors and variations between cultivated conditions should be considered. 2. In this study, large variation was found among 113 genotypes of hulless barley (Hordeum vulgare ssp.vulgare) from Qinghai-Tibetan Plateau in China, based on the sprouting index (SI), germination index (GI) and α-amylase activity (AA) which derived from sprouting test of intact spikes, germination test of threshed seeds and determination of α-amylase activity, respectively. The range of SI, GI and AA was 1.00~8.86, 0.01~0.97 and 0.00~2.76，the mean was 4.72, 0.63 and 1.22 espectively. Six resistant genotypes, including ‘XQ9-5’, ‘XQ33-9’, ‘XQ37-5’, ‘XQ42-9’, ‘XQ45-7’ and ‘JCL’, were identified based on SI. Integrating the three parameters, it was clear that both hulls and seeds involved in PHS resistance in intact spikes of hulless barley and there was no long-existent embryo dormancy found among the test genotypes. All the genotypes, except ‘XQ45-7’, had detectable α-amylase activity on the 4th day after germination. There was PHS resistance imposed by the hull and seed per se and the two factors can act together or independent of each other. Besides, landraces or Tibet hulless barley had a wider variation and relatively more PHS resistance when compared with cultivars or Sichuan hulless barley. No significant difference was found among hulless barley of different seed colors. The correlation analysis showed PHS resistance was negatively related to hundred grain weight, days to flowering, days to maturity, spike length, awn length and flag length but not related to plant height. This study provides essential information and several donor parents for breeding of resistance to PHS. 3. Alpha-amylase isozymes are encoded by a family of multigenes. They highly express in germinating seeds and is closely related to seed germination ability. In barley germinating seeds, the activity of high pI α-amylase is much higher than low pI α-amylase. The aim of this study was to determine the relationship between preharvest sprouting resistance of hulless barley and the gene structure of Amy1 gene which encodes high pI α-amylase. The coding region and cDNA of Amy1 gene of three resistant accessions, including ‘XQ32-5’ (TR1), ‘XQ37-5’ (TR2), ‘XQ45-7’ (TR3), two susceptible accessions ‘97-15’ (TS1), ‘9657’ (TS2) and one highly dormant barley accession ‘SAMSON’ (SAM) was cloned. Analysis of their DNA sequences revealed there were three exons and two introns in Amy1 gene. Thirteen variable sites were in exon2 and exon3, 2 variable sites were in intron2. SAM and TS1 had a GAACT insert segment in the same site in intron2. Only 8 variable sites caused the change of amino acid residues. There were 99% of similarity between the tested hulless barley and some of the variable sites might be related with preharvest sprouting resistance. Then, we investigated the expression level of Amy1 gene in the 7-day germination test. Results of quantitative real-time PCR indicated that the relative expression trends of Amy1 gene were the same but had significant differences in the increase fashion between hulless barleys and no detectable expression was found in SAM. Susceptible accessions had earlier expression and faster increase and reached the maximum on day 4 ~ day 5. Besides, total transcripts level was found lower in resistant accessions than susceptible accessions. This study indicated that α-amylase activity was highly related to the transcription level of Amy1 gene which not correlated to missense mutation sites. In conclusion, hulless barley, especially the landraces from Qinghai-Tibetan Plateau in China possesses high degree of variation in PHS performance, which indicates the potential of Tibetan hulless barley as a good source for breeding of resistance to PHS. This study provides several donor parents for breeding of resistance to PHS. Our results also demonstrate that agronomic traits may be used as assistants for PHS resistance selection in hulless barley. Besides, analysis of high pI α-amylase coding gene Amy1 revealed the relative high expression of was Amy1 one of the mainly reason of different PHS resistance level in hulless barley.

Differences in the rDNA-bearing chromosome divide the Asian-Pacific and Atlantic species of Crassostrea (Bivalvia, Mollusca)

Karyotype and chromosomal location of the major ribosomal RNA genes (rDNA) were studied using fluorescence in situ hybridization (FISH) in five species of Crassostrea: three Asian-Pacific species (C. gigas, C. plicatula, and C. ariakensis) and two Atlantic species (C. virginica and C. rhizophorae). FISH probes were made by PCR amplification of the intergenic transcribed spacer between the 18S and 5.8S rRNA genes, and labeled with digoxigenin-11-dUTP. All five species had a haploid number of 10 chromosomes. The Atlantic species had 1-2 submetacentric chromosomes, while the three Pacific species had none. FISH with metaphase chromosomes detected a single telomeric locus for rDNA in all five species without any variation. In all three Pacific species, rDNA was located on the long arm of Chromosome 10 (10q)-the smallest chromosome. In the two Atlantic species, rDNA was located on the short arm of Chromosome 2 (2p)-the second longest chromosome. A review of other studies reveals the same distribution of NOR sites (putative rDNA loci) in three other species: on 10q in C. sikamea and C. angulata from the Pacific Ocean and on 2p in C. gasar from the western Atlantic. All data support the conclusion that differences in size and shape of the rDNA-bearing chromosome represent a major divide between Asian-Pacific and Atlantic species of Crassostrea. This finding suggests that chromosomal divergence can occur under seemingly conserved karyotypes and may play a role in reproductive isolation and speciation.

Construction and characterization of two bacterial artificial chromosome libraries of Zhikong scallop, Chlamys farreri Jones et Preston, and identification of BAC clones containing the genes involved in its innate immune system

Large-insert bacterial artificial chromosome (BAC) libraries are necessary for advanced genetics and genomics research. To facilitate gene cloning and characterization, genome analysis, and physical mapping of scallop, two BAC libraries were constructed from nuclear DNA of Zhikong scallop, Chlamys farreri Jones et Preston. The libraries were constructed in the BamHI and MboI sites of the vector pECBAC1, respectively. The BamHI library consists of 73,728 clones, and approximately 99% of the clones contain scallop nuclear DNA inserts with an average size of 110 kb, covering 8.0x haploid genome equivalents. Similarly, the MboI library consists of 7680 clones, with an average insert of 145 kb and no insert-empty clones, thus providing a genome coverage of 1.1x. The combined libraries collectively contain a total of 81,408 BAC clones arrayed in 212 384-well microtiter plates, representing 9.1x haploid genome equivalents and having a probability of greater than 99% of discovering at least one positive clone with a single-copy sequence. High-density clone filters prepared from a subset of the two libraries were screened with nine pairs of Overgos designed from the cDNA or DNA sequences of six genes involved in the innate immune system of mollusks. Positive clones were identified for every gene, with an average of 5.3 BAC clones per gene probe. These results suggest that the two scallop BAC libraries provide useful tools for gene cloning, genome physical mapping, and large-scale sequencing in the species.

Classification of jinjiang oysters Crassostrea rivularis (Gould, 1861) from China, based on morphology and phylogenetic analysis

The jinjiang oyster Crassostrea rivularis [Gould, 1861. Descriptions of Shells collected in the North Pacific Exploring Expedition under Captains Ringgold and Rodgers. Proc. Boston Soc. Nat. Hist. 8 (April) 33-40] is one of the most important and best-known oysters in China. Based on the color of its flesh, two forms of C rivularis are recognized and referred to as the "white meat" and 11 red meat" oysters. The classification of white and red forms of this species has been a subject of confusion and debate in China. To clarify the taxonomic status of the two forms of C. rivularis, we collected and analyzed oysters from five locations along China's coast using both morphological characters and DNA sequences from mitochondrial 16S rRNA and cytochrome oxidase 1, and the nuclear 28S rRNA genes. Oysters were classified as white or red forms according to their morphological characteristics and then subjected to DNA sequencing. Both morphological and DNA sequence data suggest that the red and white oysters are two separate species. Phylogenetic analysis of DNA sequences obtained in this study and existing sequences of reference species show that the red oyster is the same species as C. ariakensis Wakiya [1929. Japanese food oysters. Jpn. J. Zool. 2, 359-367.], albeit the red oysters from north and south China are genetically distinctive. The white oyster is the same species as a newly described species from Hong Kong, C. hongkongensis Lam and Morton [2003. Mitochondrial DNA and identification of a new species of Crassostrea (Bivalvia: Ostreidae) cultured for centuries in the Pearl River Delta, Hong Kong, China. Aqua. 228, 1-13]. Although the name C. rivularis has seniority over C. ariakensis and C. hongkongensis, the original description of Ostrea rivularis by Gould [1861] does not fit shell characteristics of either the red or the white oysters. We propose that the name of C. rivularis Gould [1861] should be suspended, the red oyster should take the name C. ariakensis, and the white oyster should take the name C. hongkongensis. (C) 2004 Elsevier B.V. All rights reserved.

Classification of common oysters from North China

Oysters are commonly found on rocky shores along China's northern coast, although there is considerable confusion as to what species they are. To determine the taxonomic status of these oysters, we collected specimens from nine locations north of the Yangtze River and conducted genetic identification using DNA sequences. Fragments from three genes, mitochondrial 165 rRNA, mitochondria! cytochrome oxidase I (COI), and nuclear 285 rRNA, were sequenced in six oysters from each of the nine sites. Phylogenetic analysis of all three gene fragments clearly demonstrated that the small oysters commonly found on intertidal rocks in north China are Crassostrea gigas (Thunberg, 1793), not C. plicatula (the zhe oyster) as widely assumed. Their small size and irregular shell characteristics are reflections of the stressful intertidal environment they live in and not reliable characters for classification. Our study confirms that the oysters from Weifang, referred to as Jinjiang oysters or C. rivularis (Gould, 1861), are C. ariakensis (Wakiya, 1929). We found no evidence for the existence of C. talienwhanensis (Crosse, 1862) and other Crassostrea species in north China. Our study highlights the need for reclassifying oysters of China with molecular data.

Identification of Crassostrea ariakensis and related oysters by multiplex species-specific PCR

Genetic markers are needed for rapid and reliable identification of oysters. In this study, we developed multiplex genus- and species-specific PCR markers for the identification of oysters from China. We used the mitochondrial cytochrome oxidase I (COI) and nuclear 28S ribosomal RNA genes for marker development. DNA sequences from different species were obtained from GenBank or by direct sequencing. Sequences were aligned, and genus- and species-specific nucleotides were identified. Primers were designed for genus/species-specific amplification to generate fragments of different sizes. A multiplex set of genus- and species-specific primers from the 28S gene was able to separate C. ariakensis and C. hongkongensis from other species and assign oysters to four genera. A set of species-specific COI primers provided positive identification of all five Crassostrea species from China, C. ariakensis, C. hongkongensis, C. angulata, C. gigas, and C. sikamea in a single PCR. The multiplex PCR assays do not require fluorescence-labeling or post-PCR enzyme digestion, providing a simple, fast and reliable method for the identification of oysters from China.

Molecular authentication of the traditional Tibetan medicinal plant Swertia mussotii

Swertia mussotii is an important species in Tibetan folk medicine. However, it is quite expensive and frequently adulterated, so reliable methods for authentication of putative specimens and preparations of the species are needed to protect consumers and to support conservation measures. We show here that the chloroplast (cp) DNA rpl16 intron has limited utility for differentiating S. mussotii from closely related species, since the cpDNA rpl16 sequences are identical in S. mussotii and two other species of Swertia. However, the rDNA internal transcribed spacer (ITS) sequences differ significantly between S. mussotii and all of 13 tested potential adulterants. Thus, the ITS region provides a robust molecular marker for differentiating the medicinal S. mussotii from related adulterants. Therefore, a pair of allele-specific diagnostic primers based on the divergent ITS region was designed to distinguish S. mussotii from the other species. Authentication by allele-specific diagnostic PCR using these primers is convenient, effective and both simpler and less time-consuming than sequencing the ITS region.

Molecular phylogeny and biogeography of the Qinghai-Tibet Plateau endemic Nannoglottis (Asteraceae)

All taxa endemic to the Qinghai-Tibet Plateau are hypothesized to have originated in situ or from immediately adjacent areas because of the relatively recent formation of the plateau since the Pliocene, followed by the large-scaled biota extinction and recession caused by the Quaternary ice sheet. However, identification of specific progenitors remains difficult for some endemics, especially some endemic genera. Nannoglottis, with about eight species endemic to this region, is one such genus. Past taxonomic treatments have suggested its relationships with four different tribes of Asteraceae. We intend to identify the closest relatives of Nannoglottis by evaluating the level of monophyly, tribal delimitation, and systematic position of the genus by using molecular data from ndhF gene, trnL-F, and ITS region sequences. We find that all sampled species of Nannoglottis are a well-defined monophyly. This supports all recent taxonomic treatments of Nannoglottis, in which all sampled species were placed in one broadly re-circumscribed genus. Nannoglottis is most closely related to the Astereae, but stands as an isolated genus as the first diverging lineage of the tribe, without close relatives. A tentative relationship was suggested for Nannoglottis and the next lineage of the tribe was based on the ITS topology, the "basal group," which consists of seven genera from the Southern Hemisphere. Such a relationship is supported by some commonly shared plesiomorphic morphological characters. Despite the very early divergence of Nannoglottis in the Astereae, the tribe must be regarded to have its origin in Southern Hemisphere rather than in Asia, because based on all morphological, molecular, biogeographical, and fossil data, the Asteraceae and its major lineages (tribes) are supposed to have originated in the former area. Long-distance dispersal using Southeast Asia as a steppingstone from Southern Hemisphere to the Qinghai-Tibet Plateau is the most likely explanation for this unusual biogeographic link of Nannoglottis. The 23-32-million-year divergence time between Nannoglottis and the other Astereae estimated by DNA sequences predated the formation of the plateau. This estimation is further favored by the fossil record of the Asteraceae and the possible time of origin of the Astereae. Nannoglottis seems to have reached the Qinghai-Tibet area in the Oligocene-Eocene and then re-diversified with the uplift of the plateau. The molecular infragenetic phylogeny of the genus identifies two distinct clades, which reject the earlier infrageneric classification based on the arrangement of the involucral bracts and the length of the ligules, but agree well with the habits and ecological preferences of its current species. The "alpine shrub" vs. "coniferous forest" divergence within Nannoglottis was estimated at about 3.4 million years ago when the plateau began its first large-scale uplifting and the coniferous vegetation began to appear. Most of the current species at the "coniferous forest" clade of the genus are estimated to have originated from 1.02 to 1.94 million years ago, when the second and third uprisings of the plateau occurred, the climate oscillated and the habitats were strongly changed. The assumed evolution, speciation diversity, and radiation of Nannoglottis based on molecular phylogeny and divergence times agree well with the known geological and paleobotanical histories of the Qinghai-Tibet Plateau. (C) 2002 Elsevier Science (USA). All rights reserved.

Origin of avian genome size and structure in non-avian dinosaurs

Chris L. Organ, Andrew M. Shedlock, Andrew Meade, Mark Pagel and Scott V. Edwards (2007). Origin of avian genome size and structure in non-avian dinosaurs. Nature, 46(7132), 180-184. RAE2008

Genetic analysis of bacteriophages from clinical and environmental samples

Bacteriophages, viruses infecting bacteria, are uniformly present in any location where there are high numbers of bacteria, both in the external environment and the human body. Knowledge of their diversity is limited by the difficulty to culture the host species and by the lack of the universal marker gene present in all viruses. Metagenomics is a powerful tool that can be used to analyse viral communities in their natural environments. The aim of this study was to investigate diverse populations of uncultured viruses from clinical (a sputum of patient with cystic fibrosis, CF) and environmental samples (a sludge from a dairy food wastewater treatment plant) containing rich bacterial populations using genetic and metagenomic analyses. Metagenomic sequencing of viruses obtained from these samples revealed that the majority of the metagenomic reads (97-99%) were novel when compared to the NCBI protein database using BLAST. A large proportion of assembled contigs were assignable as novel phages or uncharacterised prophages, the next largest assignable group being single-stranded eukaryotic virus genomes. Sputum from a cystic fibrosis patient contained DNA typical of phages of bacteria that are traditionally involved in CF lung infections and other bacteria that are part of the normal oral flora. The only eukaryotic virus detected in the CF sputum was Torque Teno virus (TTV). A substantial number of assigned sequences from dairy wastewater could be affiliated with phages of bacteria that are typically found in the soil and aquatic environments, including wastewater. Eukaryotic viral sequences were dominated by plant pathogens from the Geminiviridae and Nanoviridae families, and animal pathogens from the Circoviridae family. Antibiotic resistance genes were detected in both metagenomes suggesting phages could be a source for transmissible antimicrobial resistance. Overall, diversity of viruses in the CF sputum was low, with 89 distinct viral genotypes predicted, and higher (409 genotypes) in the wastewater. Function-based screening of a metagenomic library constructed from DNA extracted from dairy food wastewater viruses revealed candidate promoter sequences that have ability to drive expression of GFP in a promoter-trap vector in Escherichia coli. The majority of the cloned DNA sequences selected by the assay were related to ssDNA circular eukaryotic viruses and phages which formed a minority of the metagenome assembly, and many lacked any significant homology to known database sequences. Natural diversity of bacteriophages in wastewater samples was also examined by PCR amplification of the major capsid protein sequences, conserved within T4-type bacteriophages from Myoviridae family. Phylogenetic analysis of capsid sequences revealed that dairy wastewater contained mainly diverse and uncharacterized phages, while some showed a high level of similarity with phages from geographically distant environments.

Bayesian variable selection in structured high-dimensional covariate spaces with applications in genomics

We consider the problem of variable selection in regression modeling in high-dimensional spaces where there is known structure among the covariates. This is an unconventional variable selection problem for two reasons: (1) The dimension of the covariate space is comparable, and often much larger, than the number of subjects in the study, and (2) the covariate space is highly structured, and in some cases it is desirable to incorporate this structural information in to the model building process. We approach this problem through the Bayesian variable selection framework, where we assume that the covariates lie on an undirected graph and formulate an Ising prior on the model space for incorporating structural information. Certain computational and statistical problems arise that are unique to such high-dimensional, structured settings, the most interesting being the phenomenon of phase transitions. We propose theoretical and computational schemes to mitigate these problems. We illustrate our methods on two different graph structures: the linear chain and the regular graph of degree k. Finally, we use our methods to study a specific application in genomics: the modeling of transcription factor binding sites in DNA sequences. © 2010 American Statistical Association.

Are the native giant tortoises from the Seychelles really extinct? A genetic perspective based on mtDNA and microsatellite data.

The extinction of the giant tortoises of the Seychelles Archipelago has long been suspected but is not beyond doubt. A recent morphological study of the giant tortoises of the western Indian Ocean concluded that specimens of two native Seychelles species survive in captivity today alongside giant tortoises of Aldabra, which are numerous in zoos as well as in the wild. This claim has been controversial because some of the morphological characters used to identify these species, several measures of carapace morphology, are reputed to be quite sensitive to captive conditions. Nonetheless, the potential survival of giant tortoise species previously thought extinct presents an exciting scenario for conservation. We used mitochondrial DNA sequences and nuclear microsatellites to examine the validity of the rediscovered species of Seychelles giant tortoises. Our results indicate that the morphotypes suspected to represent Seychelles species do not show levels of variation and genetic structuring consistent with long periods of reproductive isolation. We found no variation in the mitochondrial control region among 55 individuals examined and no genetic structuring in eight microsatellite loci, pointing to the survival of just a single lineage of Indian Ocean tortoises.

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation.

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Genomic and Functional Variation at a Normal Human Centromere

Centromeres are chromosomal loci essential for genome stability. Their malfunction can cause chromosome instability associated with cancer, infertility, and birth defects. This study focused on an intriguing centromere on human chromosome 17, which displays normal functional variation. Centromere identity can be found on either of two large arrays of repetitive DNA. We investigated inter-individual sequence variation on these two arrays and found association between array size, array variation, and centromere function. Our data suggest a functional influence of DNA sequence at this critical epigenetic locus.