Protein-coding Loci Research Articles

A novel human protein-coding locus identified using a targeted RNA enrichment technique.

Accurate and comprehensive genomic annotation, including the full list of protein-coding genes, is vital for understanding the molecular mechanisms of human biology. We have previously shown that the genome contains a multitude of yet hidden functional exons and transcripts, some of which might represent novel mRNAs. These results resonate with those from other groups and strongly argue that two decades after the completion of the first draft of the human genome sequence, the current annotation of human genes and transcripts remains far from being complete. Using a targeted RNA enrichment technique, we showed that one of the novel functional exons previously discovered by us and currently annotated as part of a long non-coding RNA, is actually a part of a novel protein-coding gene, InSETG-4, which encodes a novel human protein with no known homologs or motifs. We found that InSETG-4 is induced by various DNA-damaging agents across multiple cell types and therefore might represent a novel component of DNA damage response. Despite its low abundance in bulk cell populations, InSETG-4 exhibited expression restricted to a small fraction of cells, as demonstrated by the amplification-based single-molecule fluorescence in situ hybridization (asmFISH) analysis. This study argues that yet undiscovered human protein-coding genes exist and provides an example of how targeted RNA enrichment techniques can help to fill this major gap in our knowledge of the information encoded in the human genome.

A fully phased, chromosome-scale genome of sugar beet line FC309 enables the discovery of Fusarium yellows resistance QTL.

Sugar beet (Beta vulgaris L.) is a global source for table sugar and animal fodder. Here we report a highly contiguous, haplotype phased genome assembly and annotation for sugar beet line FC309. Both assembled haplomes for FC309 represent the largest and most contiguous assembled beet genomes reported to date, as well as gene annotations sets that capture over 1500 additional protein-coding loci compared to prior beet genome annotations. These new genomic resources were used to identify novel quantitative trait loci (QTL) for Fusarium yellows resistance from the FC309 genetic background using an F2 mapping-by-sequencing approach. The highest QTL signals were detected on Chromosome 3, spanning approximately 10Mbp in both haplomes. A parallel transcriptome profiling experiment identified candidate genes within the Chromosome 3 QTL with plausible roles in disease response, including NBS-LRR genes with expression trends supporting a role in resistance. Investigation of genetic variants in these candidate genes found one major disease resistance protein containing high effect variants of interest. Collectively, the genomic resources for FC309 presented here be foundational tool for comparative genomics, mapping other traits in the FC309 background, and as a reference genome for other beet studies due to its contiguity, completeness, and high-quality gene annotations.

Rapid Evolution of Host Repertoire and Geographic Range in a Young and Diverse Genus of Montane Butterflies.

Evolutionary changes in geographic distribution and larval host plants may promote the rapid diversification of montane insects, but this scenario has been rarely investigated. We studied rapid radiation of the butterfly genus Colias, which has diversified in mountain ecosystems in Eurasia, Africa, and the Americas. Based on a dataset of 150 nuclear protein-coding genetic loci and mitochondrial genomes, we constructed a time-calibrated phylogenetic tree of Colias species with broad taxon sampling. We then inferred their ancestral geographic ranges, historical diversification rates, and the evolution of host use. We found that the most recent common ancestor of Colias was likely geographically widespread and originated ~3.5 Ma. The group subsequently diversified in different regions across the world, often in tandem with geographic expansion events. No aspect of elevation was found to have a direct effect on diversification. The genus underwent a burst of diversification soon after the divergence of the Neotropical lineage, followed by an exponential decline in diversification rate toward the present. The ancestral host repertoire included the legume genera Astragalus and Trifolium but later expanded to include a wide range of Fabaceae genera and plants in more distantly related families, punctuated with periods of host range expansion and contraction. We suggest that the widespread distribution of the ancestor of all extant Colias lineages set the stage for diversification by isolation of populations that locally adapted to the various different environments they encountered, including different host plants. In this scenario, elevation is not the main driver but might have accelerated diversification by isolating populations.

Back together: Over 1000 single-copy nuclear loci and reproductive features support the holoendoparasitic Apodanthaceae and Rafflesiaceae as sister lineages in the order Malpighiales

The systematics of the holoendoparasitic flowering plant families Apodanthaceae and Rafflesiaceae has been discussed for over two centuries. The morphological reduction of roots, shoots and leaves in all members of both families, resulting in a cryptic mycelium-like vegetative body, has been interpreted either as a key common feature, or as a result of convergent evolution due to full dependence upon their hosts. Historically, the two families have been placed together due to similar morphological features, but recent analyses based on few mitochondrial and ribosomal gene markers placed them in the distantly related orders Cucurbitales and Malpighiales. Here we reevaluate the affinities of the Apodanthaceae and the Rafflesiaceae using a phylogenomic approach. We present (1) a historical account on their affinities over the last 200 years; (2) phylogenetic analyses reinstating their sister group relationship as part of the order Malpighiales, based on over 1000 single-copy nuclear protein-coding loci; and (3) a comprehensive list of putative morphoanatomical and developmental synapomorphies in light of the phylogenomic results, with emphasis on shared reproductive traits regardless of dramatic differences in floral size.

Phylogenomics of the major lineages of Bembidion and related ground beetles (Coleoptera: Carabidae: Bembidiini)

Abstract Bembidion Latreille (Coleoptera: Carabidae) is a genus of small ground beetles containing about 1,380 species. To test previous phylogenetic hypotheses about deeper lineages of Bembidion and near relatives, we recover and examine over 1,800 nuclear protein-coding loci from 33 species representing the main lineages of Bembidion, 10 species of other bembidiine genera, and 7 outgroups. We find that Bembidion exclusive of subgenus Phyla Motschulsky is monophyletic, and we reclassify Phyla as a separate genus. Within Bembidion we find 2 dominant clades, the Bembidion superseries (containing about 490 species in the subgenera Eupetedromus Netolitzky and Lindrochthus Maddison, the Philochthus Stephens complex, and the Bembidion series), and the Ocydromus Clairville superseries (containing almost all other Bembidion representing about 840 species). The only known lineages within Bembidion outside of these superseries are subgenus Hoquedela Müller-Motzfeld and the Desarmatocillenus Netolitzky complex, which combined contain less than 30 species. Most clades are insensitive to variations in analyses and hold up under different sets of taxa and loci, analyses at the nucleotide or amino acid levels, and different analytical methods (maximum likelihood, including posterior mean site frequency analyses, Bayesian analyses, invariant-based methods, and those that consider incomplete lineage sorting). Despite the clarity achieved in most aspects of the phylogeny, there are several unresolved regions, notably the relationships of Desarmatocillenus, Hoquedela, and Phyla to other bembidiines. A divergence dating analysis suggests that crown Bembidion is about 48 million years old (95% confidence intervals 40–58 Ma), and that the 2 large superseries are about 38 million years old (95% confidence intervals about 29–47 Ma).

Analysis of Codon Usage Bias in the Plastid Genome of Diplandrorchis sinica (Orchidaceae).

In order to understand the bias and main affecting factors of codon usage in the plastid genome of Diplandrorchis sinica, which is a rare and endangered plant species in the Orchidaceae family, the complete plastid genome sequence of D. sinica was downloaded from the GenBank database and 20 protein-coding sequences that met the analysis requirements were finally selected. The GC content, length of the amino acid (Laa), relative synonymous codon usage (RSCU), and effective number of codon (ENC) of each gene and codon were calculated using the CodonW and EMBOSS online programs. Neutral plot analysis, ENC-plot analysis, PR2-plot analysis, and correspondence analysis were performed using Origin Pro 2024 software, and correlation analysis between various indicators was performed using SPSS 23.0 software. The results showed that the third base of the codon in the plastid genome of D. sinica was rich in A and T, with a GC3 content of 27%, which was lower than that of GC1 (45%) and GC2 (39%). The ENC value ranged from 35 to 57, with an average of 47. The codon usage bias was relatively low, and there was a significant positive correlation between ENC and GC3. There were a total of 32 codons with RSCU values greater than 1, of which 30 ended with either A or U. There were a total of nine optimal codons identified, namely, UCU, UCC, UCA, GCA, UUG, AUA, CGU, CGA, and GGU. This study indicated that the dominant factor affecting codon usage bias in the plastid genome of D. sinica was natural selection pressure, while the impact of base mutations was limited. The codon usage patterns were not closely related to gene types, and the distribution of photosynthetic system genes and ribosomal protein-coding gene loci was relatively scattered, indicating significant differences in the usage patterns of these gene codons. In addition, the codon usage patterns may not be related to whether the plant is a photosynthetic autotrophic or heterotrophic nutritional type. The results of this study could provide scientific references for the genomic evolution and phylogenetic research of plant species in the family Orchidaceae.

Role of Epigenetic Factors in Determining the Biological Behavior and Prognosis of Hepatocellular Carcinoma.

The current study's objective is to evaluate the molecular genetic mechanisms influencing the biological behavior of hepatocellular carcinoma (HCC) by analyzing the transcriptomic and epigenetic signatures of the tumors. Transcriptomic data were downloaded from the NCBI GEO database. We investigated the expression differences between the GSE46444 (48 cirrhotic tissues versus 88 HCC tissues) and GSE63898 (168 cirrhotic tissues versus 228 HCC tissues) data sets using GEO2R. Differentially expressed genes were evaluated using GO and KEGG metabolic pathway analysis websites. Whole genome bisulfite sequencing (WGBS) and Methylated DNA Immunoprecipitation Sequencing (MeDIP-Seq) data sets (26 HCC tissues versus 26 adjacent non-tumoral tissues) were also downloaded from the NCBI SRA database. These data sets were analyzed using Bismark and QSEA, respectively. The methylation differences between the groups were assessed using functional enrichment analysis. In the GSE46444 data set, 80 genes were upregulated, and 315 genes were downregulated in the tumor tissue (HCC tissue) compared to the non-tumor cirrhotic tissue. In the GSE63898 data set, 1261 genes were upregulated, and 458 genes were downregulated in the cirrhotic tissue compared to the tumor tissues. WGBS revealed that 20 protein-coding loci were hypermethylated. while the hypomethylated regions were non-protein-coding. The methylated residues of the tumor tissue, non-tumorous cirrhotic tissue, and healthy tissue were comparable. MeDIP-Seq, conducted on tumoral and non-tumoral tissues, identified hypermethylated or hypomethylated areas as protein-coding regions. The functional enrichment analysis indicated that these genes were related to pathways including peroxisome, focal adhesion, mTOR, RAP1, Phospholipase D, Ras, and PI3K/AKT signal transduction. The investigation of transcriptomic and epigenetic mechanisms identified several genes significant in the biological behavior of HCC. These genes present potential targets for the development of targeted therapy.

New species of redbait from the Philippines (Teleostei, Emmelichthyidae, Emmelichthys).

We describe a new species of redbait in the genus Emmelichthys collected from fish markets on Panay and Cebu islands in the Visayas region of the Philippines. The species is externally similar to E.struhsakeri but is diagnosable by two prominent fleshy papillae associated with the cleithrum and fewer pectoral-fin rays (18-19 vs. 19-21) and gill rakers (30-33 vs. 34-41). Additionally, mitochondrial DNA differentiates this taxon from other species of Emmelichthys. We generate mitochondrial genomes for two of the three type specimens and several other emmelichthyids to place the new taxon in a phylogenetic context. Analysis of the protein-coding mitochondrial loci calls into question the monophyly of two emmelichthyid genera (Emmelichthys and Erythrocles) and highlights the need for subsequent analyses targeting the intrarelationships of the Emmelichthyidae.

Reassessment of the Phylogenetics of Two Pygmy Grasshopper Generic Groups Tetrix and Systolederus through Mitochondrial Phylogenomics Using Four New Mitochondrial Genome Assemblies.

Mitochondrial genomes offer pragmatic genetic markers to reconstruct evolutionary relationships and inform taxonomic classifications. Here, we present complete mitochondrial sequences for four Chinese pygmy grasshoppers (Tetrigidae), aiming to reevaluate phylogenetic patterns and morphological taxonomy. Our 17,643 bp, 16,274 bp, 15,086 bp, and 15,398 bp mitogenomes of Exothotettix guangxiensis, Formosatettix longwangshanensis, Euparatettix sinufemoralis and Systolederus zhengi, respectively, exhibit archetypal Tetrigidae architecture. We constructed phylogenies using 13 protein-coding loci from 39 Tetrigidae mitogenomes, revealing several genus-level clusters with statistically solid support, conflicts regarding Ex. guangxiensis, F. longwangshanensis merging into Tetrix, and two subclades of Systolederus. The dated divergence analysis indicates over 150 Mya of Tetrigidae ancestry, tracing the Systolederus generic group splits up to ~75 million years ago. Moreover, the Tetrix generic group radiated over 14 Mya across vast distributions, consistent with rapid adaptive dispersals. Our mitochondrial reconstructions suggest that Synstolederus is taxonomically overextended for a single genus, while the distinctiveness of Ex. guangxiensis and F. longwangshanensis from Tetrix appears questionable, and the Tetrix generic group comprises a potential tRNA-Ile coding region. Our integrative mitogenomic approaches will help resolve issues stemming from morphological taxonomy that is reliant on traits that are prone to convergence. This investigation enhances comprehension of Tetrigidae phylogeny and accentuates molecular systematics.

Taxonomy and pathogenicity of fungi associated with oak decline in northern and central Zagros forests of Iran with emphasis on coelomycetous species.

Oak decline is a complex disorder that seriously threatens the survival of Zagros forests. In an extensive study on taxonomy and pathology of fungi associated with oak decline in the central and northern part of Zagros forests, 462 fungal isolates were obtained from oak trees showing canker, gummosis, dieback, defoliation, and partial or total death symptoms. Based on inter-simple sequence repeat (ISSR) fingerprinting patterns, morphological characteristics, and sequences of ribosomal DNA (28S rDNA and ITS) and protein coding loci (acl1, act1, caM, tef-1α, rpb1, rpb2, and tub2), 24 fungal species corresponding to 19 genera were characterized. Forty percent of the isolates were placed in eight coelomycetous species from seven genera, namely, Alloeutypa, Botryosphaeria, Cytospora, Didymella, Gnomoniopsis, Kalmusia, and Neoscytalidium. Of these, four species are new to science, which are introduced here as taxonomic novelties: Alloeutypa iranensis sp. nov., Cytospora hedjaroudei sp. nov., Cytospora zagrosensis sp. nov., and Gnomoniopsis quercicola sp. nov. According to pathogenicity trials on leaves and stems of 2-year-old Persian oak (Quercus brantii) seedlings, Alternaria spp. (A. alternata, A. atra, and A. contlous), Chaetomium globosum, and Parachaetomium perlucidum were recognized as nonpathogenic. All coelomycetous species were determined as pathogenic in both pathogenicity trials on leaves and seedling stems, of which Gnomoniopsis quercicola sp. nov., Botryosphaeria dothidea, and Neoscytalidium dimidiatum were recognized as the most virulent species followed by Biscogniauxia rosacearum.

Isoform-resolved transcriptome of the human preimplantation embryo

Human preimplantation development involves extensive remodeling of RNA expression and splicing. However, its transcriptome has been compiled using short-read sequencing data, which fails to capture most full-length mRNAs. Here, we generate an isoform-resolved transcriptome of early human development by performing long- and short-read RNA sequencing on 73 embryos spanning the zygote to blastocyst stages. We identify 110,212 unannotated isoforms transcribed from known genes, including highly conserved protein-coding loci and key developmental regulators. We further identify 17,964 isoforms from 5,239 unannotated genes, which are largely non-coding, primate-specific, and highly associated with transposable elements. These isoforms are widely supported by the integration of published multi-omics datasets, including single-cell 8CLC and blastoid studies. Alternative splicing and gene co-expression network analyses further reveal that embryonic genome activation is associated with splicing disruption and transient upregulation of gene modules. Together, these findings show that the human embryo transcriptome is far more complex than currently known, and will act as a valuable resource to empower future studies exploring development.

Construction of a High-Density Paulownia Genetic Map and QTL Mapping of Important Phenotypic Traits Based on Genome Assembly and Whole-Genome Resequencing

Quantitative trait locus (QTL) mapping based on a genetic map is a very effective method of marker-assisted selection in breeding, and whole-genome resequencing is one of the useful methods to obtain high-density genetic maps. In this study, the hybrid assembly of Illumina, PacBio, and chromatin interaction mapping data was used to construct high-quality chromosomal genome sequences of Paulownia fortunei, with a size of 476.82 Mb, a heterozygosity of 0.52%, and a contig and scaffold N50s of 7.81 Mb and 21.81 Mb, respectively. Twenty scaffolds with a total length of 437.72 Mb were assembled into 20 pseudochromosomes. Repeat sequences with a total length of 243.96 Mb accounted for 51.16% of the entire genome. In all, 26,903 protein-coding gene loci were identified, and 26,008 (96.67%) genes had conserved functional motifs. Further comparative genomics analysis preliminarily showed that the split of P. fortunei with Tectona grandis likely occurred 38.8 (33.3–45.1) million years ago. Whole-genome resequencing was used to construct a merged genetic map of 20 linkage groups, with 2993 bin markers (3,312,780 SNPs), a total length of 1675.14 cm, and an average marker interval of 0.56 cm. In total, 73 QTLs for important phenotypic traits were identified (19 major QTLs with phenotypic variation explained ≥ 10%), including 10 for the diameter at breast height, 7 for the main trunk height, and 56 for branch-related traits. These results not only enrich P. fortunei genomic data but also form a solid foundation for fine QTL mapping and key marker/gene mining of Paulownia, which is of great significance for the directed genetic improvement of these species.

UPrimer: A Clade-Specific Primer Design Program Based on Nested-PCR Strategy and Its Applications in Amplicon Capture Phylogenomics.

Amplicon capture is a promising target sequence capture approach for phylogenomic analyses, and the design of clade-specific nuclear protein-coding locus (NPCL) amplification primers is crucial for its successful application. In this study, we developed a primer design program called UPrimer that can quickly design clade-specific NPCL amplification primers based on genome data, without requiring manual intervention. Unlike other available primer design programs, UPrimer uses a nested-PCR strategy that greatly improves the amplification success rate of the designed primers. We examined all available metazoan genome data deposited in NCBI and developed NPCL primer sets for 21 metazoan groups with UPrimer, covering a wide range of taxa, including arthropods, mollusks, cnidarians, echinoderms, and vertebrates. On average, each clade-specific NPCL primer set comprises ∼1,000 NPCLs. PCR amplification tests were performed in 6 metazoan groups, and the developed primers showed a PCR success rate exceeding 95%. Furthermore, we demonstrated a phylogenetic case study in Lepidoptera, showing how NPCL primers can be used for phylogenomic analyses with amplicon capture. Our results indicated that using 100 NPCL probes recovered robust high-level phylogenetic relationships among butterflies, highlighting the utility of the newly designed NPCL primer sets for phylogenetic studies. We anticipate that the automated tool UPrimer and the developed NPCL primer sets for 21 metazoan groups will enable researchers to obtain phylogenomic data more efficiently and cost-effectively and accelerate the resolution of various parts of the Tree of Life.

Discovery of long non-coding RNAs in the liver fluke, Fasciola hepatica.

Long non-coding (lnc)RNAs are a class of eukaryotic RNA that do not code for protein and are linked with transcriptional regulation, amongst a myriad of other functions. Using a custom in silico pipeline we have identified 6,436 putative lncRNA transcripts in the liver fluke parasite, Fasciola hepatica, none of which are conserved with those previously described from Schistosoma mansoni. F. hepatica lncRNAs were distinct from F. hepatica mRNAs in transcript length, coding probability, exon/intron composition, expression patterns, and genome distribution. RNA-Seq and digital droplet PCR measurements demonstrated developmentally regulated expression of lncRNAs between intra-mammalian life stages; a similar proportion of lncRNAs (14.2%) and mRNAs (12.8%) were differentially expressed (p<0.001), supporting a functional role for lncRNAs in F. hepatica life stages. While most lncRNAs (81%) were intergenic, we identified some that overlapped protein coding loci in antisense (13%) or intronic (6%) configurations. We found no unequivocal evidence for correlated developmental expression within positionally correlated lncRNA:mRNA pairs, but global co-expression analysis identified five lncRNA that were inversely co-regulated with 89 mRNAs, including a large number of functionally essential proteases. The presence of micro (mi)RNA binding sites in 3135 lncRNAs indicates the potential for miRNA-based post-transcriptional regulation of lncRNA, and/or their function as competing endogenous (ce)RNAs. The same annotation pipeline identified 24,141 putative lncRNAs in F. gigantica. This first description of lncRNAs in F. hepatica provides an avenue to future functional and comparative genomics studies that will provide a new perspective on a poorly understood aspect of parasite biology.

Interspecific gene flow obscures phylogenetic relationships in an important insect pest species complex

As genomic data proliferates, the prevalence of post-speciation gene flow is making species boundaries and relationships increasingly ambiguous. Although current approaches inferring fully bifurcating phylogenies based on concatenated datasets provide simple and robust answers to many species relationships, they may be inaccurate because the models ignore inter-specific gene flow and incomplete lineage sorting. To examine the potential error resulting from ignoring gene flow, we generated both a RAD-seq and a 500 protein-coding loci highly multiplexed amplicon (HiMAP) dataset for a monophyletic group of 12 species defined as the Bactrocera dorsalis sensu lato clade. With some of the world's worst agricultural pests, the taxonomy of the B. dorsalis s.l. clade is important for trade and quarantines. However, taxonomic confusion confounds resolution due to intra- and interspecific phenotypic variation and convergence, mitochondrial introgression across half of the species, and viable hybrids. We compared the topological convergence of our datasets using concatenated phylogenetic and various multispecies coalescent approaches, some of which account for gene flow. All analyses agreed on species delimitation, but there was incongruence between species relationships. Under concatenation, both datasets suggest identical species relationships with mostly high statistical support. However, multispecies coalescent and multispecies network approaches suggest markedly different hypotheses and detected significant gene flow. We suggest that the network approaches are likely more accurate because gene flow violates the assumptions of the concatenated phylogenetic analyses, but the data-reductive requirements of network approaches resulted in reduced statistical support and could not unambiguously resolve gene flow directions. Our study highlights the importance of testing for gene flow, particularly with phylogenomic datasets, even when concatenated approaches receive high statistical support.

Xylaria furcata reconsidered and nine resembling species

BackgroundXylaria collections from termite nests with dichotomously branched stromata have been identified as X. furcata. However, Léveillé’s original material is no longer available, and the modern interpretation of X. furcata is based on a 1908 collection made by von Höhnel from termite nests at Buitenzorg Botanical Garden in Java. A packet of this von Höhnel material at FH was designated as the neotype by Rogers et al. in 2005.ResultsWe reexamined the neotype from FH and its duplicates from various herbaria and found that three different species were mixed in these specimens. Despite that all of them have dichotomously branched stromata and tiny ascospores, only one fits the 2005 neotypification of X. furcata, where exposed perithecial mounds on the stromatal surface were unambiguously indicated. This portion of material is redesignated as the neotype, while the other two species with immersed perithecia are described as new: X. hoehnelii and X. robustifurcata. The ITS sequence obtained from the neotype helped us designate a specimen with cultures obtained from it as the epitype. From specimens identifiable as X. furcata, we describe four new species: X. brevifurcata, X. furcatula, X. insignifurcata, and X. tenellifurcata. Additionally, we recognize X. furcata var. hirsuta at the species level as X. hirsuta and consider X. scoparia a distinct species rather than a synonym of X. furcata. Molecular phylogenetic analyses based on three protein-coding loci showed that X. furcata and resembling species were grouped into two clusters: the X. furcata cluster with half-exposed to fully exposed perithecial mounds and the X. hoehnelii cluster with largely immersed perithecial mounds.ConclusionTen species are recognized for X. furcata and resembling species, all of which could have been identified as X. furcata in the past. Its diversity has been overlooked primarily due to the small and similar stromata. Several additional species have been confirmed to be related to X. furcata by DNA sequences but are yet to be described due to the lack of mature stromata. While the species diversity of macrotermitine termites is equally high in Africa as in Asia, all of the species are primarily found in Asia, with X. hirsuta as the only exception. This suggests that there may be many more undiscovered species for this fungal group.

Tailored machine learning models for functional RNA detection in genome-wide screens.

The insilico prediction of non-coding and protein-coding genetic loci has received considerable attention in comparative genomics aiming in particular at the identification of properties of nucleotide sequences that are informative of their biological role in the cell. We present here a software framework for the alignment-based training, evaluation and application of machine learning models with user-defined parameters. Instead of focusing on the one-size-fits-all approach of pervasive insilico annotation pipelines, we offer a framework for the structured generation and evaluation of models based on arbitrary features and input data, focusing on stable and explainable results. Furthermore, we showcase the usage of our software package in a full-genome screen of Drosophila melanogaster and evaluate our results against the well-known but much less flexible program RNAz.

Effect of Different Types of Sequence Data on Palaeognath Phylogeny.

Palaeognathae consists of five groups of extant species: flighted tinamous (1) and four flightless groups: kiwi (2), cassowaries and emu (3), rheas (4), and ostriches (5). Molecular studies supported the groupings of extinct moas with tinamous and elephant birds with kiwi as well as ostriches as the group that diverged first among the five groups. However, phylogenetic relationships among the five groups are still controversial. Previous studies showed extensive heterogeneity in estimated gene tree topologies from conserved nonexonic elements, introns, and ultraconserved elements. Using the noncoding loci together with protein-coding loci, this study investigated the factors that affected gene tree estimation error and the relationships among the five groups. Using closely related ostrich rather than distantly related chicken as the outgroup, concatenated and gene tree-based approaches supported rheas as the group that diverged first among groups (1)-(4). Whereas gene tree estimation error increased using loci with low sequence divergence and short length, topological bias in estimated trees occurred using loci with high sequence divergence and/or nucleotide composition bias and heterogeneity, which more occurred in trees estimated from coding loci than noncoding loci. Regarding the relationships of (1)-(4), the site patterns by parsimony criterion appeared less susceptible to the bias than tree construction assuming stationary time-homogeneous model and suggested the clustering of kiwi and cassowaries and emu the most likely with ∼40% support rather than the clustering of kiwi and rheas and that of kiwi and tinamous with 30% support each.

13.3 Exome Sequencing in OCD Reveals a Burden of Rare Damaging Coding Variants

The primary objective of this study was to collect nearly all genetic variants from protein-coding loci within the genome (ie, “coding variants”) across a large cohort of OCD cases and to determine if there is evidence for their contribution to overall risk.

Molecular Characterization of a Novel Relapsing Fever Borrelia Species from the Desert Cottontail (Sylvilagus audubonii) in New Mexico, USA.

The Borrelia genus comprises vector-borne, spirochete bacteria infecting vertebrates worldwide. We characterized a novel relapsing fever Borrelia species from a desert cottontail (Syvilagus audubonii) from New Mexico, US, using an established multilocus sequence analysis approach. Phylogenetic analysis of the flagellin gene (flaB) and four other protein-coding loci (clpX, pepX, recG, rplB) grouped the novel Borrelia species with hard tick relapsing fever borreliae Borrelia lonestari, Borrelia theileri, and Borrelia miyamotoi. The identity of the vectors and other vertebrate hosts, geographic distribution, and zoonotic potential of this novel Borrelia species deserve further investigation.