Ants (Hymenoptera, Formicidae) represent one of the most successful eusocial taxa in terms of both their geographic distribution and species number. The publication of seven ant genomes within the past year was a quantum leap for socio- and ant genomics. The diversity of social organization in ants makes them excellent model organisms to study the evolution of social systems. Comparing the ant genomes with those of the honeybee, a lineage that evolved eusociality independently from ants, and solitary insects suggests that there are significant differences in key aspects of genome organization between social and solitary insects, as well as among ant species. Altogether, these seven ant genomes open exciting new research avenues and opportunities for understanding the genetic basis and regulation of social species, and adaptive complex systems in general.

Background

DNA methylation is a common regulator of gene expression, including acting as a regulator of developmental events and behavioral changes in adults. Using the unique system of genetic caste determination in Pogonomyrmex barbatus, we were able to document changes in DNA methylation during development, and also across both ancient and contemporary hybridization events.

Methodology/Principal Findings

Sodium bisulfite sequencing demonstrated in vivo methylation of symmetric CG dinucleotides in P. barbatus. We also found methylation of non-CpG sequences. This validated two bioinformatics methods for predicting gene methylation, the bias in observed to expected ratio of CpG dinucleotides and the density of CpG/TpG single nucleotide polymorphisms (SNP). Frequencies of genomic DNA methylation were determined for different developmental stages and castes using ms-AFLP assays. The genetic caste determination system (GCD) is probably the product of an ancestral hybridization event between P. barbatus and P. rugosus. Two lineages obligately co-occur within a GCD population, and queens are derived from intra-lineage matings whereas workers are produced from inter-lineage matings. Relative DNA methylation levels of queens and workers from GCD lineages (contemporary hybrids) were not significantly different until adulthood. Virgin queens had significantly higher relative levels of DNA methylation compared to workers. Worker DNA methylation did not vary among developmental stages within each lineage, but was significantly different between the currently hybridizing lineages. Finally, workers of the two genetic caste determination lineages had half as many methylated cytosines as workers from the putative parental species, which have environmental caste determination.

Conclusions/Significance

These results suggest that DNA methylation may be a conserved regulatory mechanism moderating division of labor in both bees and ants. Current and historic hybridization appear to have altered genomic methylation levels suggesting a possible link between changes in overall DNA methylation and the origin and regulation of genetic caste determination in P. barbatus.

Vertebrates have experienced two rounds of whole-genome duplication (WGD) in the stem lineages of deep nodes within the group and a subsequent duplication event in the stem lineage of the teleosts—a highly diverse group of ray-finned fishes. Here, we present the first full Hox gene sequences for any member of the Acipenseriformes, the American paddlefish, and confirm that an independent WGD occurred in the paddlefish lineage, approximately 42 Ma based on sequences spanning the entire HoxA cluster and eight genes on the HoxD gene cluster. These clusters comprise different HOX loci and maintain conserved synteny relative to bichir, zebrafish, stickleback, and pufferfish, as well as human, mouse, and chick. We also provide a gene genealogy for the duplicated fzd8 gene in paddlefish and present evidence for the first Hox14 gene in any ray-finned fish. Taken together, these data demonstrate that the American paddlefish has an independently duplicated genome. Substitution patterns of the “alpha” paralogs on both the HoxA and HoxD gene clusters suggest transcriptional inactivation consistent with functional diploidization. Further, there are similarities in the pattern of sequence divergence among duplicated Hox genes in paddlefish and teleost lineages, even though they occurred independently approximately 200 Myr apart. We highlight implications on comparative analyses in the study of the “fin-limb transition” as well as gene and genome duplication in bony fishes, which includes all ray-finned fishes as well as the lobe-finned fishes and tetrapod vertebrates.

Honey bee colonies are subject to numerous pathogens and parasites. Interaction among multiple pathogens and parasites is the proposed cause for Colony Collapse Disorder (CCD), a syndrome characterized by worker bees abandoning their hive. Here we provide the first documentation that the phorid fly Apocephalus borealis, previously known to parasitize bumble bees, also infects and eventually kills honey bees and may pose an emerging threat to North American apiculture. Parasitized honey bees show hive abandonment behavior, leaving their hives at night and dying shortly thereafter. On average, seven days later up to 13 phorid larvae emerge from each dead bee and pupate away from the bee. Using DNA barcoding, we confirmed that phorids that emerged from honey bees and bumble bees were the same species. Microarray analyses of honey bees from infected hives revealed that these bees are often infected with deformed wing virus and Nosema ceranae. Larvae and adult phorids also tested positive for these pathogens, implicating the fly as a potential vector or reservoir of these honey bee pathogens. Phorid parasitism may affect hive viability since 77% of sites sampled in the San Francisco Bay Area were infected by the fly and microarray analyses detected phorids in commercial hives in South Dakota and California's Central Valley. Understanding details of phorid infection may shed light on similar hive abandonment behaviors seen in CCD.

Purpose:

Partial-thickness articular-sided rotator cuff tears have a multifactorial etiology and are associated with degeneration of the tendon. They are often described as an injury of the young athlete, although they are also found in the older population. The aim of this study was to investigate the frequency and associations of partial-thickness articular-sided tears in patients over the age of 35 years.

Design:

Retrospective

Materials and Methods:

A retrospective study of all arthroscopic procedures for rotator cuff pathology in patients over the age of 35 years over a 2-year period by a single surgeon was performed. The included patients were divided into two groups based on the arthroscopic findings: those with a partial-thickness articular-sided rotator cuff tear and those with pure tendinopathy. The groups were then compared to identify the associated pathology with the rotator cuff lesions. 2×2 contingency table analysis and unpaired Student's t-test were used for statistical analysis.

Results:

One hundred patients were included in the study of whom 62 had a partial articular-sided tear. Those with a partial articular-sided tear were older (P=0.0001), were more commonly associated with a documented injury (P=0.03), and more commonly had biceps degeneration (P=0.001) and synovitis (P=0.02) within the joint.

Conclusion:

Partial-thickness articular-sided tears are a common occurrence in patients requiring arthroscopic surgery for rotator cuff pathology over the age of 35 years. This probably reflects an injury in an already degenerate cuff. This would support the theory of intrinsic degeneration of the tendon in this age group and probably represent a different etiology to those seen in the young athletes.

Level of Evidence:

Level 3

Neurophysiological studies in monkeys show that activity of neurons in primary cortex (M1), pre-motor cortex (PMC), and cerebellum varies systematically with the direction of reaching movements. These neurons exhibit preferred direction tuning, where the level of neural activity is highest when movements are made in the preferred direction (PD), and gets progressively lower as movements are made at increasing degrees of offset from the PD. Using a functional magnetic resonance imaging adaptation (fMRI-A) paradigm, we show that PD coding does exist in regions of the human motor system that are homologous to those observed in non-human primates. Consistent with predictions of the PD model, we show adaptation (i.e., a lower level) of the blood oxygen level dependent (BOLD) time-course signal in M1, PMC, SMA, and cerebellum when consecutive wrist movements were made in the same direction (0° offset) relative to movements offset by 90° or 180°. The BOLD signal in dorsolateral prefrontal cortex adapted equally in all movement offset conditions, mitigating against the possibility that the present results are the consequence of differential task complexity or attention to action in each movement offset condition.

Background

Heterochromatin plays an important role in chromosome function and gene regulation. Despite the availability of polytene chromosomes and genome sequence, the heterochromatin of the major malaria vector Anopheles gambiae has not been mapped and characterized.

Results

To determine the extent of heterochromatin within the An. gambiae genome, genes were physically mapped to the euchromatin-heterochromatin transition zone of polytene chromosomes. The study found that a minimum of 232 genes reside in 16.6 Mb of mapped heterochromatin. Gene ontology analysis revealed that heterochromatin is enriched in genes with DNA-binding and regulatory activities. Immunostaining of the An. gambiae chromosomes with antibodies against Drosophila melanogaster heterochromatin protein 1 (HP1) and the nuclear envelope protein lamin Dm0 identified the major invariable sites of the proteins' localization in all regions of pericentric heterochromatin, diffuse intercalary heterochromatin, and euchromatic region 9C of the 2R arm, but not in the compact intercalary heterochromatin. To better understand the molecular differences among chromatin types, novel Bayesian statistical models were developed to analyze genome features. The study found that heterochromatin and euchromatin differ in gene density and the coverage of retroelements and segmental duplications. The pericentric heterochromatin had the highest coverage of retroelements and tandem repeats, while intercalary heterochromatin was enriched with segmental duplications. We also provide evidence that the diffuse intercalary heterochromatin has a higher coverage of DNA transposable elements, minisatellites, and satellites than does the compact intercalary heterochromatin. The investigation of 42-Mb assembly of unmapped genomic scaffolds showed that it has molecular characteristics similar to cytologically mapped heterochromatin.

Conclusions

Our results demonstrate that Anopheles polytene chromosomes and whole-genome shotgun assembly render the mapping and characterization of a significant part of heterochromatic scaffolds a possibility. These results reveal the strong association between characteristics of the genome features and morphological types of chromatin. Initial analysis of the An. gambiae heterochromatin provides a framework for its functional characterization and comparative genomic analyses with other organisms.

Background

Nonrandom distribution of rearrangements is a common feature of eukaryotic chromosomes that is not well understood in terms of genome organization and evolution. In the major African malaria vector Anopheles gambiae, polymorphic inversions are highly nonuniformly distributed among five chromosomal arms and are associated with epidemiologically important adaptations. However, it is not clear whether the genomic content of the chromosomal arms is associated with inversion polymorphism and fixation rates.

Methodology/Principal Findings

To better understand the evolutionary dynamics of chromosomal inversions, we created a physical map for an Asian malaria mosquito, Anopheles stephensi, and compared it with the genome of An. gambiae. We also developed and deployed novel Bayesian statistical models to analyze genome landscapes in individual chromosomal arms An. gambiae. Here, we demonstrate that, despite the paucity of inversion polymorphisms on the X chromosome, this chromosome has the fastest rate of inversion fixation and the highest density of transposable elements, simple DNA repeats, and GC content. The highly polymorphic and rapidly evolving autosomal 2R arm had overrepresentation of genes involved in cellular response to stress supporting the role of natural selection in maintaining adaptive polymorphic inversions. In addition, the 2R arm had the highest density of regions involved in segmental duplications that clustered in the breakpoint-rich zone of the arm. In contrast, the slower evolving 2L, 3R, and 3L, arms were enriched with matrix-attachment regions that potentially contribute to chromosome stability in the cell nucleus.

Conclusions/Significance

These results highlight fundamental differences in evolutionary dynamics of the sex chromosome and autosomes and revealed the strong association between characteristics of the genome landscape and rates of chromosomal evolution. We conclude that a unique combination of various classes of genes and repetitive DNA in each arm, rather than a single type of repetitive element, is likely responsible for arm-specific rates of rearrangements.

The repetitive DNA that constitutes most of the heterochromatic regions of metazoan genomes has hindered the comprehensive analysis of gene content and other functions. We have generated a detailed computational and manual annotation of 24 megabases of heterochromatic sequence in the Release 5 Drosophila melanogaster genome sequence. The heterochromatin contains a minimum of 230 to 254 protein-coding genes, which are conserved in other Drosophilids and more diverged species, as well as 32 pseudogenes and 13 noncoding RNAs. Improved methods revealed that more than 77% of this heterochromatin sequence, including introns and intergenic regions, is composed of fragmented and nested transposable elements and other repeated DNAs. Drosophila heterochromatin contains “islands” of highly conserved genes embedded in these “oceans” of complex repeats, which may require special expression and splicing mechanisms.

Annotation of an improved whole-genome shotgun assembly of the Drosophila melanogaster genome predicted 297 protein-coding genes and six non-protein-coding genes, including known heterochromatic genes, and regions of similarity to known transposable elements. Fluorescence in situ hybridization was used to correlate the genomic sequence with the cytogenetic map; the annotated euchromatic sequence extends into the centric heterochromatin on each chromosome arm.

Background

Most eukaryotic genomes include a substantial repeat-rich fraction termed heterochromatin, which is concentrated in centric and telomeric regions. The repetitive nature of heterochromatic sequence makes it difficult to assemble and analyze. To better understand the heterochromatic component of the Drosophila melanogaster genome, we characterized and annotated portions of a whole-genome shotgun sequence assembly.

Results

WGS3, an improved whole-genome shotgun assembly, includes 20.7 Mb of draft-quality sequence not represented in the Release 3 sequence spanning the euchromatin. We annotated this sequence using the methods employed in the re-annotation of the Release 3 euchromatic sequence. This analysis predicted 297 protein-coding genes and six non-protein-coding genes, including known heterochromatic genes, and regions of similarity to known transposable elements. Bacterial artificial chromosome (BAC)-based fluorescence in situ hybridization analysis was used to correlate the genomic sequence with the cytogenetic map in order to refine the genomic definition of the centric heterochromatin; on the basis of our cytological definition, the annotated Release 3 euchromatic sequence extends into the centric heterochromatin on each chromosome arm.

Conclusions

Whole-genome shotgun assembly produced a reliable draft-quality sequence of a significant part of the Drosophila heterochromatin. Annotation of this sequence defined the intron-exon structures of 30 known protein-coding genes and 267 protein-coding gene models. The cytogenetic mapping suggests that an additional 150 predicted genes are located in heterochromatin at the base of the Release 3 euchromatic sequence. Our analysis suggests strategies for improving the sequence and annotation of the heterochromatic portions of the Drosophila and other complex genomes.

The recent completion of the Drosophila melanogaster genomic sequence to high quality, and the availability of a greatly expanded set of Drosophila cDNA sequences, afforded FlyBase the opportunity to significantly improve genomic annotations.

Background

The recent completion of the Drosophila melanogaster genomic sequence to high quality and the availability of a greatly expanded set of Drosophila cDNA sequences, aligning to 78% of the predicted euchromatic genes, afforded FlyBase the opportunity to significantly improve genomic annotations. We made the annotation process more rigorous by inspecting each gene visually, utilizing a comprehensive set of curation rules, requiring traceable evidence for each gene model, and comparing each predicted peptide to SWISS-PROT and TrEMBL sequences.

Results

Although the number of predicted protein-coding genes in Drosophila remains essentially unchanged, the revised annotation significantly improves gene models, resulting in structural changes to 85% of the transcripts and 45% of the predicted proteins. We annotated transposable elements and non-protein-coding RNAs as new features, and extended the annotation of untranslated (UTR) sequences and alternative transcripts to include more than 70% and 20% of genes, respectively. Finally, cDNA sequence provided evidence for dicistronic transcripts, neighboring genes with overlapping UTRs on the same DNA sequence strand, alternatively spliced genes that encode distinct, non-overlapping peptides, and numerous nested genes.

Conclusions

Identification of so many unusual gene models not only suggests that some mechanisms for gene regulation are more prevalent than previously believed, but also underscores the complex challenges of eukaryotic gene prediction. At present, experimental data and human curation remain essential to generate high-quality genome annotations.

Telomeres are the protein–nucleic acid structures at the ends of eukaryote chromosomes. Tandem repeats of telomeric DNA are templated by the RNA component (TER1) of the ribonucleoprotein telomerase. These repeats are bound by telomere binding proteins, which are thought to interact with other factors to create a higher-order cap complex that stabilizes the chromosome end. In the budding yeast Kluyveromyces lactis, the incorporation of certain mutant DNA sequences into telomeres leads to uncapping of telomeres, manifested by dramatic telomere elongation and increased length heterogeneity (telomere deregulation). Here we show that telomere deregulation leads to enlarged, misshapen “monster” cells with increased DNA content and apparent defects in cell division. However, such deregulated telomeres became stabilized at their elongated lengths upon addition of only a few functionally wild-type telomeric repeats to their ends, after which the frequency of monster cells decreased to wild-type levels. These results provide evidence for the importance of the most terminal repeats at the telomere in maintaining the cap complex essential for normal telomere function. Analysis of uncapped and capped telomeres also show that it is the deregulation resulting from telomere uncapping, rather than excessive telomere length per se, that is associated with DNA aberrations and morphological defects.