Assessing the significance of conserved genomic aberrations using high resolution genomic microarrays

Genomic aberrations recurrent in a particular cancer type can be important prognostic markers for tumor progression. Typically in early tumorigenesis, cells incur a breakdown of the DNA replication machinery that results in an accumulation of genomic aberrations in the form of duplications, deletions, translocations, and other genomic alterations. Microarray methods allow for finer mapping of these aberrations than has previously been possible; however, data processing and analysis methods have not taken full advantage of this higher resolution. Attention has primarily been given to analysis on the single sample level, where multiple adjacent probes are necessarily used as replicates for the local region containing their target sequences. However, regions of concordant aberration can be short enough to be detected by only one, or very few, array elements. We describe a method called Multiple Sample Analysis for assessing the significance of concordant genomic aberrations across multiple experiments that does not require a-priori definition of aberration calls for each sample. If there are multiple samples, representing a class, then by exploiting the replication across samples our method can detect concordant aberrations at much higher resolution than can be derived from current single sample approaches. Additionally, this method provides a meaningful approach to addressing population-based questions such as determining important regions for a cancer subtype of interest or determining regions of copy number variation in a population. Multiple Sample Analysis also provides single sample aberration calls in the locations of significant concordance, producing high resolution calls per sample, in concordant regions. The approach is demonstrated on a dataset representing a challenging but important resource: breast tumors that have been formalin-fixed, paraffin-embedded, archived, and subsequently UV-laser capture microdissected and hybridized to two-channel BAC arrays using an amplification protocol. We demonstrate the accurate detection on simulated data, and on real datasets involving known regions of aberration within subtypes of breast cancer at a resolution consistent with that of the array. Similarly, we apply our method to previously published datasets, including a 250K SNP array, and verify known results as well as detect novel regions of concordant aberration. The algorithm has been fully implemented and tested and is freely available as a Java application at http://www.cbil.upenn.edu/MSA.     

A chimeric protein of simian immunodeficiency virus envelope glycoprotein gp140 and Escherichia coli aspartate transcarbamoylase

The envelope glycoproteins of the human immunodeficiency virus and the related simian immunodeficiency virus (SIV) mediate viral entry into host cells by fusing viral and target cell membranes. We have reported expression, purification, and characterization of gp140 (also called gp160e), the soluble, trimeric ectodomain of the SIV envelope glycoprotein, gp160 (B. Chen et al., J. Biol. Chem. 275:34946-34953, 2000). We have now expressed and purified chimeric proteins of SIV gp140 and its variants with the catalytic subunit (C) of Escherichia coli aspartate transcarbamoylase (ATCase). The fusion proteins (SIV gp140-ATC) bind viral receptor CD4 and a number of monoclonal antibodies specific for SIV gp140. The chimeric molecule also has ATCase activity, which requires trimerization of the ATCase C chains. Thus, the fusion protein is trimeric. When ATCase regulatory subunit dimers (R(2)) are added, the fusion protein assembles into dimers of trimers as expected from the structure of C(6)R(6) ATCase. Negative-stain electron microscopy reveals spikey features of both SIV gp140 and SIV gp140-ATC. The production of the fusion proteins may enhance the possibilities for structure determination of the envelope glycoprotein either by electron cryomicroscopy or X-ray crystallography.     

Biochemical Conservation and Evolution of Germacrene A Oxidase in Asteraceae

Sesquiterpene lactones are characteristic natural products in Asteraceae, which constitutes ∼8% of all plant species. Despite their physiological and pharmaceutical importance, the biochemistry and evolution of sesquiterpene lactones remain unexplored. Here we show that germacrene A oxidase (GAO), evolutionarily conserved in all major subfamilies of Asteraceae, catalyzes three consecutive oxidations of germacrene A to yield germacrene A acid. Furthermore, it is also capable of oxidizing non-natural substrate amorphadiene. Co-expression of lettuce GAO with germacrene synthase in engineered yeast synthesized aberrant products, costic acids and ilicic acid, in an acidic condition. However, cultivation in a neutral condition allowed the de novo synthesis of a single novel compound that was identified as germacrene A acid by gas and liquid chromatography and NMR analyses. To trace the evolutionary lineage of GAO in Asteraceae, homologous genes were further isolated from the representative species of three major subfamilies of Asteraceae (sunflower, chicory, and costus from Asteroideae, Cichorioideae, and Carduoideae, respectively) and also from the phylogenetically basal species, Barnadesia spinosa, from Barnadesioideae. The recombinant GAOs from these genes clearly showed germacrene A oxidase activities, suggesting that GAO activity is widely conserved in Asteraceae including the basal lineage. All GAOs could catalyze the three-step oxidation of non-natural substrate amorphadiene to artemisinic acid, whereas amorphadiene oxidase diverged from GAO displayed negligible activity for germacrene A oxidation. The observed amorphadiene oxidase activity in GAOs suggests that the catalytic plasticity is embedded in ancestral GAO enzymes that may contribute to the chemical and catalytic diversity in nature.     

RAKE and LNA-ISH reveal microRNA expression and localization in archival human brain

microRNAs (miRNAs) are small (approximately 22 nucleotide) regulatory RNAs which play fundamental roles in many biological processes. Recent studies have shown that the expression of many miRNAs is altered in various human tumors and some miRNAs may function as oncogenes or tumor suppressor genes. However, with the exception of glioblastoma multiforme, the expression of miRNAs in brain tumors is unknown. Furthermore, methods to profile miRNAs from formalin-fixed, paraffin-embedded (FFPE) archival tissues or to study their cellular and subcellular localization in FFPE tissues have been lacking. Here we report the coordinated miRNA expression analysis from the tissue level to the subcellular level, using the RAKE (RNA-primed, array-based, Klenow Enzyme) miRNA microarray platform in conjunction with Locked Nucleic Acid (LNA)-based in situ hybridization (LNA-ISH) on archival FFPE human brains and oligodendroglial tumors. The ability to profile miRNAs from archival tissues at the tissue level, by RAKE microarrays, and at the cellular level by LNA-ISH, will accelerate studies of miRNAs in human diseases.     

QTL analysis of cotton fiber length in advanced backcross populations derived from a cross between <Emphasis Type="Italic">Gossypium hirsutum</Emphasis> and <Emphasis Type="Italic">G. mustelinum</Emphasis>

Key message

QTLs for fiber length mapped in three generations of advanced backcross populations derived from crossing Gossypium hirsutum and Gossypium mustelinum showed opportunities to improve elite cottons by introgression from wild relatives.

Abstract

The molecular basis of cotton fiber length in crosses between Gossypium hirsutum and Gossypium mustelinum was dissected using 21 BC₃F₂ and 12 corresponding BC₃F_2:3 and BC₃F_2:4 families. Sixty-five quantitative trait loci (QTLs) were detected by one-way analysis of variance. The QTL numbers detected for upper-half mean length (UHM), fiber uniformity index (UI), and short fiber content (SFC) were 19, 20, and 26 respectively. Twenty-three of the 65 QTLs could be detected at least twice near adjacent markers in the same family or near the same markers across different families/generations, and 32 QTLs were detected in both one-way variance analyses and mixed model-based composite interval mapping. G. mustelinum alleles increased UHM and UI and decreased SFC for five, one, and one QTLs, respectively. In addition to the main-effect QTLs, 17 epistatic QTLs were detected which helped to elucidate the genetic basis of cotton fiber length. Significant among-family genotypic effects were detected at 18, 16, and 16 loci for UHM, UI, and SFC, respectively. Six, two, and two loci showed genotype?×?family interaction for UHM, UI and SFC, respectively, illustrating complexities that might be faced in introgression of exotic germplasm into cultivated cotton. Co-location of many QTLs for UHM, UI, and SFC accounted for correlations among these traits, and selection of these QTLs may improve the three traits simultaneously. The simple sequence repeat (SSR) markers associated with G. mustelinum QTLs will assist breeders in transferring and maintaining valuable traits from this exotic source during cultivar development.

     

Management of threatened insect species in Australia, with particular reference to butterflies

Abstract   Advances in determination of conservation status of insects through evaluation of threats (rather than simple 'rarity' alone) help to set priorities for allocation of conservation resources, and formulate effective legislative treatment and field management for the most needy species. Summaries of several recent conservation studies on butterflies exemplify the variety of practical needs, amongst which habitat protection and enhancement is paramount. However, the paucity of biological and distributional information on most Australian insects necessitates considerable flexibility in approach, so that perceptions of conservation need may be highly subjective. Listing species for protection and conservation priority must be viewed as a dynamic process, with opportunity for rapid action, either to list or de-list taxa as research and recovery actions lead to change in perceived or real status and conservation need.     

Differential expression of mycothiol pathway genes: are they affected by antituberculosis drugs?

Mycothiol (MSH) is the major cellular thiol in Mycobacterium tuberculosis (M.tb). We hypothesize that the mycothiol-dependent detoxification pathway may serve an important role during oxygen stress management in M. tuberculosis, derived from normal aerobic metabolism, the macrophage environment and through the action of anti-tubercular antibiotics, such as Isoniazid (INH). Total mRNA and DNA were isolated from M. bovis BCG at different stages of growth in 7H9 mycobacterial medium. Three genes involved in mycothiol metabolism and encoding the enzymes mycothiol S-conjugate amidase (Mca, Rv1082), NADPH dependent mycothiol reductase (mtr, Rv2855), and N-Acetyl-1-D-myo-Inosityl-2-Amino-2-Deoxy-alpha-D-Glucopyranoside Deacetylase (GlcNAc-Ins deacetylase, Rv1170 or mshB) were investigated for genomic rearrangements and expression. The results show that the genomic domains of the genes remain conserved in evolutionary diverse and unrelated M. tuberculosis isolates. The genes encoding enzymes implicated in mycothiol reduction, mtr (Rv2855) and the mycothiol-dependant detoxification of electrophilic agents, Mca (Rv1082), are shown to be actively transcribed during logarithmic M. bovis BCG growth. The gene encoding GlcNAc-Ins deacetylase (the rate limiting mycothiol biosynthesis step) shows induction in the presence of INH. Antisense oligonucleotides to both GlcNAc-Ins deacetylase (Rv1170) and mtr (Rv2855) mRNA affect mycobacterial growth. In conclusion the results presented here suggest that these enzymes are sensitive to free radical generating antituberculosis drugs and may be useful targets for new drug development.     

A cross-continental test of the Enemy Release Hypothesis: leaf herbivory on Acer platanoides (L.) is three times lower in North America than in its native Europe

Acer platanoides (Norway maple) is a widespread native tree species in Europe. It has been introduced to North America where it has often established dense stands in both secondary woodlands and relatively undisturbed mature woodlands. In Europe A. platanoides is also extending its original range, but generally seems to exist at much lower densities. One explanation for the ‘aggressiveness’ of invasive plants such as A. platanoides is that they have left behind pests and diseases which limit their population densities in their native lands (the enemy release hypothesis or ERH). To assess the ERH for Norway maple, a large network of collaborators assessed leaf herbivory rates in populations throughout Europe and North America. We found significantly lower total leaf herbivory (1.6% ± 0.19, n = 21 vs. 7.4% ± 1.94, n = 34) and lower fungal damage (1.0% ± 0.35, n = 13 vs. 3.7% ± 0.85, n = 34) in North America than in Europe over a 2 year period, which is consistent with the predictions of the Enemy Release Hypothesis. Across years, the average total leaf herbivory was significantly correlated with average annual temperature of the site (P < 0.05), although this was mostly due to sites in Europe (P < 0.001), and not sites in North America (P > 0.05). Furthermore, only populations in Europe showed very high levels of herbivory (e.g., nine sites had total leaf herbivory ranging from 10.0 to 51.2% in at least 1 year) or leaf fungal damage (only one site in North America showed high levels of fungal damage in 1 year), suggesting the possibility of more frequent episodic outbreaks in the native range. Leaf herbivory and fungal damage are only two aspects of consumer pressure and we do not know whether the differences reported here are enough to actually elicit release from top-down population control, but such large scale biogeographic differences in herbivory contribute towards understanding exotic invasions. Jonathan M. Adams and Wei Fang—equally contributed as first authors. A list of the participating members of the Transatlantic Acer platanoides Invasion Network is given in the Appendix 3.     

Neural Population Structures and Consequences for Neural Coding

Researchers studying neural coding have speculated that populations of neurons would more effectively represent the stimulus if the neurons "cooperated:" by interacting through lateral connections, the neurons would process and represent information better than if they functioned independently. We apply our new theory of information processing to determine the fidelity limits of simple population structures to encode stimulus features. We focus on noncooperative populations, which have no lateral connections. We show that they always exhibit positively correlated responses and that as population size increases, they perfectly represent the information conveyed by their inputs regardless of the individual neuron's coding scheme. Cooperative populations, which do have lateral connections, can, depending on the nature of the connections, perform better or worse than their noncooperative counterparts. We further show that common notions of synergy fail to capture the level of cooperation and to reflect the information processing properties of populations.     

Whole Exome Re-Sequencing Implicates CCDC38 and Cilia Structure and Function in Resistance to Smoking Related Airflow Obstruction

Chronic obstructive pulmonary disease (COPD) is a leading cause of global morbidity and mortality and, whilst smoking remains the single most important risk factor, COPD risk is heritable. Of 26 independent genomic regions showing association with lung function in genome-wide association studies, eleven have been reported to show association with airflow obstruction. Although the main risk factor for COPD is smoking, some individuals are observed to have a high forced expired volume in 1 second (FEV₁) despite many years of heavy smoking. We hypothesised that these “resistant smokers” may harbour variants which protect against lung function decline caused by smoking and provide insight into the genetic determinants of lung health. We undertook whole exome re-sequencing of 100 heavy smokers who had healthy lung function given their age, sex, height and smoking history and applied three complementary approaches to explore the genetic architecture of smoking resistance. Firstly, we identified novel functional variants in the “resistant smokers” and looked for enrichment of these novel variants within biological pathways. Secondly, we undertook association testing of all exonic variants individually with two independent control sets. Thirdly, we undertook gene-based association testing of all exonic variants. Our strongest signal of association with smoking resistance for a non-synonymous SNP was for rs10859974 (P = 2.34×10⁻⁴) in CCDC38, a gene which has previously been reported to show association with FEV₁/FVC, and we demonstrate moderate expression of CCDC38 in bronchial epithelial cells. We identified an enrichment of novel putatively functional variants in genes related to cilia structure and function in resistant smokers. Ciliary function abnormalities are known to be associated with both smoking and reduced mucociliary clearance in patients with COPD. We suggest that genetic influences on the development or function of cilia in the bronchial epithelium may affect growth of cilia or the extent of damage caused by tobacco smoke.