Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history.

Author Summary

Arthropods are the most abundant animals on earth. Among them, insects clearly dominate on land, whereas crustaceans hold the title for the most diverse invertebrates in the oceans. Much is known about the biology of these groups, not least because of genomic studies of the fruit fly Drosophila, the water flea Daphnia, and other species used in research. Here we report the first genome sequence from a species belonging to a lineage that has previously received very little attention—the myriapods. Myriapods were among the first arthropods to invade the land over 400 million years ago, and survive today as the herbivorous millipedes and venomous centipedes, one of which—Strigamia maritima—we have sequenced here. We find that the genome of this centipede retains more characteristics of the presumed arthropod ancestor than other sequenced insect genomes. The genome provides access to many aspects of myriapod biology that have not been studied before, suggesting, for example, that they have diversified receptors for smell that are quite different from those used by insects. In addition, it shows specific consequences of the largely subterranean life of this particular species, which seems to have lost the genes for all known light-sensing molecules, even though it still avoids light.

There has been a renewed interest in investigating the role of stabilizing selection acting on genome-wide traits such as codon usage bias. Codon bias, when synonymous codons are used at unequal frequencies, occurs in a wide variety of taxa. Standard evolutionary models explain the maintenance of codon bias through a balance of genetic drift, mutation and weak purifying selection. The efficacy of selection is expected to be reduced in regions of suppressed recombination. Contrary to observations in Drosophila melanogaster, some recent studies have failed to detect a relationship between the recombination rate, intensity of selection acting at synonymous sites, and the magnitude of codon bias as predicted under these standard models. Here, we examined codon bias in 2798 protein coding loci on the third chromosome of D. pseudoobscura using whole-genome sequences of 47 individuals, representing five common third chromosome gene arrangements. Fine-scale recombination maps were constructed using more than 1 million segregating sites. As expected, recombination was demonstrated to be significantly suppressed between chromosome arrangements, allowing for a direct examination of the relationship between recombination, selection, and codon bias. As with other Drosophila species, we observe a strong mutational bias away from the most frequently used codons. We find the rate of synonymous and nonsynonymous polymorphism is variable between different amino acids. However, we do not observe a reduction in codon bias or the strength of selection in regions of suppressed recombination as expected. Instead, we find that the interaction between weak stabilizing selection and mutational bias likely plays a role in shaping the composition of synonymous codons across the third chromosome in D. pseudoobscura.

Haplogroup (hg) H dominates present-day Western European mitochondrial (mt) DNA variability (>40%), yet was less common (~19%) amongst Early Neolithic farmers (~5450 BC) and virtually absent in Mesolithic hunter-gatherers. Here we investigate this major component of the maternal population history of modern Europeans and sequence 39 complete hg H mitochondrial genomes from ancient human remains. We then compare this ‘real-time’ genetic data with cultural changes taking place between the Early Neolithic (~5450 BC) and Bronze Age (~2200 BC) in Central Europe. Our results reveal that the current diversity and distribution of hg H were largely established by the Mid-Neolithic (~4000 BC), but with substantial genetic contributions from subsequent pan-European cultures such as the Bell Beakers expanding out of Iberia in the Late Neolithic (~2800 BC). Dated hg H genomes allow us to reconstruct the recent evolutionary history of hg H and reveal a mutation rate 45% higher than current estimates for human mitochondria.

Lubricin may be an important barrier to the development of corneal and conjunctival epitheliopathies that may occur in dry eye disease and contact lens wear.

Objective

To test the hypotheses that lubricin (ie, proteoglycan 4 [PRG4]), a boundary lubricant, is produced by ocular surface epithelia and acts to protect the cornea and conjunctiva against significant shear forces generated during an eyelid blink and that lubricin deficiency increases shear stress on the ocular surface and promotes corneal damage.

Design, Setting, and Participants

Human, porcine, and mouse tissues and cells were processed for molecular biological, immunohistochemical, and tribological studies, and wild-type and PRG4 knockout mice were evaluated for corneal damage.

Results

Our findings demonstrate that lubricin is transcribed and translated by corneal and conjunctival epithelial cells. Lubricin messenger RNA is also present in lacrimal and meibomian glands, as well as in a number of other tissues. Absence of lubricin in PRG4 knockout mice is associated with a significant increase in corneal fluorescein staining. Our studies also show that lubricin functions as an effective friction-lowering boundary lubricant at the human cornea-eyelid interface. This effect is specific and cannot be duplicated by the use of hyaluronate or bovine serum albumin solutions.

Conclusions and Relevance

Our results show that lubricin is transcribed, translated, and expressed by ocular surface epithelia. Moreover, our findings demonstrate that lubricin presence significantly reduces friction between the cornea and conjunctiva and that lubricin deficiency may play a role in promoting corneal damage.

We hypothesize that aromatase, an enzyme that controls estrogen biosynthesis, plays a major role in the sex-related differences of the meibomian gland. To begin to test this hypothesis, we examined the influence of aromatase absence, which completely eliminates estrogen production, on glandular gene expression and histology in male and female mice.

Methods.

Meibomian glands were obtained from adult, age-matched wild-type (WT) and aromatase knockout (ArKO) mice. Tissues were processed for histology or the isolation of total RNA, which was analyzed for differentially expressed mRNAs by using microarrays.

Results.

Our results show that aromatase significantly influences the expression of more than a thousand genes in the meibomian gland. The nature of this effect is primarily sex-dependent. In addition, the influence of aromatase on sex-related differences in gene expression is predominantly genotype-specific. However, many of the sex-related variations in biological process, molecular function, and cellular component ontologies, as well as in KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways, are remarkably similar between WT and ArKO mice. The loss of aromatase activity has no obvious effect on the histology of meibomian glands in male or female mice.

Conclusions.

Our findings demonstrate that aromatase has a significant impact on gene expression in the meibomian gland. The nature of this influence is sex-dependent and genotype-specific; however, many of the sex-related variations in gene ontologies and KEGG pathways are similar between WT and ArKO mice. Consequently, it appears that aromatase, and by extension estrogen, do not play a major role in the sex-related differences of the mouse meibomian gland.

Aromatase exerts a significant sex- and genotype-specific effect on meibomian gland gene expression. However, this enzyme, and by extension estrogen, do not play a major role in the sex-related differences in this tissue.

The process of generating raw genome sequence data continues to become cheaper, faster, and more accurate. However, assembly of such data into high-quality, finished genome sequences remains challenging. Many genome assembly tools are available, but they differ greatly in terms of their performance (speed, scalability, hardware requirements, acceptance of newer read technologies) and in their final output (composition of assembled sequence). More importantly, it remains largely unclear how to best assess the quality of assembled genome sequences. The Assemblathon competitions are intended to assess current state-of-the-art methods in genome assembly.

Results

In Assemblathon 2, we provided a variety of sequence data to be assembled for three vertebrate species (a bird, a fish, and snake). This resulted in a total of 43 submitted assemblies from 21 participating teams. We evaluated these assemblies using a combination of optical map data, Fosmid sequences, and several statistical methods. From over 100 different metrics, we chose ten key measures by which to assess the overall quality of the assemblies.

Conclusions

Many current genome assemblers produced useful assemblies, containing a significant representation of their genes and overall genome structure. However, the high degree of variability between the entries suggests that there is still much room for improvement in the field of genome assembly and that approaches which work well in assembling the genome of one species may not necessarily work well for another.

A major challenge of biology is understanding the relationship between molecular genetic variation and variation in quantitative traits, including fitness. This relationship determines our ability to predict phenotypes from genotypes and to understand how evolutionary forces shape variation within and between species. Previous efforts to dissect the genotype-phenotype map were based on incomplete genotypic information. Here, we describe the Drosophila melanogaster Genetic Reference Panel (DGRP), a community resource for analysis of population genomics and quantitative traits. The DGRP consists of fully sequenced inbred lines derived from a natural population. Population genomic analyses reveal reduced polymorphism in centromeric autosomal regions and the X chromosome, evidence for positive and negative selection, and rapid evolution of the X chromosome. Many variants in novel genes, most at low frequency, are associated with quantitative traits and explain a large fraction of the phenotypic variance. The DGRP facilitates genotype-phenotype mapping using the power of Drosophila genetics.

MicroRNAs (miRNAs) are small non-coding RNAs that play critical roles in regulating post transcriptional gene expression. Gall midges encompass a large group of insects that are of economic importance and also possess fascinating biological traits. The gall midge Mayetiola destructor, commonly known as the Hessian fly, is a destructive pest of wheat and model organism for studying gall midge biology and insect – host plant interactions.

Results

In this study, we systematically analyzed miRNAs from the Hessian fly. Deep-sequencing a Hessian fly larval transcriptome led to the identification of 89 miRNA species that are either identical or very similar to known miRNAs from other insects, and 184 novel miRNAs that have not been reported from other species. A genome-wide search through a draft Hessian fly genome sequence identified a total of 611 putative miRNA-encoding genes based on sequence similarity and the existence of a stem-loop structure for miRNA precursors. Analysis of the 611 putative genes revealed a striking feature: the dramatic expansion of several miRNA gene families. The largest family contained 91 genes that encoded 20 different miRNAs. Microarray analyses revealed the expression of miRNA genes was strictly regulated during Hessian fly larval development and abundance of many miRNA genes were affected by host genotypes.

Conclusion

The identification of a large number of miRNAs for the first time from a gall midge provides a foundation for further studies of miRNA functions in gall midge biology and behavior. The dramatic expansion of identical or similar miRNAs provides a unique system to study functional relations among miRNA iso-genes as well as changes in sequence specificity due to small changes in miRNAs and in their mRNA targets. These results may also facilitate the identification of miRNA genes for potential pest control through transgenic approaches.

Nasonia, a genus of four closely related parasitoid insect species, is a model system for genetic research. Their haplodiploid genetics (haploid males and diploid females) and interfertile species are advantageous for the genetic analysis of complex traits and the genetic basis of species differences. A fine-scale genomic map is an important tool for advancing genetic studies in this system. We developed and used a hybrid genotyping microarray to generate a high-resolution genetic map that covers 79% of the sequenced genome of Nasonia vitripennis. The microarray is based on differential hybridization of species-specific oligos between N. vitripennis and Nasonia giraulti at more than 20,000 markers spanning the Nasonia genome. The map places 729 scaffolds onto the five linkage groups of Nasonia, including locating many smaller scaffolds that would be difficult to map by other means. The microarray was used to characterize 26 segmental introgression lines containing chromosomal regions from one species in the genetic background of another. These segmental introgression lines have been used for rapid screening and mapping of quantitative trait loci involved in species differences. Finally, the microarray is extended to bulk-segregant analysis and genotyping of other Nasonia species combinations. These resources should further expand the usefulness of Nasonia for studies of the genetic basis and architecture of complex traits and speciation.

The evolutionary importance of hybridization and introgression has long been debated1. We used genomic tools to investigate introgression in Heliconius, a rapidly radiating genus of neotropical butterflies widely used in studies of ecology, behaviour, mimicry and speciation2-5 . We sequenced the genome of Heliconius melpomene and compared it with other taxa to investigate chromosomal evolution in Lepidoptera and gene flow among multiple Heliconius species and races. Among 12,657 predicted genes for Heliconius, biologically important expansions of families of chemosensory and Hox genes are particularly noteworthy. Chromosomal organisation has remained broadly conserved since the Cretaceous, when butterflies split from the silkmoth lineage. Using genomic resequencing, we show hybrid exchange of genes between three co-mimics, H. melpomene, H. timareta, and H. elevatus, especially at two genomic regions that control mimicry pattern. Closely related Heliconius species clearly exchange protective colour pattern genes promiscuously, implying a major role for hybridization in adaptive radiation.

Over the last ten years we have seen great efforts focused on revising amphibian systematics. Phylogenetic reconstructions derived from DNA sequence data have played a central role in these revisionary studies but have typically under-sampled the diverse frog family Microhylidae. Here, we present a detailed phylogenetic study focused on expanding previous hypotheses of relationships within this cosmopolitan family. Specifically, we placed an emphasis on assessing relationships among New World genera and those taxa with uncertain phylogenetic affinities (i.e., incertae sedis).

Results

One mitochondrial and three nuclear genes (about 2.8 kb) were sequenced to assess phylogenetic relationships. We utilized an unprecedented sampling of 200 microhylid taxa representing 91% of currently recognized subfamilies and 95% of New World genera. Our analyses do not fully resolve relationships among subfamilies supporting previous studies that have suggested a rapid early diversification of this clade. We observed a close relationship between Synapturanus and Otophryne of the subfamily Otophryninae. Within the subfamily Gastrophryninae relationships between genera were well resolved.

Conclusion

Otophryninae is distantly related to all other New World microhylids that were recovered as a monophyletic group, Gastrophryninae. Within Gastrophryninae, five genera were recovered as non-monophyletic; we propose taxonomic re-arrangements to render all genera monophyletic. This hypothesis of relationships and updated classification for New World microhylids may serve as a guide to better understand the evolutionary history of this group that is apparently subject to convergent morphological evolution and chromosome reduction. Based on a divergence analysis calibrated with hypotheses from previous studies and fossil data, it appears that microhylid genera inhabiting the New World originated during a period of gradual cooling from the late Oligocene to mid Miocene.

Many genomes have been sequenced to high-quality draft status using Sanger capillary electrophoresis and/or newer short-read sequence data and whole genome assembly techniques. However, even the best draft genomes contain gaps and other imperfections due to limitations in the input data and the techniques used to build draft assemblies. Sequencing biases, repetitive genomic features, genomic polymorphism, and other complicating factors all come together to make some regions difficult or impossible to assemble. Traditionally, draft genomes were upgraded to “phase 3 finished” status using time-consuming and expensive Sanger-based manual finishing processes. For more facile assembly and automated finishing of draft genomes, we present here an automated approach to finishing using long-reads from the Pacific Biosciences RS (PacBio) platform. Our algorithm and associated software tool, PBJelly, (publicly available at https://sourceforge.net/projects/pb-jelly/) automates the finishing process using long sequence reads in a reference-guided assembly process. PBJelly also provides “lift-over” co-ordinate tables to easily port existing annotations to the upgraded assembly. Using PBJelly and long PacBio reads, we upgraded the draft genome sequences of a simulated Drosophila melanogaster, the version 2 draft Drosophila pseudoobscura, an assembly of the Assemblathon 2.0 budgerigar dataset, and a preliminary assembly of the Sooty mangabey. With 24× mapped coverage of PacBio long-reads, we addressed 99% of gaps and were able to close 69% and improve 12% of all gaps in D. pseudoobscura. With 4× mapped coverage of PacBio long-reads we saw reads address 63% of gaps in our budgerigar assembly, of which 32% were closed and 63% improved. With 6.8× mapped coverage of mangabey PacBio long-reads we addressed 97% of gaps and closed 66% of addressed gaps and improved 19%. The accuracy of gap closure was validated by comparison to Sanger sequencing on gaps from the original D. pseudoobscura draft assembly and shown to be dependent on initial reference quality.

Predicting organismal phenotypes from genotype data is important for plant and animal breeding, medicine, and evolutionary biology. Genomic-based phenotype prediction has been applied for single-nucleotide polymorphism (SNP) genotyping platforms, but not using complete genome sequences. Here, we report genomic prediction for starvation stress resistance and startle response in Drosophila melanogaster, using ∼2.5 million SNPs determined by sequencing the Drosophila Genetic Reference Panel population of inbred lines. We constructed a genomic relationship matrix from the SNP data and used it in a genomic best linear unbiased prediction (GBLUP) model. We assessed predictive ability as the correlation between predicted genetic values and observed phenotypes by cross-validation, and found a predictive ability of 0.239±0.008 (0.230±0.012) for starvation resistance (startle response). The predictive ability of BayesB, a Bayesian method with internal SNP selection, was not greater than GBLUP. Selection of the 5% SNPs with either the highest absolute effect or variance explained did not improve predictive ability. Predictive ability decreased only when fewer than 150,000 SNPs were used to construct the genomic relationship matrix. We hypothesize that predictive power in this population stems from the SNP–based modeling of the subtle relationship structure caused by long-range linkage disequilibrium and not from population structure or SNPs in linkage disequilibrium with causal variants. We discuss the implications of these results for genomic prediction in other organisms.

Author Summary

The ability to accurately predict values of complex phenotypes from genotype data will revolutionize plant and animal breeding, personalized medicine, and evolutionary biology. To date, genomic prediction has utilized high-density single-nucleotide polymorphism (SNP) genotyping arrays, but the availability of sequence data opens new frontiers for genomic prediction methods. This article is the first application of genomic phenotype prediction using whole-genome sequence data in a substantial sample of a higher eukaryote. We use ∼2.5 million SNPs with minor allele frequency greater than 2.5% derived from genomic sequences of the “Drosophila Genetic Reference Panel” to predict phenotypes for two traits, starvation resistance and startle-induced locomotor behavior. We systematically address prediction within versus across sexes, genomic best linear unbiased prediction (GBLUP) versus a Bayesian approach, and the effect of SNP density. We find that (i) genomic prediction can be efficiently implemented using sequence data via GBLUP, (ii) there is little gain in predictive ability if the number of SNPs is increased above 150,000, and (iii) neither implicit nor explicit marker selection substantially improves the predictive ability. Although the findings must be seen against the background of small sample sizes, the results illustrate both the potential of the approach and the challenges ahead.

Genome sequencing technologies promise to revolutionize our understanding of genetics, evolution, and disease by making it feasible to survey a broad spectrum of sequence variation on a population scale. However, this potential can only be realized to the extent that methods for extracting and interpreting distinct forms of variation can be established. The error profiles and read length limitations of early versions of next-generation sequencing technologies rendered them ineffective for some sequence variant types, particularly microsatellites and other tandem repeats, and fostered the general misconception that such variants are inherently inaccessible to these platforms. At the same time, tandem repeats have emerged as important sources of functional variation. Tandem repeats are often located in and around genes, and frequent mutations in their lengths exert quantitative effects on gene function and phenotype, rapidly degrading linkage disequilibrium between markers and traits. Sensitive identification of these variants in large-scale next-gen sequencing efforts will enable more comprehensive association studies capable of revealing previously invisible associations. We present a population-scale analysis of microsatellite repeats using whole-genome data from 158 inbred isolates from the Drosophila Genetics Reference Panel, a collection of over 200 extensively phenotypically characterized isolates from a single natural population, to uncover processes underlying repeat mutation and to enable associations with behavioral, morphological, and life-history traits. Analysis of repeat variation from next-generation sequence data will also enhance studies of genome stability and neurodegenerative diseases.

Human meibomian gland dysfunction, a leading cause of dry eye, is accompanied by numerous changes in glandular gene expression. The nature of these alterations suggests that keratinization plays an important role in this disease.

Purpose.

Meibomian gland dysfunction (MGD) may be the leading cause of dry eye syndrome throughout the world. However, the precise mechanism(s) underlying the pathogenesis of this disease is unclear. This study was conducted to identify meibomian gland genes that may promote the development and/or progression of human MGD.

Methods.

Lid tissues were obtained from male and female MGD patients and age-matched controls after eyelid surgeries (e.g., to correct entropion or ectropion). Meibomian glands were isolated and processed for RNA extraction and the analysis of gene expression.

Results.

The results show that MGD is associated with significant alterations in the expression of almost 400 genes in the human meibomian gland. The levels of 197 transcripts, including those encoding various small proline-rich proteins and S100 calcium-binding proteins, are significantly increased, whereas the expression of 194 genes, such as claudin 3 and cell adhesion molecule 1, is significantly decreased. These changes, which cannot be accounted for by sex differences, are accompanied by alterations in many gene ontologies (e.g., keratinization, cell cycle, and DNA repair). The findings also show that the human meibomian gland contains several highly expressed genes that are distinct from those in an adjacent tissue (i.e., conjunctival epithelium).

Conclusions.

The results demonstrate that MGD is accompanied by multiple changes in gene expression in the meibomian gland. The nature of these alterations, including the upregulation of genes encoding small proline-rich proteins and S100 calcium-binding proteins, suggest that keratinization plays an important role in the pathogenesis of MGD.

Fieldwork conducted throughout Timor-Leste in September 2004 and July 2009 resulted in a collection or recording of 263 herpetological specimens (100 amphibians, 163 reptiles), comprising at least seven species of frogs and toads, 20 species of lizards, seven species of snakes, two species of turtles, and one species of crocodile. Among the amphibians, the most frequently encountered species were toads (Duttaphrynus melanostictus), rice paddy frogs (genus Fejervarya), and rhacophorid treefrogs (Polypedates cf. leucomystax). All three variants of rice paddy frogs encountered represent undescribed species similar to Fejervarya verruculosa from neighboring Wetar Island. Records of Fejervarya cancrivora and Fejervarya limnocharis for Timor Island are apparently errors based on misidentification. We obtained voucher specimens for a total of 147 lizards and voucher photographs only for four specimens of Varanus timorensis. Aside from geckos frequently associated with human habitations (e.g., Gehyra mutilata, Gekko gecko, Hemidactylus frenatus, Hemidactylus platyurus), we discovered an as yet undescribed species of bent-toed gecko, genus Cyrtodactylus, in the Same valley. Our specimens of Hemidactylus platyurus are the first record of this species from Timor-Leste. Commonly encountered skinks included four-fingered skinks (genus Carlia), wedge skinks (genus Sphenomorphus), and night skinks (genus Eremiascincus). Notable among the 15 snakes collected was the frequency of pitvipers (Cryptelytrops insularis), which amounted to over 25% of all snakes. Our specimen of the wolfsnake Lycodon subcinctus is the first record of this species for Timor-Leste. Based on these findings, it appears that the biodiversity of amphibians and reptiles in this remote corner of Wallacea is much greater than previously thought, particularly with respect to scincid lizards. The detail we provide in the species accounts is designed to allow the use of this report as a preliminary field guide to the amphibians and reptiles of Timor-Leste. However, survey work is ongoing.

Defects in sex steroid receptors have been linked to the onset, progression and severity, as well as the sex-related prevalence, of a variety of autoimmune disorders, including lupus, rheumatoid arthritis, multiple sclerosis and diabetes. We hypothesize that defects in estrogen receptor α (ESR1), estrogen receptor β (ESR2) and/or the androgen receptor (AR) may also contribute to the development of lacrimal gland autoimmune sequelae in Sjögren’s syndrome. To begin to test this hypothesis, we examined whether mutations exist in the coding regions of ESR1, ESR2 and AR transcripts in lacrimal tissues of mouse models of Sjögren’s syndrome.

Methods

Lacrimal and submandibular glands were collected from adult MRL/MpJ-Tnfrsf6lpr, nonobese diabetic and/or BALB/c mice. Tissues were pooled according to sex and experiment and processed for cDNA generation. PCR primers were designed to amplify 566–875 base pair segments of the entire open reading frame of each receptor. Segments were amplified, purified and then sequenced. Receptor sequences were assembled and compared to each other and to known NCBI sequences.

Results

Our results show that almost all ESR1, ESR2 and AR sequences in exocrine tissues of male and female autoimmune and non-autoimmune mice were identical to those of NCBI standards. There was a G→A shift at position 998 of the ESR2 complete coding sequence in all tissue samples when compared to NCBI reference sequence U81451.1, but this polymorphism was not found in other ESR2 reference sequences.

Conclusions

Our findings indicate that defects in the coding region of sex steroid receptors do not contribute to the pathogenesis of lacrimal gland disease in mouse models of Sjögren’s syndrome.

Paroxysmal nocturnal hemoglobinuria is a rare disorder of hemopoietic stem cells. Affected individuals have a triad of clinical associations – intravascular hemolysis, an increased risk of thromboembolism, and bone marrow failure. Most of the symptoms experienced in this disease occur due to the absence of complement regulatory proteins on the surface of the red blood cells. Complement activation is thus not checked and causes destruction of these cells. Eculizumab is a monoclonal antibody treatment which specifically binds to the complement protein C5, preventing its cleavage, and so halts the complement cascade and prevents the formation of the terminal complement proteins. Eculizumab prevents intravascular hemolysis, stabilizes hemoglobin levels, reduces or stops the need for blood transfusions, and improves fatigue and patient quality of life as well as reducing pulmonary hypertension, decreasing the risk of thrombosis and protecting against worsening renal function. It is not a curative therapy but has a great benefit on those with this rare debilitating condition.

The human X chromosome has a unique biology that was shaped by its evolution as the sex chromosome shared by males and females. We have determined 99.3% of the euchromatic sequence of the X chromosome. Our analysis illustrates the autosomal origin of the mammalian sex chromosomes, the stepwise process that led to the progressive loss of recombination between X and Y, and the extent of subsequent degradation of the Y chromosome. LINE1 repeat elements cover one-third of the X chromosome, with a distribution that is consistent with their proposed role as way stations in the process of X-chromosome inactivation. We found 1,098 genes in the sequence, of which 99 encode proteins expressed in testis and in various tumour types. A disproportionately high number of mendelian diseases are documented for the X chromosome. Of this number, 168 have been explained by mutations in 113 X-linked genes, which in many cases were characterized with the aid of the DNA sequence.

The whole genome sequence of Tribolium castaneum, a worldwide coleopteran pest of stored products, has recently been determined. In order to facilitate accurate annotation and detailed functional analysis of this genome, we have compiled and analyzed all available expressed sequence tag (EST) data. The raw data consist of 61,228 ESTs, including 10,704 obtained from NCBI and an additional 50,524 derived from 32,544 clones generated in our laboratories. These sequences were amassed from cDNA libraries representing six different tissues or stages, namely: whole embryos; whole larvae; larval hindguts and Malpighian tubules; larval fat bodies and carcasses; adult ovaries; and adult heads. Assembly of the 61,228 sequences collapsed into 12,269 clusters (groups of overlapping ESTs representing single genes), of which 10,134 mapped onto 6,463 (39%) of the 16,422 GLEAN gene models (i.e. official Tribolium gene list). Approximately 1,600 clusters (13% of the total) lack corresponding GLEAN models, despite high matches to the genome, suggesting that a considerable number of transcribed sequences were missed by the gene prediction programs or were removed by GLEAN. We conservatively estimate that the current EST set represents more than 7,500 transcription units.

We hypothesize that sex steroids induce sex-specific and/or opposite effects in the lacrimal and meibomian glands and that these actions may influence the prevalence of dry eye syndrome. The objective of this study was to begin to test this hypothesis.

Methods

Lacrimal and meibomian glands were obtained from ovariectomized mice that had been treated with testosterone or control vehicle for 14 days. Samples were processed for the isolation of RNA, and analyzed for differentially expressed mRNAs using CodeLink Bioarrays and quantitative real-time PCR (qPCR) techniques. Data were compared to those obtained following testosterone treatment of orchiectomized mice, as well as after the administration of 17β-estradiol and/or progesterone to ovariectomized mice.

Results

Our findings demonstrate that testosterone regulates the expression of thousands of genes in the lacrimal and meibomian glands of ovariectomized mice. The magnitude and extent of these hormonal effects, which encompassed numerous biological, molecular, and cellular ontologies, was tissue-dependent. Particularly notable was the androgen stimulation of meibomian gland genes related to lipid metabolic pathways, and the suppression of genes associated with keratinization. Many of the genes regulated by testosterone in female tissues were identical to those controlled by androgens in male lacrimal and meibomian glands. However, some genes were modulated in a sex-specific manner. In addition, a number of the androgen-regulated genes in female glands were altered in the opposite direction by 17β-estradiol and/or progesterone.

Conclusions

Our results support our hypothesis that sex steroids may induce sex-specific and/or opposite effects in the lacrimal and meibomian glands. Whether these actions contribute to the prevalence of dry eye remains to be determined.

Sex-associated differences have been identified in the anatomy, physiology and pathophysiology of the human cornea. We hypothesize that many of these differences are due to fundamental variations in gene expression. Our objective in this study was to determine whether such differences exist in human corneal epithelial cells both in vivo and in vitro.

Methods

Human corneal epithelial cells were isolated from the corneoscleral rims of male and female donors. Cells were processed either directly for RNA extraction, or first cultured in phenol red-free keratinocyte serum-free media. The RNA samples were examined for differentially expressed mRNAs by using of CodeLink Bioarrays and Affymetrix GeneChips. Data were analyzed with GeneSifter.Net software.

Results

Our results demonstrate that sex significantly influences the expression of over 600 genes in human corneal epithelial cells in vivo. These genes are involved in a broad spectrum of biologic processes, molecular functions and cellular components, such as metabolic processes, DNA replication, cell migration, RNA binding, oxidoreductase activity and nucleoli. We also identified significant, sex-related effects on gene expression in human corneal epithelial cells in vitro. However, with few exceptions (e.g., X- and Y-linked genes), these sex-related differences in gene expression in vitro were typically different than those in vivo.

Conclusions

Our findings support our hypothesis that sex-related differences exist in the gene expression of human corneal epithelial cells. Variations in gene expression may contribute to sex-related differences in the prevalence of certain corneal diseases.

As a first step to creating a comprehensive atlas of gene-expression patterns during Drosophila embryogenesis, 2,179 genes have been examinded by in situ hybridization to fixed Drosophila embryos. Of the genes assayed, 63.7% displayed dynamic expression patterns that were documented with 25,690 digital photomicrographs of individual embryos.

Background

Cell-fate specification and tissue differentiation during development are largely achieved by the regulation of gene transcription.

Results

As a first step to creating a comprehensive atlas of gene-expression patterns during Drosophila embryogenesis, we examined 2,179 genes by in situ hybridization to fixed Drosophila embryos. Of the genes assayed, 63.7% displayed dynamic expression patterns that were documented with 25,690 digital photomicrographs of individual embryos. The photomicrographs were annotated using controlled vocabularies for anatomical structures that are organized into a developmental hierarchy. We also generated a detailed time course of gene expression during embryogenesis using microarrays to provide an independent corroboration of the in situ hybridization results. All image, annotation and microarray data are stored in publicly available database. We found that the RNA transcripts of about 1% of genes show clear subcellular localization. Nearly all the annotated expression patterns are distinct. We present an approach for organizing the data by hierarchical clustering of annotation terms that allows us to group tissues that express similar sets of genes as well as genes displaying similar expression patterns.

Conclusions

Analyzing gene-expression patterns by in situ hybridization to whole-mount embryos provides an extremely rich dataset that can be used to identify genes involved in developmental processes that have been missed by traditional genetic analysis. Systematic analysis of rigorously annotated patterns of gene expression will complement and extend the types of analyses carried out using expression microarrays.