Whole genome sequencing of enriched chloroplast DNA using the Illumina GAII platform

<p>Abstract</p> <p>Background</p> <p>Complete chloroplast genome sequences provide a valuable source of molecular markers for studies in molecular ecology and evolution of plants. To obtain complete genome sequences, recent studies have made use of the polymerase chain reaction to amplify overlapping fragments from conserved gene loci. However, this approach is time consuming and can be more difficult to implement where gene organisation differs among plants. An alternative approach is to first isolate chloroplasts and then use the capacity of high-throughput sequencing to obtain complete genome sequences. We report our findings from studies of the latter approach, which used a simple chloroplast isolation procedure, multiply-primed rolling circle amplification of chloroplast DNA, Illumina Genome Analyzer II sequencing, and de novo assembly of paired-end sequence reads.</p> <p>Results</p> <p>A modified rapid chloroplast isolation protocol was used to obtain plant DNA that was enriched for chloroplast DNA, but nevertheless contained nuclear and mitochondrial DNA. Multiply-primed rolling circle amplification of this mixed template produced sufficient quantities of chloroplast DNA, even when the amount of starting material was small, and improved the template quality for Illumina Genome Analyzer II (hereafter Illumina GAII) sequencing. We demonstrate, using independent samples of karaka (<it>Corynocarpus laevigatus</it>), that there is high fidelity in the sequence obtained from this template. Although less than 20% of our sequenced reads could be mapped to chloroplast genome, it was relatively easy to assemble complete chloroplast genome sequences from the mixture of nuclear, mitochondrial and chloroplast reads.</p> <p>Conclusions</p> <p>We report successful whole genome sequencing of chloroplast DNA from karaka, obtained efficiently and with high fidelity.</p

The Association Between Vitamin D and Multiple Sclerosis Risk: 1,25(OH)2D3 Induces Super-Enhancers Bound by VDR

A super-enhancer (SE) is a cluster of enhancers with a relatively high density of particular chromatin features. SEs typically regulate key genes that can determine cell identity and differentiation. Identifying SEs and their effects may be critical in predicting key regulatory genes, such as master transcription factor genes or oncogenes. Signal inducible SEs are dense stretches of signal terminal transcription factor (TF) binding regions, and may modulate the interaction between environmental factors (e.g., Vitamin D) and genetic factors (i.e., risk variants) in complex diseases such as multiple sclerosis (MS). As a complex autoimmune disease, the etiology and progression of MS, including the interaction between Vitamin D and MS risk variants, is still unclear and can be explored from the aspect of signal SEs. Vitamin D [with its active form: 1,25(OH)2D3], is an environmental risk factor for MS. It binds the Vitamin D receptor (VDR) and regulates gene expression. This study explores the association between VDR super-enhancers (VSEs) and MS risk variants. Firstly, we reanalyse public ChIP-seq and RNA-seq data to classify VSEs into three categories according to their combinations of persistent and secondary VDR binding. Secondly, we indicate the genes with VSE regions that are near MS risk variants. Furthermore, we find that MS risk variants are enriched in VSE regions, and we indicate some genes with a VSE overlapping MS risk variant for further exploration. We also find two clusters of genes from the set of genes showing correlation of expression patterns with the MS risk gene ZMIZ1 that appear to be regulated by VSEs in THP-1 cells. It is the first time that VSEs have been analyzed, and we directly connect the genetic risk factors for MS risk with Vitamin D based on VSEs

Association of Genetic Variation with Keratoconus

Importance: Keratoconus is a condition in which the cornea progressively thins and protrudes in a conical shape, severely affecting refraction and vision. It is a major indication for corneal transplant. To discover new genetic loci associated with keratoconus and better understand the causative mechanism of this disease, we performed a genome-wide association study on patients with keratoconus.Objective: To identify genetic susceptibility regions for keratoconus in the human genome.Design, Setting, and Participants: This study was conducted with data from eye clinics in Australia, the United States, and Northern Ireland. The discovery cohort of individuals with keratoconus and control participants from Australia was genotyped using the Illumina HumanCoreExome single-nucleotide polymorphism array. After quality control and data cleaning, genotypes were imputed against the 1000 Genomes Project reference panel (phase III; version 5), and association analyses were completed using PLINK version 1.90. Single-nucleotide polymorphisms with P -6 were assessed for replication in 3 additional cohorts. Control participants were drawn from the cohorts of the Blue Mountains Eye Study and a previous study of glaucoma. Replication cohorts were from a previous keratoconus genome-wide association study data set from the United States, a cohort of affected and control participants from Australia and Northern Ireland, and a case-control cohort from Victoria, Australia. Data were collected from January 2006 to March 2019.Main Outcomes and Measures: Associations between keratoconus and 6 252 612 genetic variants were estimated using logistic regression after adjusting for ancestry using the first 3 principal components.Results: The discovery cohort included 522 affected individuals and 655 control participants, while the replication cohorts included 818 affected individuals (222 from the United States, 331 from Australia and Northern Ireland, and 265 from Victoria, Australia) and 3858 control participants (2927 from the United States, 229 from Australia and Northern Ireland, and 702 from Victoria, Australia). Two novel loci reached genome-wide significance (defined as P -8), with a P value of 7.46 × 10-9 at rs61876744 in patatin-like phospholipase domain-containing 2 gene (PNPLA2) on chromosome 11 and a P value of 6.35 × 10-12 at rs138380, 2.2 kb upstream of casein kinase I isoform epsilon gene (CSNK1E) on chromosome 22. One additional locus was identified with a P value less than 1.00 × 10-6 in mastermind-like transcriptional coactivator 2 (MAML2) on chromosome 11 (P = 3.91 × 10-7). The novel locus in PNPLA2 reached genome-wide significance in an analysis of all 4 cohorts (P = 2.45 × 10-8).Conclusions and Relevance: In this relatively large keratoconus genome-wide association study, we identified a genome-wide significant locus for keratoconus in the region of PNPLA2 on chromosome 11

Evolutionary relationships and divergence times among the native rats of Australia

Background The genus Rattus is highly speciose and has a complex taxonomy that is not fully resolved. As shown previously there are two major groups within the genus, an Asian and an Australo-Papuan group. This study focuses on the Australo-Papuan group and particularly on the Australian rats. There are uncertainties regarding the number of species within the group and the relationships among them. We analysed 16 mitochondrial genomes, including seven novel genomes from six species, to help elucidate the evolutionary history of the Australian rats. We also demonstrate, from a larger dataset, the usefulness of short regions of the mitochondrial genome in identifying these rats at the species level. Results Analyses of 16 mitochondrial genomes representing species sampled from Australo-Papuan and Asian clades of Rattus indicate divergence of these two groups ~2.7 million years ago (Mya). Subsequent diversification of at least 4 lineages within the Australo-Papuan clade was rapid and occurred over the period from ~ 0.9-1.7 Mya, a finding that explains the difficulty in resolving some relationships within this clade. Phylogenetic analyses of our 126 taxon, but shorter sequence (1952 nucleotides long), Rattus database generally give well supported species clades. Conclusions Our whole mitochondrial genome analyses are concordant with a taxonomic division that places the native Australian rats into the Rattus fuscipes species group. We suggest the following order of divergence of the Australian species. R. fuscipes is the oldest lineage among the Australian rats and is not part of a New Guinean radiation. R. lutreolus is also within this Australian clade and shallower than R. tunneyi while the R. sordidus group is the shallowest lineage in the clade. The divergences within the R. sordidus and R. leucopus lineages occurring about half a million years ago support the hypotheses of more recent interchanges of rats between Australia and New Guinea. While problematic for inference of deeper divergences, we report that the analysis of shorter mitochondrial sequences is very useful for species identification in rats

A multi-ethnic genome-wide association study implicates collagen matrix integrity and cell differentiation pathways in keratoconus

Keratoconus is characterised by reduced rigidity of the cornea with distortion and focal thinning that causes blurred vision, however, the pathogenetic mechanisms are unknown. It can lead to severe visual morbidity in children and young adults and is a common indication for corneal transplantation worldwide. Here we report the first large scale genome-wide association study of keratoconus including 4,669 cases and 116,547 controls. We have identified significant association with 36 genomic loci that, for the first time, implicate both dysregulation of corneal collagen matrix integrity and cell differentiation pathways as primary disease-causing mechanisms. The results also suggest pleiotropy, with some disease mechanisms shared with other corneal diseases, such as Fuchs endothelial corneal dystrophy. The common variants associated with keratoconus explain 12.5% of the genetic variance, which shows potential for the future development of a diagnostic test to detect susceptibility to disease

Comparing vision and macular thickness in neovascular age-related macular degeneration, diabetic macular oedema and retinal vein occlusion patients treated with intravitreal antivascular endothelial growth factor injections in clinical practice

Objective To compare the visual outcomes of intravitreal antivascular endothelial growth factor (anti-VEGF) injections in neovascular age-related macular degeneration (nAMD), diabetic macular oedema (DMO) and retinal vein occlusion (RVO) in a real-world setting.Methods and analysis Retrospective analysis of data from the Tasmanian Ophthalmic Biobank database. The median change in best-corrected visual acuity (BCVA) between baseline and 12 months post initiating intravitreal anti-VEGF treatment were compared between the three diseases. Final BCVA, central macular thickness (CMT), cumulative number of injections and overall predictors of change in BCVA and CMT were also determined.Results At 12 months, change in BCVA was significantly different between nAMD, DMO and RVO cohorts (p=0.032), with lower median change for DMO (2 letters, range −5 to 20) than for RVO (11 letters, range −20 to 35). Likewise, CMT change was significantly different between the three cohorts (p=0.022), with a smaller reduction in CMT in DMO (−54 µm, range −482 to 50) than RVO patients (−137 µm, range −478 to 43; p=0.033). Total number of injections received (p=0.028) and final BCVA score (p=0.024) were also significantly different between the groups. Baseline BCVA was a negative predictor (p=0.042) and baseline CMT a positive predictor (p&lt;0.001) of outcome. After adjusting for baseline BCVA and CMT, diagnosis of nAMD or RVO was a predictor of visual improvement compared with the DMO.Conclusions At the end of 12 months, nAMD and RVO cohorts had the greatest improvement in BCVA, however the final BCVA for DMO was significantly better than for nAMD

The complete mitochondrial genome of the eastern grey kangaroo (Macropus giganteus)

We present the complete mitochondrial genome (accession number: LK995454) of an iconic Australian species, the eastern grey kangaroo (Macropus giganteus). The mitogenomic organization is consistent with other marsupials, encoding 13 protein-coding genes, 22 tRNA genes, 2 ribosomal RNA genes, an origin of light strand replication and a control region or Dloop. No repetitive sequences were detected in the control region. The M. giganteus mitogenome exemplifies a combination of tRNA gene order and structural peculiarities that appear to be unique to marsupials. We present a maximum likelihood phylogeny based on complete mitochondrial protein and RNA coding sequences that confirms the phylogenetic position of the grey kangaroo among macropodids