Haplotype analysis of vernalization loci in European barley germplasm reveals novel VRN-H1 alleles and a predominant winter VRN-H1/VRN-H2 multi-locus haplotype

In barley, variation in the requirement for vernalization (an extended period of low temperature before flowering can occur) is determined by the VRN-H1, -H2 and -H3 loci. In European cultivated germplasm, most variation in vernalization requirement is accounted for by alleles at VRN-H1 and VRN-H2 only, but the range of allelic variation is largely unexplored. Here we characterise VRN-H1 and VRN-H2 haplotypes in 429 varieties representing a large portion of the acreage sown to barley in Western Europe over the last 60 years. Analysis of genotype, intron I sequencing data and growth habit tests identified three novel VRN-H1 alleles and determined the most frequent VRN-H1 intron I rearrangements. Combined analysis of VRN-H1 and VRN-H2 alleles resulted in the classification of seventeen VRN-H1/VRN-H2 multi-locus haplotypes, three of which account for 79% of varieties. The molecular markers employed here represent powerful diagnostic tools for prediction of growth habit and assessment of varietal purity. These markers will also allow development of germplasm to test the behaviour of individual alleles with the aim of understanding the relationship between allelic variation and adaptation to specific agri-environments.

Epidemiology and Genetic Variability of Human Metapneumovirus During a 4-Year-Long Study in Southeastern Brazil

Epidemiological and molecular characteristics of human metapneumovirus (hMPV) were compared with human respiratory syncytial virus (hRSV) in infants and young children admitted for acute lower respiratory tract infections in a prospective study during four consecutive years in subtropical Brazil. GeneScan polymerase chain assays (GeneScan RT-PCR) were used to detect hMPV and hRSV in nasopharyngeal aspirates of 1,670 children during January 2003 to December 2006. hMPV and hRSV were detected, respectively, in 191 (11.4%) and in 702 (42%) of the children admitted with acute lower respiratory tract infections at the Sao Paulo University Hospital. Sequencing data of the hMPV F gene revealed that two groups of the virus, each divided into two subgroups, co-circulated during three consecutive years. It was also shown that a clear dominance of genotype B1 occurred during the years 2004 and 2005, followed by genotype A2 during 2006. J. Med. Virol. 81:915-921,2009. (C) 2009 Wiley-Liss, Inc.

Yeasts isolated from a fungus-growing ant nest, including the description of Trichosporon chiarellii sp nov., an anamorphic basidiomycetous yeast

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Chromosomal similarity between the Scaly-headed parrot (Pionus maximiliani), the Short-tailed parrot (Graydidascalus brachyurus) and the Yellow-faced parrot (Salvatoria xanthops) (Psittaciformes: Aves): a cytotaxonomic analysis

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Estudo moleculares das regiões cromossômicas 18p e 18q proximal em portadores de gagueira persistente da familial

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Avaliação do potencial biotecnológico de linhagens mutantes de Rhizobium tropici

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Narrow Fungal Mycorrhizal Diversity in a Population of the Orchid Coppensia doniana

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Estudo de transcriptomas por RNAseq em tecidos de cabeça e glândula salivar de Rhodnius montenegrensis e Rhodnius robustus (Hemiptera, Reduviidade, Triatominae)

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Low-Cost Ultra-Wide Genotyping Using Roche/454 Pyrosequencing for Surveillance of HIV Drug Resistance

Background: Great efforts have been made to increase accessibility of HIV antiretroviral therapy (ART) in low and middle-income countries. The threat of wide-scale emergence of drug resistance could severely hamper ART scale-up efforts. Population-based surveillance of transmitted HIV drug resistance ensures the use of appropriate first-line regimens to maximize efficacy of ART programs where drug options are limited. However, traditional HIV genotyping is extremely expensive, providing a cost barrier to wide-scale and frequent HIV drug resistance surveillance. Methods/Results: We have developed a low-cost laboratory-scale next-generation sequencing-based genotyping method to monitor drug resistance. We designed primers specifically to amplify protease and reverse transcriptase from Brazilian HIV subtypes and developed a multiplexing scheme using multiplex identifier tags to minimize cost while providing more robust data than traditional genotyping techniques. Using this approach, we characterized drug resistance from plasma in 81 HIV infected individuals collected in Sao Paulo, Brazil. We describe the complexities of analyzing next-generation sequencing data and present a simplified open-source workflow to analyze drug resistance data. From this data, we identified drug resistance mutations in 20% of treatment naive individuals in our cohort, which is similar to frequencies identified using traditional genotyping in Brazilian patient samples. Conclusion: The developed ultra-wide sequencing approach described here allows multiplexing of at least 48 patient samples per sequencing run, 4 times more than the current genotyping method. This method is also 4-fold more sensitive (5% minimal detection frequency vs. 20%) at a cost 3-5 x less than the traditional Sanger-based genotyping method. Lastly, by using a benchtop next-generation sequencer (Roche/454 GS Junior), this approach can be more easily implemented in low-resource settings. This data provides proof-of-concept that next-generation HIV drug resistance genotyping is a feasible and low-cost alternative to current genotyping methods and may be particularly beneficial for in-country surveillance of transmitted drug resistance.

Multicentric Carpotarsal Osteolysis Is Caused by Mutations Clustering in the Amino-Terminal Transcriptional Activation Domain of MAFB

Multicentric carpotarsal osteolysis (MCTO) is a rare skeletal dysplasia characterized by aggressive osteolysis, particularly affecting the carpal and tarsal bones, and is frequently associated with progressive renal failure. Using exome capture and next-generation sequencing in five unrelated simplex cases of MCTO, we identified previously unreported missense mutations clustering within a 51 base pair region of the single exon of MAFB, validated by Sanger sequencing. A further six unrelated simplex cases with MCTO were also heterozygous for previously unreported mutations within this same region, as were affected members of two families with autosomal-dominant MCTO. MAFB encodes a transcription factor that negatively regulates RANKL-induced osteoclastogenesis and is essential for normal renal development. Identification of this gene paves the way for development of novel therapeutic approaches for this crippling disease and provides insight into normal bone and kidney development.

A comparative transcriptome analysis reveals expression profiles conserved across three Eimeria spp. of domestic fowl and associated with multiple developmental stages

Coccidiosis of the domestic fowl is a worldwide disease caused by seven species of protozoan parasites of the genus Eimeria. The genome of the model species, Eimeria tenella, presents a complexity of 55-60 MB distributed in 14 chromosomes. Relatively few studies have been undertaken to unravel the complexity of the transcriptome of Eimeria parasites. We report here the generation of more than 45,000 open reading frame expressed sequence tag (ORESTES) cDNA reads of E. tenella, Eimeria maxima and Eimeria acervulina, covering several developmental stages: unsporulated oocysts, sporoblastic oocysts, sporulated oocysts, sporozoites and second generation merozoites. All reads were assembled to constitute gene indices and submitted to a comprehensive functional annotation pipeline. In the case of E. tenella, we also incorporated publicly available ESTs to generate an integrated body of information. Orthology analyses have identified genes conserved across different apicomplexan parasites, as well as genes restricted to the genus Eimeria. Digital expression profiles obtained from ORESTES/EST countings, submitted to clustering analyses, revealed a high conservation pattern across the three Eimeria spp. Distance trees showed that unsporulated and sporoblastic oocysts constitute a distinct clade in all species, with sporulated oocysts forming a more external branch. This latter stage also shows a close relationship with sporozoites, whereas first and second generation merozoites are more closely related to each other than to sporozoites. The profiles were unambiguously associated with the distinct developmental stages and strongly correlated with the order of the stages in the parasite life cycle. Finally, we present The Eimeria Transcript Database (http://www.coccidia.icb.usp.br/eimeriatdb), a website that provides open access to all sequencing data, annotation and comparative analysis. We expect this repository to represent a useful resource to the Eimeria scientific community, helping to define potential candidates for the development of new strategies to control coccidiosis of the domestic fowl. (C) 2011 Australian Society for Parasitology Inc. Published by Elsevier Ltd. All rights reserved.

Transposable elements: structure and dynamic of the autonomous retrotransposon R2 in Arthropoda with non-canonical reproduction

Eukaryotic ribosomal DNA constitutes a multi gene family organized in a cluster called nucleolar organizer region (NOR); this region is composed usually by hundreds to thousands of tandemly repeated units. Ribosomal genes, being repeated sequences, evolve following the typical pattern of concerted evolution. The autonomous retroelement R2 inserts in the ribosomal gene 28S, leading to defective 28S rDNA genes. R2 element, being a retrotransposon, performs its activity in the genome multiplying its copy number through a “copy and paste” mechanism called target primed reverse transcription. It consists in the retrotranscription of the element’s mRNA into DNA, then the DNA is integrated in the target site. Since the retrotranscription can be interrupted, but the integration will be carried out anyway, truncated copies of the element will also be present in the genome. The study of these truncated variants is a tool to examine the activity of the element. R2 phylogeny appears, in general, not consistent with that of its hosts, except some cases (e.g. Drosophila spp. and Reticulitermes spp.); moreover R2 is absent in some species (Fugu rubripes, human, mouse, etc.), while other species have more R2 lineages in their genome (the turtle Mauremys reevesii, the Japanese beetle Popilia japonica, etc). R2 elements here presented are isolated in 4 species of notostracan branchiopods and in two species of stick insects, whose reproductive strategies range from strict gonochorism to unisexuality. From sequencing data emerges that in Triops cancriformis (Spanish gonochoric population), in Lepidurus arcticus (two putatively unisexual populations from Iceland) and in Bacillus rossius (gonochoric population from Capalbio) the R2 elements are complete and encode functional proteins, reflecting the general features of this family of transposable elements. On the other hand, R2 from Italian and Austrian populations of T. cancriformis (respectively unisexual and hermaphroditic), Lepidurus lubbocki (two elements within the same Italian population, gonochoric but with unfunctional males) and Bacillus grandii grandii (gonochoric population from Ponte Manghisi) have sequences that encode incomplete or non-functional proteins in which it is possible to recognize only part of the characteristic domains. In Lepidurus couesii (Italian gonochoric populations) different elements were found as in L. lubbocki, and the sequencing is still in progress. Two hypothesis are given to explain the inconsistency of R2/host phylogeny: vertical inheritance of the element followed by extinction/diversification or horizontal transmission. My data support previous study that state the vertical transmission as the most likely explanation; nevertheless horizontal transfer events can’t be excluded. I also studied the element’s activity in Spanish populations of T. cancriformis, in L. lubbocki, in L. arcticus and in gonochoric and parthenogenetic populations of B. rossius. In gonochoric populations of T. cancriformis and B. rossius I found that each individual has its own private set of truncated variants. The situation is the opposite for the remaining hermaphroditic/parthenogenetic species and populations, all individuals sharing – in the so far analyzed samples - the majority of variants. This situation is very interesting, because it isn’t concordant with the Muller’s ratchet theory that hypothesizes the parthenogenetic populations being either devoided of transposable elements or TEs overloaded. My data suggest a possible epigenetic mechanism that can block the retrotransposon activity, and in this way deleterious mutations don’t accumulate.

STATISTICAL METHODS FOR RARE AND COMMON VARIANT ASSOCIATION STUDIES

Complex human diseases are a major challenge for biological research. The goal of my research is to develop effective methods for biostatistics in order to create more opportunities for the prevention and cure of human diseases. This dissertation proposes statistical technologies that have the ability of being adapted to sequencing data in family-based designs, and that account for joint effects as well as gene-gene and gene-environment interactions in the GWA studies. The framework includes statistical methods for rare and common variant association studies. Although next-generation DNA sequencing technologies have made rare variant association studies feasible, the development of powerful statistical methods for rare variant association studies is still underway. Chapter 2 demonstrates two adaptive weighting methods for rare variant association studies based on family data for quantitative traits. The results show that both proposed methods are robust to population stratification, robust to the direction and magnitude of the effects of causal variants, and more powerful than the methods using weights suggested by Madsen and Browning [2009]. In Chapter 3, I extended the previously proposed test for Testing the effect of an Optimally Weighted combination of variants (TOW) [Sha et al., 2012] for unrelated individuals to TOW &ndash F, TOW for Family &ndash based design. Simulation results show that TOW &ndash F can control for population stratification in wide range of population structures including spatially structured populations, is robust to the directions of effect of causal variants, and is relatively robust to percentage of neutral variants. In GWA studies, this dissertation consists of a two &ndash locus joint effect analysis and a two-stage approach accounting for gene &ndash gene and gene &ndash environment interaction. Chapter 4 proposes a novel two &ndash stage approach, which is promising to identify joint effects, especially for monotonic models. The proposed approach outperforms a single &ndash marker method and a regular two &ndash stage analysis based on the two &ndash locus genotypic test. In Chapter 5, I proposed a gene &ndash based two &ndash stage approach to identify gene &ndash gene and gene &ndash environment interactions in GWA studies which can include rare variants. The two &ndash stage approach is applied to the GAW 17 dataset to identify the interaction between KDR gene and smoking status.

"Hot spot" in the PROP1 gene responsible for combined pituitary hormone deficiency.

As pituitary function depends on the integrity of the hypothalamic-pituitary axis, any defect in the development and organogenesis of this gland may account for a form of combined pituitary hormone deficiency (CPHD). Although pit-1 was 1 of the first factors identified as a cause of CPHD in mice, many other homeodomain and transcription factors have been characterized as being involved in different developmental stages of pituitary gland development, such as prophet of pit-1 (prop-1), P-Lim, ETS-1, and Brn 4. The aims of the present study were first to screen families and patients suffering from different forms of CPHD for PROP1 gene alterations, and second to define possible hot spots and the frequency of the different gene alterations found. Of 73 subjects (36 families) analyzed, we found 35 patients, belonging to 18 unrelated families, with CPHD caused by a PROP1 gene defect. The PROP1 gene alterations included 3 missense mutations, 2 frameshift mutations, and 1 splice site mutation. The 2 reported frameshift mutations could be caused by any 2-bp GA or AG deletion at either the 148-GGA-GGG-153 or 295-CGA-GAG-AGT-303 position. As any combination of a GA or AG deletion yields the same sequencing data, the frameshift mutations were called 149delGA and 296delGA, respectively. All but 1 mutation were located in the PROP1 gene encoding the homeodomain. Importantly, 3 tandem repeats of the dinucleotides GA at location 296-302 in the PROP1 gene represent a hot spot for CPHD. In conclusion, the PROP1 gene seems to be a major candidate gene for CPHD; however, further studies are needed to evaluate other genetic defects involved in pituitary development.