An Analytical Workflow for Metagenomic Data and its Application to the Study of Chronic Obstructive Pulmonary Disease (COPD)

Metagenomics is the culture-independent study of genetic material obtained directly from environmental samples. It has become a realistic approach to understanding microbial communities thanks to advances in high-throughput DNA sequencing technologies over the past decade. Current research has shown that different sites of the human body house varied bacterial communities. There is a strong correlation between an individual’s microbial community profile at a given site and disease. Metagenomics is being applied more often as a means of comparing microbial profiles in biomedical studies. The analysis of the data collected using metagenomics can be quite challenging and there exist a plethora of tools for interpreting the results. An automatic analytical workflow for metagenomic analyses has been implemented and tested using synthetic datasets of varying quality. It is able to accurately classify bacteria by taxa and correctly estimate the richness and diversity of each set. The workflow was then applied to the study of the airways microbiome in Chronic Obstructive Pulmonary Disease (COPD). COPD is a progressive lung disease resulting in narrowing of the airways and restricted airflow. Despite being the third leading cause of death in the United States, little is known about the differences in the lung microbial community profiles of healthy individuals and COPD patients. Bronchoalveolar lavage (BAL) samples were collected from COPD patients, active or ex-smokers, and never smokers and sequenced by 454 pyrosequencing. A total of 56 individuals were recruited for the study. Substantial colonization of the lungs was found in all subjects and differentially abundant genera in each group were identified. These discoveries are promising and may further our understanding of how the structure of the lung microbiome is modified as COPD progresses. It is also anticipated that the results will eventually lead to improved treatments for COPD.

Genome-wide association of multiple complex traits in outbred mice by ultra low-coverage sequencing

We thank the High-Throughput Genomics Group at the Wellcome Trust Centre for Human Genetics and the Wellcome Trust Sanger Institute for the generation of the sequencing data. This work was funded by Wellcome Trust grant 090532/Z/09/Z (J.F.). Primary phenotyping of the mice was supported by the Mary Lyon Centre and Mammalian Genetics Unit (Medical Research Council, UK Hub grant G0900747 91070 and Medical Research Council, UK grant MC U142684172). D.A.B acknowledges support from NIH R01AR056280. The sleep work was supported by the state of Vaud (Switzerland) and the Swiss National Science Foundation (SNF 14694 and 136201 to P.F.). The ECG work was supported by the Netherlands CardioVascular Research Initiative (Dutch Heart Foundation, Dutch Federation of University Medical Centres, the Netherlands Organization for Health Research and Development, and the Royal Netherlands Academy of Sciences) PREDICT project, InterUniversity Cardiology Institute of the Netherlands (ICIN; 061.02; C.A.R., C.R.B). Na Cai is supported by the Agency of Science, Technology and Research (A*STAR) Graduate Academy. The authors wish to acknowledge excellent technical assistance from: Ayako Kurioka, Leo Swadling, Catherine de Lara, James Ussher, Rachel Townsend, Sima Lionikaite, Ausra S. Lionikiene, Rianne Wolswinkel and Inge van der Made. We would like to thank Thomas M Keane and Anthony G Doran for their help in annotating variants and adding the FVB/NJ strain to the Mouse Genomes Project.

A metagenomic comparison of endemic viruses from broiler chickens with runting stunting syndrome and from normal birds

Runting-stunting syndrome (RSS) in broiler chickens is an enteric disease that causes significant economic losses to poultry producers worldwide due to elevated feed conversion ratios, decreased body weight during growth, and excessive culling. Of specific interest are the viral agents associated with RSS which have been difficult to fully characterise to date. Past research into the aetiology of RSS has implicated a wide variety of RNA and DNA viruses however, to date, no individual virus has been identified as the main agent of RSS and the current opinion is that it may be caused by a community of viruses, collectively known as the virome. This paper attempts to characterise the viral pathogens associated with 2 – 3 week old RSS-affected and unaffected broiler chickens using next-generation sequencing and comparative metagenomics. Analysis of the viromes identified a total of 20 DNA & RNA viral families, along with 2 unidentified categories, comprised of 31 distinct viral genera and 7 unclassified genera. The most abundant viral families identified in this study were the Astroviridae, Caliciviridae, Picornaviridae, Parvoviridae, Coronaviridae, Siphoviridae, and Myoviridae. This study has identified historically significant viruses associated with the disease such as chicken astrovirus, avian nephritis virus, chicken parvovirus, and chicken calicivirus along with relatively novel viruses such as chicken megrivirus and sicinivirus 1 and will help expand the knowledge related to enteric disease in broiler chickens, provide insights into the viral constituents of a healthy avian gut, and identify a variety of enteric viruses and viral communities appropriate for further study.

High diversity of skin-associated bacterial communities of marine fishes is promoted by their high variability among body parts, individuals and species

Animal-associated microbiotas form complex communities, which are suspected to play crucial functions for their host fitness. However, the biodiversity of these communities, including their differences between host species and individuals, has been scarcely studied, especially in case of skin-associated communities. In addition, the intraindividual variability (i.e. between body parts) has never been assessed to date. The objective of this study was to characterize skin bacterial communities of two teleostean fish species, namely the European seabass (Dicentrarchus labrax) and gilthead seabream (Sparus aurata), using a high-throughput DNA sequencing method. In order to focus on intrinsic factors of host-associated bacterial community variability, individuals of the two species were raised in controlled conditions. Bacterial diversity was assessed using a set of four complementary indices, describing the taxonomic and phylogenetic facets of biodiversity and their respective composition (based on presence/absence data) and structure (based on species relative abundances) components. Variability of bacterial diversity was quantified at the interspecific, interindividual and intraindividual scales. We demonstrated that fish surfaces host highly diverse bacterial communities, whose composition was very different from that of surrounding bacterioplankton. This high total biodiversity of skin-associated communities was supported by the important variability, between host species, individuals and the different body parts (dorsal, anal, pectoral and caudal fins).

Sistematic, population genetics, propagation and conservation of Prunus azorica (Rosaceae)

Tese de Doutoramento, Biologia (Ecologia Vegetal), 25 de Junho de 2013, Universidade dos Açores.

Sistematics, population genetics, propagation and conservation of Picconia azorica (Tutin) Knobl. (Oleaceae)

Tese de Doutoramento, Biologia (Ecologia Vegetal), 24 de Junho de 2013, Universidade dos Açores.

Effects of physical exercise training in DNA damage and repair activity in humans with different genetic polymorphisms ofhOGG1(Ser326Cys)

The main purpose of this pilot study was to investigate the possible influence of genetic polymorphisms of the hOGG1 (Ser326Cys) gene in DNA damage and repair activity by 8-oxoguanine DNA glycosylase 1 (OGG1 enzyme) in response to 16 weeks of combined physical exercise training. Thirty-two healthy Caucasian men (40-74 years old) were enrolled in this study. All the subjects were submitted to a training of 16 weeks of combined physical exercise. The subjects with Ser/Ser genotype were considered as wild-type group (WTG), and Ser/Cys and Cys/Cys genotype were analysed together as mutant group (MG). We used comet assay in conjunction with formamidopyrimidine DNA glycoslyase (FPG) to analyse both strand breaks and FPG-sensitive sites. DNA repair activity were also analysed with the comet assay technique. Our results showed no differences between DNA damage (both strand breaks and FPG-sensitive sites) and repair activity (OGG1) between genotype groups (in the pre-training condition). Regarding the possible influence of genotype in the response to 16 weeks of physical exercise training, the results revealed a decrease in DNA strand breaks in both groups, a decrease in FPG-sensitive sites and an increase in total antioxidant capacity in the WTG, but no changes were found in MG. No significant changes in DNA repair activity was observed in both genotype groups with physical exercise training. This preliminary study suggests the possibility of different responses in DNA damage to the physical exercise training, considering the hOGG1 Ser326Cys polymorphism.

Colletotrichum species in Australia

Forty-four species of Colletotrichum are confirmed as present in Australia based on DNA sequencing analyses. Many of these species were identified directly as a result of two workshops organised by the Subcommittee on Plant Health Diagnostics in Australia in 2015 that covered morphological and molecular approaches to identification of Colletotrichum. There are several other species of Colletotrichum reported from Australia that remain to be substantiated by DNA sequence-based methods. This body of work aims to provide a basis from which to critically examine a number of isolates of Colletotrichum deposited in Australian culture collections.

Tecnología Nanopore para la secuenciación de DNA

Nanopore-based sequencer will open the path to the fourth-generation DNA sequencing technology. The main differences between this technique and the previous ones are: DNA molecule that will be sequenced does not need a previous amplification step, is not necessary any type of specific label both molecular adaptors, and it has been abolished enzymatic process in the nucleotide sequence identification event. These differences have as result a more economic method since don’t spend the necessary reagents for the previous techniques, furthermore it lets to sequence samples with a low DNA concentration. This technique is based in the use of a membrane with a biologic nanopore inserted in it whereby the molecule to analyze (analyte) it made to pass, this membrane is placed between two reservoirs containing ions, when an external volatage is applied in both sides this lead to an ion current through the nanopore. When an analyte cross the nanopore the ion current is modified, that modification in the amplitude and duration of ion current determine the physical and chemical properties of that analyte. By means of subsequent statistical analyzes it can be determined to what sequence own this ion current blockade patterns. More used nanopores are the biologic ones, although they are working to develop synthetic nanopores. The main biologic nanopores are: α-Hemolysin from Staphylococcus aureus (α-HL), Mycobacterium smegmatis porin A (MspA) and bacteriophage phi29 pore (phi29). Α-HL and MspA have in their narrowest point a diameter similar to nucleotide size, they are functional at high temperature both wide range of pH (2-12) but MspA is able to read four nucleotide at the same time while α- HL just can read one by one. Finally, phi29 present a bigger diameter what let to get information about DNA spatial conformation and their interaction with proteins (Feng et al., 2015). Nowaday Oxford Nanopore Technologies (ONT) is the only company which has developed Nanopore technology; they have two devices available to sequencing (PromethION and MinION). The MinION is a single-use DNA sequencing device with the size of a USB memory with a total of 3000 nanopores that can sequence until 200kb. The PrometheION is big size sequencer that own 48 different cells, what let to sequence different samples at the same time, with a total of 144.000 nanopores and reading of several megabases (https://www.nanoporetech.com/). The high processivity and low cost become this technique in a great option to massive- sequencing.

Evolution of genome sequencing techniques

The quality and the speed for genome sequencing has advanced at the same time that technology boundaries are stretched. This advancement has been divided so far in three generations. The first-generation methods enabled sequencing of clonal DNA populations. The second-generation massively increased throughput by parallelizing many reactions while the third-generation methods allow direct sequencing of single DNA molecules. The first techniques to sequence DNA were not developed until the mid-1970s, when two distinct sequencing methods were developed almost simultaneously, one by Alan Maxam and Walter Gilbert, and the other one by Frederick Sanger. The first one is a chemical method to cleave DNA at specific points and the second one uses ddNTPs, which synthesizes a copy from the DNA chain template. Nevertheless, both methods generate fragments of varying lengths that are further electrophoresed. Moreover, it is important to say that until the 1990s, the sequencing of DNA was relatively expensive and it was seen as a long process. Besides, using radiolabeled nucleotides also compounded the problem through safety concerns and prevented the automation. Some advancements within the first generation include the replacement of radioactive labels by fluorescent labeled ddNTPs and cycle sequencing with thermostable DNA polymerase, which allows automation and signal amplification, making the process cheaper, safer and faster. Another method is Pyrosequencing, which is based on the “sequencing by synthesis” principle. It differs from Sanger sequencing, in that it relies on the detection of pyrophosphate release on nucleotide incorporation. By the end of the last millennia, parallelization of this method started the Next Generation Sequencing (NGS) with 454 as the first of many methods that can process multiple samples, calling it the 2º generation sequencing. Here electrophoresis was completely eliminated. One of the methods that is sometimes used is SOLiD, based on sequencing by ligation of fluorescently dye-labeled di-base probes which competes to ligate to the sequencing primer. Specificity of the di-base probe is achieved by interrogating every 1st and 2nd base in each ligation reaction. The widely used Solexa/Illumina method uses modified dNTPs containing so called “reversible terminators” which blocks further polymerization. The terminator also contains a fluorescent label, which can be detected by a camera. Now, the previous step towards the third generation was in charge of Ion Torrent, who developed a technique that is based in a method of “sequencing-by-synthesis”. Its main feature is the detection of hydrogen ions that are released during base incorporation. Likewise, the third generation takes into account nanotechnology advancements for the processing of unique DNA molecules to a real time synthesis sequencing system like PacBio; and finally, the NANOPORE, projected since 1995, also uses Nano-sensors forming channels obtained from bacteria that conducts the sample to a sensor that allows the detection of each nucleotide residue in the DNA strand. The advancements in terms of technology that we have nowadays have been so quick, that it makes wonder: ¿How do we imagine the next generation?

Shaping the Dicot Fruit: Molecular and Genomic Approaches to Fruit Development

The fruit is one of the most complex and important structures produced by flowering plants, and understanding the development and maturation process of fruits in different angiosperm species with diverse fruit structures is of immense interest. In the work presented here, molecular genetics and genomic analysis are used to explore the processes that form the fruit in two species: The model organism Arabidopsis and the diploid strawberry Fragaria vesca. One important basic question concerns the molecular genetic basis of fruit patterning. A long-standing model of Arabidopsis fruit (the gynoecium) patterning holds that auxin produced at the apex diffuses downward, forming a gradient that provides apical-basal positional information to specify different tissue types along the gynoecium’s length. The proposed gradient, however, has never been observed and the model appears inconsistent with a number of observations. I present a new, alternative model, wherein auxin acts to establish the adaxial-abaxial domains of the carpel primordia, which then ensures proper development of the final gynoecium. A second project utilizes genomics to identify genes that regulate fruit color by analyzing the genome sequences of Fragaria vesca, a species of wild strawberry. Shared and distinct SNPs among three F. vesca accessions were identified, providing a foundation for locating candidate mutations underlying phenotypic variations among different F. vesca accessions. Through systematic analysis of relevant SNP variants, a candidate SNP in FveMYB10 was identified that may underlie the fruit color in the yellow-fruited accessions, which was subsequently confirmed by functional assays. Our lab has previously generated extensive RNA-sequencing data that depict genome-scale gene expression profiles in F. vesca fruit and flower tissues at different developmental stages. To enhance the accessibility of this dataset, the web-based eFP software was adapted for this dataset, allowing visualization of gene expression in any tissues by user-initiated queries. Together, this thesis work proposes a well-supported new model of fruit patterning in Arabidopsis and provides further resources for F. vesca, including genome-wide variant lists and the ability to visualize gene expression. This work will facilitate future work linking traits of economic importance to specific genes and gaining novel insights into fruit patterning and development.

DNA prime-protein boost based vaccination with a conserved region of leptospiral immunoglobulin-like A and B proteins enhances protection against leptospirosis

Leptospirosis is a zoonotic disease caused by pathogenic spirochetes of the Leptospira genus. Vaccination with bacterins has severe limitations. Here, we evaluated the N-terminal region of the leptospiral immunoglobulin-like B protein (LigBrep) as a vaccine candidate against leptospirosis using immunisation strategies based on DNA primeprotein boost, DNA vaccine, and subunit vaccine. Upon challenge with a virulent strain of Leptospira interrogans , the prime-boost and DNA vaccine approaches induced significant protection in hamsters, as well as a specific IgG antibody response and sterilising immunity. Although vaccination with recombinant fragment of LigBrep also produced a strong antibody response, it was not immunoprotective. These results highlight the potential of LigBrep as a candidate antigen for an effective vaccine against leptospirosis and emphasise the use of the DNA prime-protein boost as an important strategy for vaccine development.

Genetic diversity of Cryptosporidium identified in clinical samples from cities in Brazil and Argentina

The identification and characterisation of Cryptosporidium genotypes and subtypes are fundamental to the study of cryptosporidiosis epidemiology, aiding in prevention and control strategies. The objective was to determine the genetic diversity of Cryptosporidium in samples obtained from hospitals of Rio de Janeiro, Brazil, and Buenos Aires, Argentina. Samples were analysed by microscopy and TaqMan polymerase chain reaction (PCR) assays for Cryptosporidium detection, genotyped by nested-PCR-restriction fragment length polymorphism (RFLP) analysis of the 18S rRNA gene and subtyped by DNA sequencing of the gp60 gene. Among the 89 samples from Rio de Janeiro, Cryptosporidium spp were detected in 26 by microscopy/TaqMan PCR. In samples from Buenos Aires, Cryptosporidium was diagnosed in 15 patients of the 132 studied. The TaqMan PCR and the nested-PCR-RFLP detected Cryptosporidium parvum , Cryptosporidium hominis , and co-infections of both species. In Brazilian samples, the subtypes IbA10G2 and IIcA5G3 were observed. The subtypes found in Argentinean samples were IbA10G2, IaA10G1R4, IaA11G1R4, and IeA11G3T3, and mixed subtypes of Ia and IIa families were detected in the co-infections. C. hominis was the species more frequently detected, and subtype family Ib was reported in both countries. Subtype diversity was higher in Buenos Aires than in Rio de Janeiro and two new subtypes were described for the first time.

Beech bark disease distribution and resistance in Michigan and fungal endophyte ecology of resistant and susceptible beech (Fagus grandifolia Ehrh.)

Beech bark disease (BBD), a non-native association of the fungal pathogen Neonectria faginata and the beech scale insect Cryptococcus fagisuga, has dramatically affected American beech within North American forests. To monitor the spread and effects of BBD in Michigan, a network of forest health monitoring plots was established in 2001 following the disease discovery in Ludington State Park (Mason County). Forest health canopy condition and basic forestry measurements including basal area were reassessed on beech trees in these plots in 2011 and 2012. The influence of bark-inhabiting fungal endophytes on BBD resistance was investigated by collecting cambium tissue from apparently resistant and susceptible beech. Vigor rating showed significant influences of BBD in sample beech resulting in reduced health and substantiated by significant increases of dead beech basal area over time. C. fagisuga distribution was found to be spatially clustered and widespread in the 22 counties in Michigan's Lower Peninsula which contained monitoring plots. Neonectria has been found in Emmet, Cheboygan and Wexford in the Lower Peninsula which may coincide with additional BBD introduction locations. Surveys for BBD resistance resulted in five apparently resistant beech which were added to a BBD resistance database. The most frequently isolated endophytes from cambium tissue were identified by DNA sequencing primarily as Deuteromycetes and Ascomycetes including Chaetomium globosum, Neohendersonia kickxii and Fusarium flocciferum. N. faginata in antagonism trials showed significant growth reduction when paired with three beech fungal endophytes. The results of the antagonism trial and decay tests indicate that N. faginata may be a relatively poor competitor in vivo with limited ability to degrade cellulose.

Understanding ten-eleven translocation-2 in hematological and nervous systems

I proposed the study of two distinct aspects of Ten-Eleven Translocation 2 (TET2) protein for understanding specific functions in different body systems. ^ In Part I, I characterized the molecular mechanisms of Tet2 in the hematological system. As the second member of Ten-Eleven Translocation protein family, TET2 is frequently mutated in leukemic patients. Previous studies have shown that the TET2 mutations frequently occur in 20% myelodysplastic syndrome/myeloproliferative neoplasm (MDS/MPN), 10% T-cell lymphoma leukemia and 2% B-cell lymphoma leukemia. Genetic mouse models also display distinct phenotypes of various types of hematological malignancies. I performed 5-hydroxymethylcytosine (5hmC) chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA sequencing (RNA-Seq) of hematopoietic stem/progenitor cells to determine whether the deletion of Tet2 can affect the abundance of 5hmC at myeloid, T-cell and B-cell specific gene transcription start sites, which ultimately result in various hematological malignancies. Subsequent Exome sequencing (Exome-Seq) showed that disease-specific genes are mutated in different types of tumors, which suggests that TET2 may protect the genome from being mutated. The direct interaction between TET2 and Mutator S Homolog 6 (MSH6) protein suggests TET2 is involved in DNA mismatch repair. Finally, in vivo mismatch repair studies show that the loss of Tet2 causes a mutator phenotype. Taken together, my data indicate that TET2 binds to MSH6 to protect genome integrity. ^ In Part II, I intended to better understand the role of Tet2 in the nervous system. 5-hydroxymethylcytosine regulates epigenetic modification during neurodevelopment and aging. Thus, Tet2 may play a critical role in regulating adult neurogenesis. To examine the physiological significance of Tet2 in the nervous system, I first showed that the deletion of Tet2 reduces the 5hmC levels in neural stem cells. Mice lacking Tet2 show abnormal hippocampal neurogenesis along with 5hmC alternations at different gene promoters and corresponding gene expression downregulation. Through the luciferase reporter assay, two neural factors Neurogenic differentiation 1 (NeuroD1) and Glial fibrillary acidic protein (Gfap) were down-regulated in Tet2 knockout cells. My results suggest that Tet2 regulates neural stem/progenitor cell proliferation and differentiation in adult brain.^