BCOR Internal Tandem Duplication Associated Uterine Sarcoma: Expanding the Clinicopathologic Spectrum.

The diagnosis of high-grade endometrial stromal sarcoma has become more refined following molecular characterization of these tumors. Recently BCOR internal tandem duplications (ITD) have been identified in a small number of high-grade endometrial stromal sarcoma. Here we present an additional case of this rare entity in a young woman in her late teens. She presented with menorrhagia and underwent resection of 2 uterine lesions. The tumor was a spindle cell neoplasm composed of long fascicles with low to moderate cellularity, mild to moderate cytologic atypia, and up to 2 mitotic figures per 10 high power fields. Necrosis was not identified. Immunohistochemical stains showed the tumor to be positive for cyclin D1 in >50% of tumor cells, focally positive for CD10, and negative for SMA, desmin, h-caldesmon, and ALK1. BCOR ITD was confirmed by polymerase chain reaction with subsequent Sanger sequencing. Clues to the diagnosis of BCOR ITD uterine sarcoma include young patient age, uniform nuclear features, and diffuse positivity for cyclin D1. These features should prompt further molecular interrogation for definitive diagnosis, which is important for prognostication.

Assessment of the Implementation of Pharmacogenomic Testing in a Pediatric Tertiary Care Setting.

Importance: Pharmacogenomic (PGx) testing provides preemptive pharmacotherapeutic guidance regarding the lack of therapeutic benefit or adverse drug reactions of PGx targeted drugs. Pharmacogenomic information is of particular value among children with complex medical conditions who receive multiple medications and are at higher risk of developing adverse drug reactions. Objectives: To assess the implementation outcomes of a PGx testing program comprising both a point-of-care model that examined targeted drugs and a preemptive model informed by whole-genome sequencing that evaluated a broad range of drugs for potential therapy among children in a pediatric tertiary care setting. Design, Setting, and Participants: This cohort study was conducted at The Hospital for Sick Children in Toronto, Ontario, from January 2017 to September 2020. Pharmacogenomic analyses were performed among 172 children who were categorized into 2 groups: a point-of-care cohort and a preemptive cohort. The point-of-care cohort comprised 57 patients referred to the consultation clinic for planned therapy with PGx targeted drugs and/or for adverse drug reactions, including lack of therapeutic benefit, after the receipt of current or past medications. The preemptive cohort comprised 115 patients who received exploratory whole-genome sequencing-guided PGx testing for their heart conditions from the cardiac genome clinic at the Ted Rogers Centre for Heart Research. Exposures: Patients received PGx analysis of whole-genome sequencing data and/or multiplex genotyping of 6 pharmacogenes (CYP2C19, CYP2C9, CYP2D6, CYP3A5, VKORC1, and TPMT) that have established PGx clinical guidelines. Main Outcomes and Measures: The number of patients for whom PGx test results warranted deviation from standard dosing regimens. Results: A total of 172 children (mean [SD] age, 8.5 [5.6] years; 108 boys [62.8%]) were enrolled in the study. In the point-of-care cohort, a median of 2 target genes (range, 1-5 genes) were investigated per individual, with CYP2C19 being the most frequently examined; genotypes in 21 of 57 children (36.8%) were incompatible with standard treatment regimens. As expected from population allelic frequencies, among the 115 children in the whole-genome sequencing-guided preemptive cohort, 92 children (80.0%) were recommended to receive nonstandard treatment regimens for potential drug therapies based on their 6-gene pharmacogenetic profile. Conclusions and Relevance: In this cohort study, among both the point-of-care and preemptive cohorts, the multiplex PGx testing program provided dosing recommendations that deviated from standard regimens at an overall rate that was similar to the population frequencies of relevant variants.

Structural neuroimaging correlates of social deficits are similar in autism spectrum disorder and attention-deficit/hyperactivity disorder: analysis from the POND Network.

Autism spectrum disorder (ASD), attention-deficit/hyperactivity disorder (ADHD), and obsessive-compulsive disorder (OCD) have been associated with difficulties recognizing and responding to social cues. Neuroimaging studies have begun to map the social brain; however, the specific neural substrates contributing to social deficits in neurodevelopmental disorders remain unclear. Three hundred and twelve children underwent structural magnetic resonance imaging of the brain (controls = 32, OCD = 44, ADHD = 77, ASD = 159; mean age = 11). Their social deficits were quantified on the Social Communication Questionnaire (SCQ) and the Reading the Mind in the Eyes Test (RMET). Multivariable regression models were used to examine the structural neuroimaging correlates of social deficits, with both a region of interest and a whole-brain vertex-wise approach. For the region of interest analysis, social brain regions were grouped into three networks: (1) lateral mentalization (e.g., temporal-parietal junction), (2) frontal cognitive (e.g., orbitofrontal cortex), and (3) subcortical affective (e.g., limbic system) regions. Overall, social communication deficits on the SCQ were associated with thinner cortices in the left lateral regions and the right insula, and decreased volume in the ventral striatum, across diagnostic groups (p = 0.006 to <0.0001). Smaller subcortical volumes were associated with more severe social deficits on the SCQ in ASD and ADHD, and less severe deficits in OCD. On the RMET, larger amygdala/hippocampal volumes were associated with fewer deficits across groups. Overall, patterns of associations were similar in ASD and ADHD, supporting a common underlying biology and the blurring of the diagnostic boundaries between these disorders.

Epstein-Barr virus latent gene EBNA-1 genetic diversity among transplant patients compared with patients with infectious mononucleosis.

INTRODUCTION: As a step toward evaluating the association between Epstein-Barr virus genetic diversity and post-transplant lymphoproliferative disorder (PTLD), we conducted a preliminary study to compare the genetic diversity of the EBNA-1 gene among transplant patients and patients with infectious mononucleosis (IM). METHODS: We sequenced the EBNA-1 gene in blood samples from study subjects using Sanger methodology. The sequences were aligned with a reference strain and compared with publicly available sequences. RESULTS: We analyzed 33 study samples and 25 publicly available sequences along with the reference strain B95-8. The evaluable samples were from sixteen patients with IM (median age 14.0 years, range 2-24) and 17 transplant patients. There were six children without PTLD (median age 1.93 years, range 0.79-7.46) and 11 who developed PTLD (median age 5.67 years, range 0.96-17.45). A predominant EBNA-1 variant (P-thr) was identified across the study groups. Differences were observed between the samples from the IM patients compared with the transplant samples. CONCLUSION: The predominant EBNA-1 strain is in contrast to reports of the predominant strain in North America. The results suggest differences between the EBNA-1 strains among the study groups. Further studies will examine the relationship between EBNA-1 strains and PTLD occurrence and outcomes.

The Personal Genome Project Canada: findings from whole genome sequences of the inaugural 56 participants.

BACKGROUND: The Personal Genome Project Canada is a comprehensive public data resource that integrates whole genome sequencing data and health information. We describe genomic variation identified in the initial recruitment cohort of 56 volunteers. METHODS: Volunteers were screened for eligibility and provided informed consent for open data sharing. Using blood DNA, we performed whole genome sequencing and identified all possible classes of DNA variants. A genetic counsellor explained the implication of the results to each participant. RESULTS: Whole genome sequencing of the first 56 participants identified 207 662 805 sequence variants and 27 494 copy number variations. We analyzed a prioritized disease-associated data set (n = 1606 variants) according to standardized guidelines, and interpreted 19 variants in 14 participants (25%) as having obvious health implications. Six of these variants (e.g., in BRCA1 or mosaic loss of an X chromosome) were pathogenic or likely pathogenic. Seven were risk factors for cancer, cardiovascular or neurobehavioural conditions. Four other variants - associated with cancer, cardiac or neurodegenerative phenotypes - remained of uncertain significance because of discrepancies among databases. We also identified a large structural chromosome aberration and a likely pathogenic mitochondrial variant. There were 172 recessive disease alleles (e.g., 5 individuals carried mutations for cystic fibrosis). Pharmacogenomics analyses revealed another 3.9 potentially relevant genotypes per individual. INTERPRETATION: Our analyses identified a spectrum of genetic variants with potential health impact in 25% of participants. When also considering recessive alleles and variants with potential pharmacologic relevance, all 56 participants had medically relevant findings. Although access is mostly limited to research, whole genome sequencing can provide specific and novel information with the potential of major impact for health care.

Genome sequencing as a platform for pharmacogenetic genotyping: a pediatric cohort study.

Whole-genome sequencing and whole-exome sequencing have proven valuable for diagnosing inherited diseases, particularly in children. However, usage of sequencing data as a pharmacogenetic screening tool to ensure medication safety and effectiveness remains to be explored. Sixty-seven variants in 19 genes with known effects on drug response were compared between genome sequencing and targeted genotyping data for coverage and concordance in 98 pediatric patients. We used targeted genotyping data as a benchmark to assess accuracy of variant calling, and to identify copy number variations of the CYP2D6 gene. We then predicted clinical impact of these variants on drug therapy. We find genotype concordance across those panels to be > 97%. Concordance of CYP2D6 predicted phenotype between estimates of whole-genome sequencing and targeted genotyping panel were 90%; a result from a lower coverage depth or variant calling difficulties in our whole-genome sequencing data when copy number variation and/or the CYP2D6*4 haplotype were present. Importantly, 95 children had at least one clinically actionable pharmacogenetic variant. Diagnostic genomic sequencing data can be used for pre-emptive pharmacogenetic screening. However, concordance between genome-wide sequencing and target genotyping needs to be characterized for each of the pharmacologically important genes.

Oxytocin Receptor Polymorphisms are Differentially Associated with Social Abilities across Neurodevelopmental Disorders.

Oxytocin is a pituitary neuropeptide that affects social behaviour. Single nucleotide polymorphisms (SNPs) in the oxytocin receptor gene (OXTR) have been shown to explain some variability in social abilities in control populations. Whether these variants similarly contribute to the severity of social deficits experienced by children with neurodevelopmental disorders is unclear. Social abilities were assessed in a group of children with autism spectrum disorder (ASD, n = 341) or attention deficit hyperactivity disorder (ADHD, n = 276) using two established social measures. Scores were compared by OXTR genotype (rs53576, rs237887, rs13316193, rs2254298). Unexpectedly, the two most frequently studied OXTR SNPs in the general population (rs53576 and rs2254298) were associated with an increased severity of social deficits in ASD (p < 0.0001 and p = 0.0005), yet fewer social deficits in ADHD (p = 0.007 and p < 0.0001). We conclude that these genetic modifier alleles are not inherently risk-conferring with respect to their impact on social abilities; molecular investigations are greatly needed.

Novel 25 kb Deletion of MERTK Causes Retinitis Pigmentosa With Severe Progression.

Purpose: Retinitis pigmentosa (RP) describes a complex group of inherited retinal dystrophies with almost 300 reported genes and loci. We investigated the genetic etiology of autosomal recessive RP (arRP) in a large kindred with 5 affected family members, who reside on the island of Newfoundland, Canada. Methods: Genetic linkage analysis was performed on 12 family members (Infinium HumanOmni2.5-8 BeadChip). Whole exome sequencing analysis (Illumina HiSeq) was performed on one affected individual. A custom pipeline was applied to call, annotate, and filter variants. FishingCNV was used to scan the exome for rare copy number variants (CNVs). Candidate CNVs subsequently were visualized from microarray data (CNVPartition v.3.1.6.). MERTK breakpoints were mapped and familial cosegregation was tested using Sanger Sequencing. Results: We found strong evidence of linkage to a locus on chromosome 2 (logarithm of the odds [LOD] 4.89 [θ = 0]), at an interval encompassing the MERTK gene. Whole exome sequencing did not uncover candidate point mutations in MERTK, or other known RP genes. Subsequently, CNV analysis of the exome data and breakpoint mapping revealed a 25,218 bp deletion of MERTK, encompassing exons 6 to 8, with breakpoints in introns 5 (chr2:112,725,292) and 8 (chr2:112,750,421). A 48 bp insertion sequence was buried within the breakpoint; 18 bps shared homology to MIR4435-2HG and LINC00152, and 30 bp mapped to MERTK. The deletion cosegregated with arRP in the family. Conclusions: This study describes the molecular and clinical characterization of an arRP family segregating a novel 25 kb deletion of MERTK. These findings may assist clinicians in providing a diagnosis for other unsolved RP cases.

De Novo Genome and Transcriptome Assembly of the Canadian Beaver (Castor canadensis).

The Canadian beaver (Castor canadensis) is the largest indigenous rodent in North America. We report a draft annotated assembly of the beaver genome, the first for a large rodent and the first mammalian genome assembled directly from uncorrected and moderate coverage (< 30 ×) long reads generated by single-molecule sequencing. The genome size is 2.7 Gb estimated by k-mer analysis. We assembled the beaver genome using the new Canu assembler optimized for noisy reads. The resulting assembly was refined using Pilon supported by short reads (80 ×) and checked for accuracy by congruency against an independent short read assembly. We scaffolded the assembly using the exon-gene models derived from 9805 full-length open reading frames (FL-ORFs) constructed from the beaver leukocyte and muscle transcriptomes. The final assembly comprised 22,515 contigs with an N50 of 278,680 bp and an N50-scaffold of 317,558 bp. Maximum contig and scaffold lengths were 3.3 and 4.2 Mb, respectively, with a combined scaffold length representing 92% of the estimated genome size. The completeness and accuracy of the scaffold assembly was demonstrated by the precise exon placement for 91.1% of the 9805 assembled FL-ORFs and 83.1% of the BUSCO (Benchmarking Universal Single-Copy Orthologs) gene set used to assess the quality of genome assemblies. Well-represented were genes involved in dentition and enamel deposition, defining characteristics of rodents with which the beaver is well-endowed. The study provides insights for genome assembly and an important genomics resource for Castoridae and rodent evolutionary biology.

The genetic diversity of Epstein-Barr virus in the setting of transplantation relative to non-transplant settings: A feasibility study.

This study examines EBV strains from transplant patients and patients with IM by sequencing major EBV genes. We also used NGS to detect EBV DNA within total genomic DNA, and to evaluate its genetic variation. Sanger sequencing of major EBV genes was used to compare SNVs from samples taken from transplant patients vs. patients with IM. We sequenced EBV DNA from a healthy EBV-seropositive individual on a HiSeq 2000 instrument. Data were mapped to the EBV reference genomes (AG876 and B95-8). The number of EBNA2 SNVs was higher than for EBNA1 and the other genes sequenced within comparable reference coordinates. For EBNA2, there was a median of 15 SNV among transplant samples compared with 10 among IM samples (p = 0.036). EBNA1 showed little variation between samples. For NGS, we identified 640 and 892 variants at an unadjusted p value of 5 × 10(-8) for AG876 and B95-8 genomes, respectively. We used complementary sequence strategies to examine EBV genetic diversity and its application to transplantation. The results provide the framework for further characterization of EBV strains and related outcomes after organ transplantation.

GRIN1 polymorphisms do not affect susceptibility or phenotype in NMDA receptor encephalitis.

OBJECTIVE: To determine whether distinct single nucleotide polymorphisms (SNPs) within the glutamate receptor ionotropic NMDA 1 gene (GRIN1) are associated with NMDA receptor (NMDAR) encephalitis and whether these same variants are associated with variability in the clinical presentation and course of affected patients. METHODS: We performed clinical follow-up on 48 patients with NMDAR encephalitis and NMDAR autoantibodies detected in serum or CSF. All RefSeq GRIN1 coding exons were sequenced in 39 Caucasian-European patients, and the frequencies of SNPs were compared with those of an ethnically similar population using a case-control study design. Predetermined clinical variables were compared between patients with and without identified SNPs. RESULTS: Two SNPs were identified in GRIN1: 24 (62%) Caucasian-European patients with NMDAR encephalitis had alternate alleles at both rs6293 (exon 6) and rs1126442 (exon 7; exon numbering according to NM_001185090). The SNPs were in complete linkage disequilibrium. The frequency of these variants did not differ between patients with NMDAR encephalitis and ethnically matched individuals in the general population. No differences in clinical presentation, measures of disease severity, clinical course, or outcomes were observed between patients with different genotypes at these SNPs. CONCLUSION: Disease susceptibility or course in patients with NMDAR encephalitis was not strongly affected by SNPs in GRIN1. This study provides an estimate of the frequency of SNPs in GRIN1 in patients with NMDAR encephalitis and emphasizes the need for multisite collaborative studies enrolling larger numbers of patients to identify the genetic contributions to NMDAR encephalitis.

Whole-Exome Sequencing and Targeted Copy Number Analysis in Primary Ciliary Dyskinesia.

Primary ciliary dyskinesia (PCD) is an autosomal-recessive disorder resulting from loss of normal ciliary function. Symptoms include neonatal respiratory distress, chronic sinusitis, bronchiectasis, situs inversus, and infertility. Clinical features may be subtle and highly variable, making the diagnosis of PCD challenging. The diagnosis can be confirmed with ciliary ultrastructure analysis and/or molecular genetic testing of 32 PCD-associated genes. However, because of this genetic heterogeneity, comprehensive molecular genetic testing is not considered the standard of care, and the most efficient molecular approach has yet to be elucidated. Here, we propose a cost-effective and time-efficient molecular genetic algorithm to solve cases of PCD. We conducted targeted copy number variation (CNV) analysis and/or whole-exome sequencing on 20 families (22 patients) from a subset of 45 families (52 patients) with a clinical diagnosis of PCD who did not have a molecular genetic diagnosis after Sanger sequencing of 12 PCD-associated genes. This combined molecular genetic approach led to the identification of 4 of 20 (20%) families with clinically significant CNVs and 7 of 20 (35%) families with biallelic pathogenic mutations in recently identified PCD genes, resulting in an increased molecular genetic diagnostic rate of 55% (11/20). In patients with a clinical diagnosis of PCD, whole-exome sequencing followed by targeted CNV analysis results in an overall molecular genetic yield of 76% (34/45).

Long read nanopore sequencing for detection of HLA and CYP2D6 variants and haplotypes.

Haplotypes are often critical for the interpretation of genetic laboratory observations into medically actionable findings. Current massively parallel DNA sequencing technologies produce short sequence reads that are often unable to resolve haplotype information. Phasing short read data typically requires supplemental statistical phasing based on known haplotype structure in the population or parental genotypic data. Here we demonstrate that the MinION nanopore sequencer is capable of producing very long reads to resolve both variants and haplotypes of HLA-A, HLA-B and CYP2D6 genes important in determining patient drug response in sample NA12878 of CEPH/UTAH pedigree 1463, without the need for statistical phasing. Long read data from a single 24-hour nanopore sequencing run was used to reconstruct haplotypes, which were confirmed by HapMap data and statistically phased Complete Genomics and Sequenom genotypes. Our results demonstrate that nanopore sequencing is an emerging standalone technology with potential utility in a clinical environment to aid in medical decision-making.

Discrimination of SNPs in GC-rich regions using a modified hydrolysis probe chemistry protocol.

Allelic discrimination using TaqMan 5'-nuclease assay chemistry has been in routine use for many years, and the catalog of Life Technologies' predesigned SNP genotyping assays now exceeds 4 million entries. However, predesigned assays are often not available for genomic regions with a high GC content, nor can an assay necessarily be designed in this type of region using the manufacturer's design pipelines. Additionally, when an assay is available, the performance can be poor when using standard protocols. Here we report a modified allelic discrimination protocol for variants that reside in extremely GC-rich (GC > 75%) regions. The approach resolves fluorescent signal from reference and variant alleles, allowing all samples to be successfully assigned a genotype call. This protocol modification adds an extra step to the standard workflow, but the increased time is a productive compromise to generate high-quality data.

OTX2 mutations cause autosomal dominant pattern dystrophy of the retinal pigment epithelium.

PURPOSE: To identify the genetic cause of autosomal-dominant pattern dystrophy (PD) of the retinal pigment epithelium (RPE) in two families. METHODS AND RESULTS: Two families with autosomal-dominant PD were identified. Eight members of family 1 (five affected) were subjected to whole-genome SNP genotyping; multipoint genome-wide linkage analysis identified 7 regions of potential linkage, and genotyping four additional individuals from family 1 resulted in a maximum logarithm of odds score of 2.09 observed across four chromosomal regions. Exome sequencing of two affected family 1 members identified 15 shared non-synonymous rare coding sequence variants within the linked regions; candidate genes were prioritised and further analysed. Sanger sequencing confirmed a novel heterozygous missense variant (E79K) in orthodenticle homeobox 2 (OTX2) that segregated with the disease phenotype. Family 2 with PD (two affected) harboured the same missense variant in OTX2. A shared haplotype of 19.68 cM encompassing OTX2 was identified between affected individuals in the two families. Within the two families, all except one affected demonstrated distinct 'patterns' at the macula. In vivo structural retinal imaging showed discrete areas of RPE-photoreceptor separation at the macula in all cases. Electroretinogram testing showed generalised photoreceptor degeneration in three cases. Mild developmental anomalies were observed, including optic nerve head dysplasia (four cases), microcornea (one case) and Rathke's cleft cyst (one case); pituitary hormone levels were normal. CONCLUSIONS: This is the first report implicating OTX2 to underlie PD. The retinal disease resembles conditional mice models that show slow photoreceptor degeneration secondary to loss of Otx2 function in the adult RPE.

Genome-wide investigation of DNA methylation marks associated with FV Leiden mutation.

In order to investigate whether DNA methylation marks could contribute to the incomplete penetrance of the FV Leiden mutation, a major genetic risk factor for venous thrombosis (VT), we measured genome-wide DNA methylation levels in peripheral blood samples of 98 VT patients carrying the mutation and 251 VT patients without the mutation using the dedicated Illumina HumanMethylation450 array. The genome-wide analysis of 388,120 CpG probes identified three sites mapping to the SLC19A2 locus whose DNA methylation levels differed significantly (p<3 10-8) between carriers and non-carriers. The three sites replicated (p<2 10-7) in an independent sample of 214 individuals from five large families ascertained on VT and FV Leiden mutation among which 53 were carriers and 161 were non-carriers of the mutation. In both studies, these three CpG sites were also associated (2.33 10-110.05). In conclusion, our work clearly illustrates some promises and pitfalls of DNA methylation investigations on peripheral blood DNA in large epidemiological cohorts. DNA methylation levels at SLC19A2 are influenced by SNPs in LD with FV Leiden, but these DNA methylation marks do not explain the incomplete penetrance of the FV Leiden mutation.

Reprever: resolving low-copy duplicated sequences using template driven assembly.

Genomic sequence duplication is an important mechanism for genome evolution, often resulting in large sequence variations with implications for disease progression. Although paired-end sequencing technologies are commonly used for structural variation discovery, the discovery of novel duplicated sequences remains an unmet challenge. We analyze duplicons starting from identified high-copy number variants. Given paired-end mapped reads, and a candidate high-copy region, our tool, Reprever, identifies (a) the insertion breakpoints where the extra duplicons inserted into the donor genome and (b) the actual sequence of the duplicon. Reprever resolves ambiguous mapping signatures from existing homologs, repetitive elements and sequencing errors to identify breakpoint. At each breakpoint, Reprever reconstructs the inserted sequence using profile hidden Markov model (PHMM)-based guided assembly. In a test on 1000 artificial genomes with simulated duplication, Reprever could identify novel duplicates up to 97% of genomes within 3 bp positional and 1% sequence errors. Validation on 680 fosmid sequences identified and reconstructed eight duplicated sequences with high accuracy. We applied Reprever to reanalyzing a re-sequenced data set from the African individual NA18507 to identify >800 novel duplicates, including insertions in genes and insertions with additional variation. polymerase chain reaction followed by capillary sequencing validated both the insertion locations of the strongest predictions and their predicted sequence.

Phylogenetic relationships and divergence times of Charadriiformes genera: multigene evidence for the Cretaceous origin of at least 14 clades of shorebirds.

Comparative study of character evolution in the shorebirds is presently limited because the phylogenetic placement of some enigmatic genera remains unclear. We therefore used Bayesian methods to obtain a well-supported phylogeny of 90 recognized genera using 5 kb of mitochondrial and nuclear sequences. The tree comprised three major clades: Lari (gulls, auks and allies plus buttonquails) as sister to Scolopaci (sandpipers, jacanas and allies), and in turn sister to Charadrii (plovers, oystercatchers and allies), as in previous molecular studies. Plovers and noddies were not recovered as monophyletic assemblages, and the Egyptian plover Pluvianus is apparently not a plover. Molecular dating using multiple fossil constraints suggests that the three suborders originated in the late Cretaceous between 79 and 102 Mya, and at least 14 lineages of modern shorebirds survived the mass extinction at the K/T boundary. Previous difficulties in determining the phylogenetic relationships of enigmatic taxa reflect the fact that they are well-differentiated relicts of old, genus-poor lineages. We refrain from suggesting systematic revisions for shorebirds at this time because gene trees may fail to recover the species tree when long branches are connected to deep, shorter branches, as is the case for some of the enigmatic taxa.

Sequences from 14 mitochondrial genes provide a well-supported phylogeny of the Charadriiform birds congruent with the nuclear RAG-1 tree.

Because of the difficulties of constructing a robust phylogeny for Charadriiform birds using morphological characters, recent studies have turned to DNA sequences to resolve the systematic uncertainties of family-level relationships in this group. However, trees constructed using nuclear genes or the mitochondrial Cytochrome b gene suggest deep-level relationships of shorebirds that differ from previous studies based on morphology or DNA-DNA hybridization distances. To test phylogenetic hypotheses based on nuclear genes (RAG-1, myoglobin intron-2) and single mitochondrial genes (Cytochrome b), approximately 13,000 bp of mitochondrial sequence was collected for one exemplar species of 17 families of Charadriiformes plus potential outgroups. Maximum likelihood and Bayesian analyses show that trees constructed from long mitochondrial sequences are congruent with the nuclear gene topologies [Chardrii (Lari, Scolopaci)]. Unlike short mitochondrial sequences (such as Cytochrome b alone), longer sequences yield a well-supported phylogeny for shorebirds across various taxonomic levels. Examination of substitution patterns among mitochondrial genes reveals specific genes (especially ND5, ND4, ND2, and COI) that are better suited for phylogenetic analyses among shorebird families because of their relatively homogeneous nucleotide composition among lineages, slower accumulation of substitutions at third codon positions, and phylogenetic utility in both closely and distantly related lineages. For systematic studies of birds in which family and generic levels are examined simultaneously, we recommend the use of both nuclear and mitochondrial sequences as the best strategy to recover relationships that most likely reflect the phylogenetic history of these lineages.

RAG-1 sequences resolve phylogenetic relationships within Charadriiform birds.

The Charadriiformes is a large and diverse order of shorebirds currently classified into 19 families, including morphologically aberrant forms that are of uncertain phylogenetic placement within non-passerine birds in general. Recent attempts using morphological characters have failed to recover a well-supported phylogeny depicting higher level relationships within Charadriiformes and the limits to the order, primarily because of inconsistency and homoplasy in these data. Moreover, these trees are incongruent with the relationships presented in the DNA hybridization tapestry of, including the location of the root and the branching order of major clades within the shorebirds. To help clarify this systematic confusion we therefore sequenced the large RAG-1 nuclear exon (2850 bp) from 36 species representing 17 families of shorebirds for which DNA was available. Trees built with maximum parsimony, maximum likelihood or Bayesian methods are topologically identical and fully resolved, with high support at basal nodes. This further attests to the phylogenetic utility of the RAG-1 sequences at higher taxonomic levels within birds. The RAG-1 tree is topologically similar to the DNA hybridization tree in depicting three major subordinal clades of shorebirds, the Charadrii (thick-knees, sheathbills, plovers, oystercatchers, and allies), Scolopaci (sandpipers and jacanas) and the Lari (coursers, pratincoles, gulls, terns, skimmers, and skuas). However, the basal split in the RAG-1 tree is between Charadrii and (Scolopaci+Lari), whereas in the DNA hybridization tree Scolopaci is the sister group to the (Charadrii+Lari). Thus in both of these DNA-based trees the Alcidae (auks, murres, and allies) are not basal among shorebirds as hypothesized in morphological trees, but instead are placed as a tip clade within Lari. The enigmatic buttonquails (Turnicidae), variously hypothesized as being allied to either the Galliformes, Gruiformes, or Charadriiformes, are shown to be a basal lineage in the more conventional Lari clade. Divergence times estimated with rate-smoothing methods and minimum time constraints imposed at nodes with key fossils suggest that Charadriiformes originated in Gondwanaland.